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Linear  damping  factor  in  CX  in  Equation  [1] 

Fourier  cosine  transform  value  at  co 
Fourier  sine  transform  value  at  a 
Frequency  in  cycles  per  second 
An  increment  of  time 

Mass 

Damped  natural  frequency  of  linear  system 
First  forward  difference  at  n 
Second  forward  difference  at  n 
Time  in  seconds 

Relative  displacement  (Y-Z)  or  response  displacement 
P8eudo*response  velocity 
Pseudo-response  acceleration 

Absolute  displacement  of  the  linear  oscillator  from  its 
equilibrium  position 

Displacement  of  foundation  from  its  equilibrium  position 

Ratio  of  damping  to  critical  linear  damping  or  damping  coefficient 

Undamped  natural  frequency  of  linear  system  in  radians  per  second 
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ABSTRACT 


This  paper  describes  a  computer  program  and  a  mathematical  model 
(developed  by  G.J.  O’Hara  of  the  Naval  Research  Laboratory)  for  computing 
shock  and  residual  shock  spectra  of  linear  elastic  structures  responding  to 
transient  forces  and  foundation  motions.  A  wide  range  of  output  options, 
from  tabulated  results  to  a  four-coordinate  log-log  grid,  is  possible. 

Included  in  the  mathematical  analysis  of  the  shock  spectra  is  a  refinement 
in  determining  maximum  responses,  suggested  by  O’Hara,  resulting  in  a 
more  accurate  shock  spectrum.  The  program  is  written  in  FORTRAN  IV 
for  the  IBM  7090  and  includes  plotting  subroutines  for  the  General 
Dynamics  S-C  4020  printer-plotter.  The  program  is  also  coded  without  the 
S-C  4020  subroutines. 


ADMINISTRATIVE  INFORMATION 

The  Naval  Ship  Systems  Command  (formerly  Bureau  of  Ships)  authorized  and  sponsored 
this  work  by  Bureau  of  Ships  letter  10460  Serial  732-7  of  5  Feb  1965. 

INTRODUCTION 

The  Naval  Ship  Systems  Command  requested  the  Applied  Mathematics  Laboratory  of 
the  David  Taylor  Model  Basin  to  develop  a  series  of  programs  widely  applicable  in  scope  for 
static  and  dynamic  structural  analysis.  The  code  presented  is  a  computerization  of  O’Hara’s1 
techniques  for  shock  analysis.  The  program  generates  various  forms  of  shock  and  residual 
shock  spectra  for  the  response  of  a  structure  to  a  transient  force  or  foundation  motion  when 
the  mathematical  n\odel  used  for  the  structure  is  a  single-degree-of-freedom  linear  oscillator. 
The  spectra  are  used  to  find  the  displacement  of  a  structure  caused  by  an  arbitrarily  pre¬ 
scribed  transient  motion  to  its  base,  such  as  response  of  buried  structures  to  an  underground 
blast  or  response  of  structures  within  a  vehicle  during  ballistic  impact.  In  addition,  if  the 
input  shock  function  to  the  program  is  an  acceleration,  the  residual  shock  spectrum  that  is 
generated  is  also  the  Fourier  spectrum  of  the  system. 

The  program  is  written  in  FORTRAN  IV  and  is  executed  on  an  IBM  7090  using  the 
IBSYS,  Version  13,  Monitor  System  in  conjunction  with  a  Stromberg-Carlson  4020  printer- 
plotter.  In  the  appendix,  coding  is  furnished  without  reference  to  the  plotter,  thus  making 
the  program  more  adaptable  to  other  FORTRAN  IV  compilers. 


t 


References  are  listed  on  page  117. 
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MATHEMATICAL  ANALYSIS 


The  mathematical  model  for  shock  analysis  here  is  the  simple  linear  oscillator.  The 
general  differential  equation  for  the  motion  in  Figure  1  is 


MY  +  CX  +  KX  =  F(t) 


[1] 


where  X  is  the  relative  displacement  Y-Z,  CX  is  the  damping  term,  and  KX  is  the  restoring 
force,  based  on  Hooke’s  law.  Now  using  standard  notation,  where  C  is  less  than  the  critical 


Figure  1  -  Linear  Oscillator 


1 


damping  factor,  we  let 


K 


*>*  «  77  ,  «  =  -7—  ,  and  p  =  w 

M  2  Mcj 


Substituting  this  into  Equation  [1],  we  obtain  the  following  equation 


Y  +  2  aco  Y  +  co2  Y  =  2  aco  Z  +  co2  Z  + 


M 


[2] 


There  are  two  types  of  motion  that  are  of  interest  that  result  from,  either  no  applied  force, 
i.e.,  F(t)  =  0  and  foundation  motion,  or  no  foundation  motion,  i.e.,  Z  =  0  and  an  applied 
force.  The  remainder  of  this  section  will  mainly  describe  the  first  case,  where  F(t)  =  0. 
The  latter  case,  where  Z  -  0,  will  be  dealt  with  only  at  the  end.  The  word  shock  will  be 
used  interchangeably  with  the  term  input  function,  namely,  Z,  or  F(t).  When  F(t)  «=  0, 
Equation  [2]  becomes 
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Y  +  2a  mY  +  g>2Y  =  2a  coZ  +  a2  Z, 
or  since  Y  -  Z  =  X, 

X+2ao)X  +  <o2X  =  -Z 


[3] 


Using  Laplace  transforms,  we  find  the  solution  to  Equation  [3]  to  be 


X(t)  =  X (0)  e~aoit  (cos  pt  +  —,~2--=z.  ain  pA  +  X (0) 

V  fl-a2  ) 


e~aco  ‘  sin  pt 


t 

-U 


Z{T)  e-aco^-T^  sin  p(t-T)dT 


[4] 


whereas  the  derivative  of  Equation  [4]  is 


X(t)  ■  e~  a0i  1  X(0)  ^cos  pt  +  ^  a._  sin  pt p  sin  pt 


p  a  cos  pt  "I  ^(0)e“a&)t 

+  J+  ; 


o)  sin  pt  +  p  cos  pt) 


+ 


Z(T)  e 


ao)[t-  r]  j  a(o  sj n  p  (t  -  T)  -  p  cos  p  (t  -  T)]  dT 

[5] 


Equation  [4]  gives  the  relative  response  displacement  of  the  mass  to  the  shock  Z(t)  in  terms 

of  given  initial  displacements  and  velocities.  The  solution  of  the  differential  equation  is  not 

useful  yet,  for  the  integral  must  be  evaluated.  Ndw,  if  the  solution  were  known  at  t  =.<x, 

the  values  of  relative  displacement  (jf)  and  relative  velocity  (^)  could  be  considered  initial 

values  at  t j  for  the  solution  at  tx  +  \t,  using  Equations  [4]  and  [5]  again.  A  convenient 

procedure  followed  in  the  numerical  solution  of  the  differential  equation  consists  of  sub* 

dividing  the  time  record  into  equal  increments,  each  of  duration  At  =  h  seconds,  so  the 

method  of  evaluating  X(t)  at  each  point  will  require  the  same  equations.  Of  course,  the 

values  of  X  and  X  at  the  beginning  of  the  time  record  must  be  given  to  use  this  procedure 

for  solving  the  differential  equation,  and  also  the  input  function  2  must  be  digitized  into 

•* 

equal  increments.  A  parabolic  approximation  for  Z(t)  between  0  and  h  is  substituted  in 
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Equations  [4]  and  [5],  then  the  integrals  are  evaluated.  In  addition,  the  approximation  for 

Z(t)  at  any  point  in  the  time  record  will  go  through  its  preceding  point  as  well  as  its  succes- 

•• 

sive  point.  If  the  time  is  subdivided  into  enough  intervals,  A,  no  sudden  peaks  of  Z  should 

•• 

be  missed  by  the  approximating  parabola  for  Z. 

A  useful  change  in  notation  follows:  X(ih)  =  X{,  Z(ih )  =  Zt,  where  Zi  is  the  ( i  +  1)  th 
point  in  the  time  record,  i  =  0,  1,  2,  3,  . . .  ,  k,  and  k  +  1  is  the  number  of  points  in  the  time 
record. 

••  ••  •• 

The  parabola  passing  through  Zn-V  Zn'  Zn  +  i  yields  the  equation 


m  = 


(Zn  +  1  ~2Zn  +  Zn-l) 

2  A’2 


*2  + 


(Zn  +  l-Zn-l) 

2A 


t  +  z . 


[6] 


for  (n  -  1)  A  -  t  -  (n  +  1)  A.  A  few  notational  changes  will  simplify  Equation  [6]:  let 


Zn+l~Zn  and  Sn 


- 1 


s  -  s  . 

n  n  —  1 


'n  +  1 


-2 


so  Equation  [6]  becomes 


S  t 

n 


=  Zn  +  *T~ 


[7] 


Now,  using  Equation  [7]  in  the  integral  in  place  of  Z(t)  in  Equation  [4],  integrating  by  parts 
from  0  to  A,  and  finally  substituting  A  for  t  in  the  other  terms  of  Equation  [4]  lead  to  the  final 
equation  for  the  response  at  each  point  in  the  time  record: 


6>An  +  1  =  oXne 


-  ao)  h 


(cos  ph  +  -  °L_  sin  ph  ] 

v/TTT2  / 


Y  e 
+  Xn  " 


-  a  cj  h 


sin  ph 


JT, 


£ 

-  —  ^1  -  e~  a  h  ^  cos  ph  +  a  ^  sin  pA^  j 


- —  [l  -  i^.(i_e~a<uA  cos  ph)  - 

cj  L  cjA 


(1  -  2a2)  e~a(oh  sin  ph 
ph 


] 
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and  using  the  same  method  in  developing  Equation  [8]  transforms  Equation  [5]  into 


The  response  X  is  scaled  by  <o  in  Equation  [8],  because  vX  is  a  meaningful  quantity,  as 

•  • 

will  be  explained  later.  Consider  a  time  record  of  Z  versus  time;  the  initial  values  must  be 
given  before  Equations  [8]  and  [9]  can  be  used,  i.e.,  Z(20)  and  Z(fQ)  must  be  known.  41so, 

SQ  for  t  -  tQ  +  h  is  Z(tQ  +  h)  -  Z(tQ),  and  S%_  j  where  n  =  1,  becomes  Z(2  h  +  tQ) 

-2  Z(A  +  £0)  +  Z(tQ). 

Of  more  interest  to  the  engineer  is  not  the  response  ( X )  at  each  point  in  the  time 

record  for  a  particular  frequency,  but  the  maximum  or  minimum  response  for  each  frequency 

over  a  given  range  of  frequencies.  Essentially,  for  the  given  range  of  frequencies  one 

obtains  the  maximum  response  as  a  function  of  frequency.  The  maximum  and  minimum 

response  (X  a)  for  each  frequency  is  selected  from  the  values  of  X  <o  for  the  complete  time 

record.  The  refinement  of  the  maximum  and  minimum  (a  more  accurate  estimate)  is  achieved 

•  • 
by  first  noting  that  the  product  of  Jfmax,  the  derivative  of  the  chosen  -Vmax,  and  Zmax_A, 

the  derivative  of  Z  is  positive  if  reaches  its  true  maximum  between  and 

Xm ax  +  /i  and  is  negative  if  the  true  maximum  occurs  between  ^max  and  ^max_A;  this  is 

illustrated  by  the  two  possibilities  that  follow 
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Here,  let  Xmxx  -  X(nh),  Xmxx  +  h  =  *((n  +  1)A),  and  Xmxx_h  =  ^((n  -  1)A) 

Possibility  1 

AT((n  - 1) A)  •  i(nA)  is  negative,  so  true  maximum  occurs  between  Xmxx_h  and  Xmmx 


(n  —  ).)h  nh  (n  +  l)h 

TIME 


Possibility  2 

X((n  -  1)A)  •  X(nh)  is  positive,  so  true  maximum  occurs  between  Xmmx  and  Xmmx  +  h 


(n-l)h  nh  (n  +  l)h 
TIME 
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Equation  [8]  is  used  after  it  is  determined  whether  the  true  maximum  or  minimum  is 

within  the  range  of  or  .  *.  Xais  then  evaluated  at  intervals  of 

0.2  A,  and  the  largest  X  to  in  the  subinterval  is  taken  as  the  more  accurate  (refined)  maximum. 

The  method  for  computing  the  refined  minimum  value  of  Xto  is  the  same.  The  calculation 

for  X  <u  at  a  given  frequency  is  stopped  at  a  time  which  corresponds  to  the  start  of  residual 

time  (when  the  shock  or  impulse  ceases),  which  is  usually  determined  by  the  evnerimenter 

•  • 

from  examining  the  plot  of  the  input  function  Z(t). 

Some  definitions  of  variables  appearing  in  the  mathematical  analysis  are  in  order  at 
this  point. 

X  =  response  displacement 
Xa  =  pseudo-response  velocity 
Xco2  =  pseudo-response  acceleration 
X  =  true  response  velocity 

y/(X o)2  +  X2  =  residual  response  (in  velocity  units)  at  the  end  of  the 
time  record,  i.e.,  at  the  beginning  of  residual  time.  The  residual  response  is  computed  only 
when  there  is  no  damping. 

The  shock  spectrum  is  usually  the  plot  of  maximum  !Zo>|  versus  f,  where  f  is  2  no,  and 
the  residual  shock  spectrum  is  a  plot  of  J(X  co)2  +  X2  versus  f. 

When  o)  =  0,  equations  for  X  co  and  X  as  already  presented  become  meaningless,  even 
though  X  at  u>  =  0  has  physical  meaning.  New  equations  are  derived  from  [8]  and  [9]  by 
using  L’Hospital’s  rule.  In  Equations  [8]  and  [9]  o  is  set  equal  to  zero,  and  when  an  inde- 
terrrinant  form  is  encountered,  L’Hospital’s  rule  is  used;  therefore,  yielding 


X 


n  +  1 


Xn  +  Xnh  - 

n  fi 


Snh2  h2 
6  24 


[10] 


X 


n  +  1 


-  Zh  -  s 

FI  FI 


S„2-lA 


12 


[11] 


It  is  observed  that  damping  does  not  play  any  role  in  the  foregoing  equations. 

An  added  feature  of  the  residual  shock  spectrum  solution  is  that,  if  the  input  shock 
function  is  an  acceleration  (Z),  then  the  Fourier  spectrum  is  identical  to  the  residual  shock 
spectrum.  The  Fourier  spectrum  magnitude  at  a>  is  defined  as  ( F 2  +  F2)i/2t  where  it  is  found 
that 


Fc  «=  -  X  co  sin  a  t  -  X  cos  o>  t 
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and 


Fs  =  A  a  cos  at  -  X  sin  oo  t 

By  substitution,  ( Fe 2  +  Fs2)1/2  becomes  [(A'  a)2  +  A2]1/2,  which  is  the  residual  shock  spec¬ 
trum  magnitude  for  a.  The  Fourier  spectrum  is  not  necessarily  equivalent  to  the  residual 
shock  spectrum  when  the  input  shock  spectrum  is  Z.  A  thorough  discussion  of  Fourier  spec¬ 
tra  can  be  found  in  Reference  1.  The  author  is  presently  finishing  a  program  dealing  with 
Fourier  spectra  such  as  computing  the  Fourier  spectrum  magnitude,  the  phase  angle,  the 
Fourier  cosine  transform,  and  the  Fourier  sine  transform,  so  a  detailed  description  of  Fourier 
spectra  will  be  contained  in  a  subsequent  report  from  the  Applied  Mathematics  Laboratory  at 
the  David  Taylor  Model  Basin. 

•• 

Until  now  the  input  shock  function  was  assumed  to  be  Z  (foundation  acceleration),  so 
the  necessary  equations  for  an  input  function  which  is  the  foundation  velocity  (Z)  will  now 
be  considered. 


5n  =  *„+1 


S2  ,  =  sn  -Sn  .  =  Z„  ,  .  -2Z„  +  Z„  , 

n  —  1  n  n  —  1  n  +  1  n  n  - 1 


Z(t)  ~zn  +  sn 


t 

h 


[12] 


•  •  • 

where  Z  is  replaced  by  Z. 


Xn  +  i  w  =  Xn  e 


acti  h 


Xn  e  ~  aa>  h  sin  ph 

+  J=72 


e-a«ih 


2a 

Oih 


+  e 


-  a  co  h 


[(? 


/I  -  2a2 

V  Oih 


[13] 
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Equation  [13]  is  obtained  from  Equation  [4]  by  first  integrating  by  parts,  secondly,  substituting 
[12]  for  Z(t)  in  the  new  integral,  integrating  again  by  parts,  and  finally  substituting  A  for  t 
in  the  other  terms  of  Equation  [4].  Equation  [14]  is  obtained  in  a  similar  manner,  using 
Equation  [5]  in  place  of  Equation  [4].  When  eo  equals  0, 


4 


and 


X 


n  +  1 


The  rest  of  the  calculations  using  Z  are  the  same  as  those  using  Z,  except  that  for 
the  Z  input,  the  residual  shock  spectrum  is  not  necessarily  identical  to  the  Fourier  spectrum. 

Finally,  to  solve  a  problem  where  there  is  no  foundation  motion  but  an  applied  force 
F(t),  we  must  solve  the  differential  equation. 


X  +  2  a  (oX  +  a2  X  = 


m 

M 


[15] 


Equations  [15]  and  [3]  are  the  same  except  for  their  right-hand  members,  so  if  the  input  func* 

..  Fit) 

tion  is  an  applied  force,  viz.,  F(t),  then,  since  Z(t)  =  -  — —  ,  all  of  the  F' s  are  scaled  by 

n  * 

1 

-  —  and  the  input  function  is  treated  as  if  it  were  Z(t).  The  output  would  correctly  be  the 
response  to  F(t). 
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OUTLINE  OF  THE  PROGRAM 


The  Roman  numeral  headings  below  correspond  to  the  comment  cards  in  the  program 
listing.  Tables  1  and  2  and  the  flow  charts  in  the  following  section  should  aid  in  under¬ 
standing  the  ensuing  outline.  In  addition,  0  is  used  to  denote  an  alphabetic  0  in  FORTRAN. 
The  Roman  numeral  above  any  figure  in  the  flow  chart  corresponds  to  a  numeral  in  the  fol¬ 
lowing  description. 

I.  Input 

•  *0 

Data  is  read,  such  as  the  forcing  function  (Z  or  Z),  output  options,  and  other  constants. 

IPS  (four-coordinate  grid  option)  is  checked  for  zero,  if  not  zero,  IP  1  and  IP2  (other 
plots)  are  deleted  even  if  they  are  requested. 

II.  Write  and  Plot  Input  Data 

The  time  T  corresponding  to  each  Z  is  generated  by  use  of  the  time  interval  H. 

Each  Z  and  corresponding  T  are  written.  From  statements  110  to  123,  Z  versus  T  is 
plotted. 

III.  Compute  S(N)  and  S2ND(N)  for  Each  Time  T 

Each  S  and  S2ND  are  computed. 

A  frequency  card  is  read,  i.e.,  FREQl,  FREQ2 ,  DELTAF,  then  if  FREQ1  -  FREQ2, 
the  blank  card  was  read,  indicating  the  run  is  finished.  If  FREQl  t  FREQ2  the  program 
continues. 

The  number  of  frequencies  is  computed  from  the  card  just  read,  i.e.,  N0FREQ  - 
(FREQ2  -  FREQl)/ DELTAF. 

IV.  Damped  or  Undamped;  if  Undamped,  Find  Number  of  Residual  Points 

The  large  D0-loop ,  which  is  indexed  for  each  alpha,  is  executed. 

I  ALPHA  is  set  equal  to  one  if  no  damping  is  present  ( a  ■  0),  or  set  equal  to  two  if 
there  is  damping  (a  £  0).  If  no  damping  is  present,  the  start  of  residual  time  is  computed, 
viz.,  T(KH). 

V.  Compute  Constants 

Many  mathematical  expressions  used  repeatedly  in  the  program  are  computed. 

IFREQ  is  set  equal  to  one  if  the  frequency  is  zero  or  is  set  equal  to  two  if  the 
frequency  is  not  zero. 

Finally,  there  is  a  major  branch  in  the  program,  depending  on  whether  (1)  Z  *>  Z  and 
frequency  -  0,  (2)  Z  -  Z  and  frequency  +  0,  (3)  Z  -  Z  and  frequency  =  0,  and  (4)  Z  -  2  and 
frequency  /  0.  Of  the  four  possibilities  just  mentioned  only  possibility  2  will  be  assumed 
in  the  rest  of  the  discussion. 
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VI.  Compute  X0MEGA  (N),  XD0T  For  Velocity  Input 

More  constants  and  mathematical  expressions  are  evaluated.  N  XOMEGA's  and 
X D0T's  are  computed  in  a  D0-loop  ending  at  statement  210. 

Then  there  is  transfer  to  preliminary  calculations. 

VII.  Preliminary  Calculations  for  Minimum  and  Maximum  Response 

Maximum  and  minimum  X0MEGA' s  are  determined;  then  control  is  transfered  to 
statement  225. 

Statements  225  to  235  contain  the  refinement  calculations  for  maximum  and  minimum 

X  co. 

Then  from  statements  211  to  215  the  different  output  features,  which  are  selected  in 
the  input,  are  executed. 

VIII.  Either  Restart  With  New  Frequency  or  Restart  With  New  Alpha  or  Stop 

FREQ  and  FREQ2  are  compared;  if  equal,  the  end  of  the  alpha  D0-loop  is  reached- 
statement  800-and  either  the  entire  calculation  is  restarted  for  a  new  alpha,  or  the  next 
frequency  card  is  read. 

If  FREQ  4  FREQ2,  then  FREQ  is  incremented  by  DELTAF,  and  the  calculations  for 
X0MEGA ,  XD0T,  etc.  are  repeated. 
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TABLE  1 


Glossary  of  Important  FORTRAN  Variables  (Excluding  Input  Variables) 


Variable 

Math  Symbol  and/or  Explanation 

1  Dimensioned 

S 

Sp  first  forward  difference  where  St  -  Zt  +  j  -  Zf.. 

S2ND 

S2,  second  forward  difference  where  S2  ■  Zi+  t  -  2Zf  + 

X0MEGA 

X  co  (and  A’  when  co  ■■  0)  response  displacement  times  natural 
frequency. 

X 

X,  response  displacement,  used  for  plotting  in  /PI,  also  used  for 
plotting  -co2  X  in  /P2. 

XD0T 

X ,  response  velocity. 

0SCIL 

P,  frequency. 

VEL 

Maximum  absolute  value  of  Xc o  for  each  <u. 

RES1D 

y/(X  co)2  +  k2 ,  at  the  start  of  residual  time. 

Undimensioned 

41 

1Qlop,  which  is  the  value  of  the  first  4-line  drawn  on  the  four- 
coordinate  graph  (lower  left-hand  corner). 

Dl 

10L®P,  which  is  the  value  of  the  first  0-line  drawn  on  the  four- 
coordinate  graph  (upper  left-hand  corner). 

IER 

Number  of  points  out  of  range  when  plotting  response  (X)  versus 
time  in  IP  1. 

IERA 

Number  of  points  out  of  range  when  plotting  response  velocity 
( X )  versus  time  in  /Pi. 

1ERB 

Number  of  points  out  of  range  when  plotting  - X  co2  versus  time 
in  IP2. 

IERR 

Number  of  points  out  of  range  when  plotting  the  input  function  Z. 
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FLOW  CHART 
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NOTE:  Aj  11  Initially  *»L 


FLOW  CHART  (Continued) 
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INPUT 


The  input  is  read  into  the  computer  mostly  in  the  form  of  FORTRAN  NAMELISTS; 
consequently,  the  format  for  the  data  is  not  very  restrictive.  The  data  are  usually  on 
punched  cards.  If  the  shock  function  Z  is  available  on  an  analog  tape,  a  conversion  to 
digitized  tape  or  to  cards  must  be  accomplished  before  the  program  can  be  used.  Of  course, 
the  NAMELISTS  can  be  replaced  by  regular  format  statements,  which  would  be  mandatory  if 
the  program  was  converted  to  FORTRAN  II. 

See  Table  2  for  the  list  of  input  variables  with  an  explanation  of  each  variable. 

ZAZA1  to  ZAZAl  are  NAMELIST  names,  each  containing  one  or  more  variables.  (See  any 
IBM  FORTRAN  IV  manual  for  an  explanation  of  NAMELISTS.)  The  data  cards  are  arranged 
as  follows:  SDATA  card  first,  then  JSZAZA1  data,  SZAZA2  data,  SZAZA3  data,  SZAZA1 
data,  frequency  cards,  and  finally  the  system  cards:  END  OF  FILE ,  81 B SYS,  8REST0RE, 
and  END  OF  FILE. 

There  is  no  limit  to  the  number  of  frequency  ranges  that  can  be  handled,  at  one  card 
per  range.  The  last  card  must  be  a  blank  card.  The  frequency  cards  are  read  under  format 
statement  920,  viz.,  3F20.8,  so  the  decimal  point  of  the  first  number  FREQ1  is  in  column  12, 
that  of  the  second  number  FREQ2  in  column  32,  and  that  of  the  last  number  DELTAF  in 
column  52.  FREQl  must  not  equal  FREQ2\  therefore,  the  first  number  on  the  frequency  card 
must  be  less  than  the  second  number.  The  number  of  frequencies,  ( FREQ2  -  FREQl)  -f 
DELTAF,  must  not  exceed  the  storage  capacity  of  the  variable  0SCIL  when  using  IPS,  which 
is  now  dimensioned  for  1000  words. 

If  IP  1  to  /P5  are  set  equal  to  zero,  the  only  output  that  will  result  is  a  plot  of  the 
input  function,  a  tabulation  of  the  input  function,  and  a  tabulation  of  S(  and  S ?. 

If  the  shock  function  Z  is  not  on  cards,  but  generated  from  a  formula  inserted  into  the 
program,  then  the  READ  (5,  ZAZA  1)  statement  must  be  deleted  from  the  program. 

The  physical  units  for  most  of  the  quantities  used  in  this  program  are  assumed 
consistent.  The  units  for  the  input  quantities  determine  the  units  for  the  output  values  ex¬ 
cept  where  specified;  for  example,  the  input  frequencies  should  be  in  cycles  per  second. 
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TABLE  2 


Input  Variables 


FORTRAN 

Name 

Comment 

Z 

Shock  function  which  is  dimensioned,  either  acceleration  or  velocity. 

T 

Initial  time,  really  T (1),  and  is  usually  zero. 

TZER0 

Starting  time  of  residual  shock,  i.e.,  end  of  transient. 

H 

Time  increment  into  which  record  is  subdivided. 

XO 

Response  displacement  at  time  T ,  usually  zero. 

N 

Number  of  points,  i.e.,  number  of  Z’s,  must  not  exceed  the  dimen¬ 
sioned  value  for  Z  in  the  DIMENSION  statement. 

XD0T0 

Response  velocity  at  time  T ,  usually  zero. 

ISZ 

If  ISZ  *=  0,  then  Z  ■»  Z;  if  ISZ  £  0,  then  Z  ■  Z. 

IP1 

Output  Option  1,  if  /Pi  -  0,  ignored;  if  /Pi  ^  0,  plots  of  response 
displacement  versus  time  and  response  velocity  versus  time  are 
generated  for  each  frequency  and  each  alpha. 

IP2 

Output  Option  2,  if  /P2  -  0,  ignored;  if  IP2  A  0,  a  plot  of  negative 
pseudo-response  acceleration  versus  time  is  generated  for  each 
frequency  and  each  alpha. 

IPS 

Output  Option  3,  if  IPS  ■  0,  ignored;  if  IPS  t  0,  a  four-coordinate, 
log-log  graph  is  produced  for  each  range  of  frequencies  submitted, 
as  well  as  for  each  alpha,  consisting  of  response  displacement 
versus  frequency,  pseudo-response  velocity  versus  frequency, 
pseudo-response  acceleration  versus  frequency,  and,  if  no  damping 
(alpha  «  0),  residual  response  velocity  versus  frequency.  If  IPS  is 
used,  i.e.,  IPS  f  0,  the  program  will  inhibit  /PI  and  IP2. 

IPi 

Output  Option  4  is  ignored  if  /P4  -  0;  if  /P4  +  0,  tabulation  of 
results  is  generated  (see  output  for  explanation). 

IPS 

Output  Option  5  is  ignored,  if  IPS  -  0;  if  IPS  A  0,  tabulation  of 
results,  more  extensive  than  /P4,  is  generated  (see  output  for 
explanation). 

\  ALPHA 

Damping  coefficient,  which  is  dimensioned. 

NALPHA 

Number  of  damping  coefficients  submitted  and  must  not  be  greater 
than  the  dimension  of  ALPHA  in  the  DIMENSION  statement. 

FREQ1 

Initial  frequency  in  cycles  per  second. 

FREQ2 

Last  frequency  in  cycles  per  second. 

DEL TAF 

Frequency  increment  in  cycles  per  second. 
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OUTPUT 


The  output  from  the  program  is  composed  of  a  plot  and  tabulation  of  the  input  shock 
function,  a  printout  of  each  Sk  and  Sk,  and  the  input  options.  If  the  printed  Sk  and  S2  are 
not  desired,  the  statement  immediately  following  statement  200  in  the  program  should  be 
removed. 

See  Figures  3  and  8  for  an  example  graph  of  an  input  function  and  Sample  Problems  1 
and  2  for  the  tabulation  of  all  points  plotted  for  the  input  functions.  Sample  Problems  1  and 
2  (Appendixes  A  and  B)  also  include  sample  printouts  of  Sk  and  Sk. 

Figures  9,  10,  12,  and  13  are  examples  of  l Pi,  which  is  a  plot  of  X  versus  time  and 
X  versus  time  for  each  frequency  and  each  alpha  (damping  coefficient). 

IP2,  illustrated  in  Figures  11  and  14,  is  a  plot  of  -o2  X  versus  time  for  each  frequency 
and  each  alpha.  The  quantity  -a>2X  is  equivalent  to  the  absolute  acceleration,  T,  when  Z  =  0. 

IPS,  illustrated  in  Figures  4-6  and  15-18,  is  a  four-coordinate,  log-log  grid  which  is 
a  plot  of  velocity  (ordinate)  versus  frequency  (abscissa),  with  acceleration  lines  having 
negative  slope  and  velocity  lines  having  positive  slope  superimposed  on  an  originally  two- 
coordinate  grid.  The  value  of  the  A  lines  increases  logarithmically  from  left  to  right,  and 
the  line  labeled  A  in  the  lower  left  corner  of  the  grid  corresponds  to  the  value  for  A  printed 
in  the  lower  right  margin.  The  value  of  the  D  lines  decreases  downward  logarithmically,  and 
the  line  labeled  D  in  the  upper  left  corner  of  the  grid  corresponds  to  the  value  for  D  printed 
in  the  lower  right  margin.  The  shock  spectrum  is  the  starred  plot  (X  «  versus  /),  and  the 
residual  shock  spectrum  (\/(X  a)2  +  X2  versus  /),  plotted  only  when  a  -0,  is  the  plot  of 
zero  points.  One  complete  four-coordinate  graph  is  generated  for  each  frequency  range  and 
damping  coefficient  a  .  Hence,  for  any  given  frequency  within  the  plotted  range,  one  could 
read  from  the  starred  curve  the  pseudo-acceleration  A,  the  pseudo-velocity,  and  the  displace¬ 
ment  D;  one  could  also  read  from  the  zero  curve,  provided  there  is  no  damping,  the  residual 
shock  spectrum  magnitude. 

IPi  (Short  Form)  is,  for  each  frequency  range  and  damping  coefficient,  a  tabulation  of 
frequency,  the  maximum  value  between  the  pseudo-velocity  (A  o>)  and  the  residual  response 
(J(X  a>)2  +  X2),  and,  if  no  damping  is  present,  the  residual  response. 

IPS  (Long  Form)  is,  for  each  frequency  range  and  damping  coefficient,  a  tabulation  of 
frequency,  maximum  positive  pseudo-velocity,  minimum  negative  pseudo-velocity,  and,  if  no 
damping  is  present,  the  residual  response.  See  the  sample  problems  for  examples  of  IPi  and 
IPS  (Short  and  Long  Forms). 

One  method  for  utilizing  the  various  graphical  outputs  is  first  to  generate  the  log-log 
grid  (shock  spectrum)  from  IPS,  then  select  from  the  shock  spectrum  the  important  frequencies. 
Next,  these  frequencies  are  used  as  input  in  running  the  program  a  second  time  for  IP  1  and/or 
IP2  to  produce  response  versus  time  plots  for  each  desired  frequency. 
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OPERATING  INSTRUCTIONS 


The  subroutines  required  are  from  the  IBM  7090  plotting  package  written  by  North 
American  Aviation  for  the  Stromberg-Carlson  4020  (charactron  plotter-printer).  If  the 
S-C  4020  is  not  available  to  the  user,  the  alternate  deck  is  used,  in  which  no  S-C  4020 
subroutines  are  called,  and  there  are  only  two  output  options,  /P4  and  IP 5  (see  Appendix  C). 

No  input  tapes,  except  for  unit  5,  are  required;  the  program  is  run  under  IBSYS  version 
13.  Three  output  tapes  are  needed  for  the  program,  one  tape  for  the  S-C  4020,  if  used;  logical 
tape  unit  6;  and  logical  tape  unit  8.  /P4  is  written  on  unit  8,  and  IP5  is'written  on  unit  6. 

An  approximation  of  the  running  time  of  the  program,  using  the  input  setup  from 
Sample  Problem  1  (see  Figure  2),  is  about  7  minutes,  including  compiling  time.  The  running 
time  for  the  program  is  increased  considerably  if  more  than  just  a  few  frequencies  are  used 
with  /PI  and  IP2. 

Most  of  the  memory  is  used  (32tf);  see  the  dimension  statement  of  the  program-listing 
for  the  array  sizes  which  influence  the  amount  of  memory  needed. 
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APPENDIX  A 

SAMPLE  PROBLEM  1 

Arrangement  of  Input  Data 
and 

All  the  Various  Forms  of  Output  That  Can  Be 
Generated  Except  /Pi  and  IP2 
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SAMPLE  PROBLEM  1 


The  input  shock  function  is  a  foundation  acceleration,  Z  =  100  [e“12  9t  sin  200  wt 
+  e~206t  sin  600  nt],  which  is  generated  in  the  program  for  each  desired  point.  The  follow¬ 
ing  cards  are  inserted  at  the  beginning  of  the  program  to  generate  the  Z's. 


1 02  Z(l)=100  *(EXP(— 129.  *T(I))*SIN(TIZ)+EXP(-208.*T(I))*SIN(3.*TIZ))  \ 


I 


TIZ=2  *PI  *100  *T(I) 

11  i  c  u 


00102  1=1,  N 


OOOUOOOOOO  00000  000000000000000000000000000000000000000001)000000000000000000000 
t  :  i  <  ■.  <  •  i  !  >c  "C  * -t  c  *  -s :«  n  .-in*  i- 1:  uw  n»j)  »»<o  <I<I<1<4<1H4I  u«««  !MMi«  i;»rsi  HHmc !.’«««  u  mi  :t  n  t  3  i  :i  »• 

i  ii  i  m  ii  inn  1111111111111111111111111111111111111111111111111111111111111111 

22222223222  22222222222222222222222222222222222222222222222222222222222222222222 
33333333333333  3  3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 33 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 33 3 3 
414(44  4444444444444444444444444444444444444444444444444444444444444444444444444 
55555555555555555  5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 
SS6S6SS  GEESESt8SG666SC6&GC6S6E6SSE6666BSSG6G6S(SSGS6S66SfiCSS6EC66SSS65SSS66ES£S 
77717777777777777777777777772177777777777777777777777777772777777777777777777777 
Minimum  i  iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiitiiiiimiitiiiiiiiiiiiiiiiin 

9390939993099  999999993999399999999999999999993999999999999939999999999999999999 

<  ;  i  <  s  i  i  i  i  n  n  i  i  «  ii  n  i-  a  i  n  n  j:n:<  n.i  r  n.-t jo  :i  a  ii  is*:;  jait<a<iuu<i<!H<?  <i«»r  uumkh  r  uiinu  u  umuhi' min  n  u  m  n  ;mic 


See  Figure  2  for  the  arrangement  of  the  input  data. 
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END  0F  FILE 


r  fRESTORE 


flBSYS 
^  END  0F  FILE 


0.  100. 

o7  1000. 

r  fZAZM  ALPHA  (1)-  0.  NALPHA  -  1| 


^  JZAZA3  IP1  =  0,  IP2=1,  IP3=1,  IP4  =  1,IP5-1» 


H  =  .3000E  —  04,  X0=0,  N  =  1200,  XOOT  0  =  0,  ISZ  =  If 
JZAZA2  T(l)  =  0,  TZERO  =  .03, 


JDATA 


1. 
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(SYSTEM  CARD) 


(SYSTEM  CARD) 


(SYSTEM  CARD) 
(SYSTEM  CARD) 


(BLANK  CARD) 
(FREQUENCY  CARD[ 
(FREQUENCY  CARD) 


(3RD  NAMELIST) 


(2ND  NAMELIST) 


(1ST  NAMELIST) 


(SYSTEM  CARD) 


III.OIIOIKKaOOOliaigOMIMigillllltlOOOMMMIIMMKDtRIIIIIMOIlltMIBSIIMIIIII 

I  I  I  •  »  •  I  •  I  ti  K  tj  I!  M  n  «  (I  •  «  1*  H  ,!,!»  7,  »  II  X!t  M  1!  1, 11  M  14  ml  HU  <•  till  1144  44  H  <1  41  44  41  41  47  »  «  44  49  41  »  MM  II II  9194  IMMIM II  IfftltlHnmiltlltf 

I  I  ,  I  .  II  I  I  I  II  I  I  I  I  I  I  II  II  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  II  I  I  I  I  I  II  I 
77771177777711177777777711777171777717171117711777727777777771777777717777777777 
1 1  1  3  J  11 11 1 J 1 J 1 J J 1 3  3  J  3  1  )  n  n  1  U J  J  J  J I )  J  3  J 1 J J  3  3  J  ))  J  J  3  I J  J  J J  3  3  1 3 ))  I ))  1 3 11 )) )) 1 11 )  11 
4  444444444444**44*4*4*44444444444444*4444444444****444444*44444444*444444444444 
55555555555555555555 S55S555S5555 5555555555555555555 5 555555555555 5555555555555355 

11117111117111117171711117711177711111)71711711111111111111171111777111111711177 

'IIIIIIIIIMIMIIIMMIIIMMIIIIIIMIIIIIMIMIIMIIIIMIMIIIIHIIIIMIIIIIIII 

HmiiiimnninnMMiininnmmimMmnMimiinmiimiiiimmi 

I  I  I  I  9  •  I  #  9  19  H  17  II  M  14  If  17  If  If  79  71  77 11 74  14 14  71  If  74  If  .11  II 11  94  14  14  II  V  J4 19 41 47  4144  44  It  4'  1944  99  91  47  41  44  94  99  91  99  99  4911  If  91 44  ft  49  4199  99  If  It  11 11 14  14  If  II  If  14  tf 
I  ill.  •  ••(Ill 
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Figure  2  -  Arrangement  of  Input  Data  for  Sample  Problem  1 

NOTE:  NAMELIST  0ZAZA1  it  omitted  became  the  Zj’a  are  generated  within  the  program. 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  1 


SNERD  -  EFN  SOURCE  STATEMENT  -  IFN(S) 


06/10 


DIMENSION  2(2000 >.T! 2000). S( 2000). S2N0I 2000) . XOMEGA (2000 ) tX (2000) * 
IXOOT ( 2000 I • ALPHA ( 25  ) •OSCIL( 10C0 t  .VEL( 1000) .RES ID ( 1000) 

EQUIVALENCE  (RESIO.X(n).(VEL.X(lOOI  It 

NAMEL IST/ZAZAI /Z/ZAZA2/T . TZERC .  H.  XO.N* 

1 XOOTO . ISZ/ZAZA3/IPI.IP2.IP3.IP4. I P5/ZAZA4/ ALPHA .N ALPHA 
MR ITE( 6. 900 ) 

P  1*3. 1*15927 
REA0(5.ZAZA2 ) 

REAC(5./AZA3) 

REAOI S, ZAZA4 ) 

IF ( IP3tS0.Sl.S0 

50  IP1*0 
IP2*0 

51  CONTINUE 

— - - WRITE  AND  PLOT  INPUT  OATA - - 

00  100  1*2. N 
100  T( I )*T ( 1-1 )4H 
00  102  1*1. N 

TIZ*2.«PIM00.«T(  I  } 

102  Z ( I )*10C.*(EXP(-1 29  •*  T ( I ) ) *SI N( TI Z )4EXP ( -206. *T ( I ) ) *S IN(3.*TIZ)I 
IF ( ISZI104.10S.104 

104  MR ITE( 6.90 1  ) 

CO  TO  106 

105  WRITE ( 6 . 902 ) 

106  CONTINUE 

MR ITE( 6.903  1  (  Z ( I ).T( I ) • I *1 .N ) 

ZMIN*Z( 1 ) 

ZMAX*Z ( 1 > 

00  110  1*2. N 

IF( ZMAX— Z( I ) tioe. 109.100 
toe  zmax*z(  1 1 
CO  TO  110 

109  IF(ZM|N-Z( I ) 1110.110.107 

107  ZM IN*2 til 

110  CONTINUE 

CALL  CAMRAVOSI 
ZN*N 

CX*H«ZN/IO. 

0V*( ZMAX— ZMIN ) / 10. 

CALL  GRI01V( l.T( 1 )• T ( N ). ZMIN. ZMAX, OX, OV. 1 .1 .1 .1 .6.6) 

CALL  APLOTVIN.T.Z, 1.1, 1,42. IERR ) 

IF ( IERRtllS.116.llS 
US  MR  I  TE(  6.904  I  IERR 
116  CONTINUE 
NN»N— 1 

00  120  1*1 »NN 

CALL  LINEV(NXV(T(  I  }  I.NVV(Z( 1  I l,NXV( TI I  41  I t ,NYV(Z( 141)1 t 

120  CONTINUE 

CALL  PR 1NTV(— IS, 1SHT I ME  IN  SEC0N0S.4S2 .6) 

IF( ISZ)121.122.121 

121  CALL  APRNTV (0 14.-23, 23HFOUNDATION  ACCELERATION. 4.6961 
GO  TO  123 

122  CALL  APRNTV 1 0 ,  —  *4,  — 19. 1 9HF0UNDATI0N  VELOCITY. 4 .664) 

123  CONTINUE 
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C 

C - COMPUTE  SC  N  )  AND  S2NDCN)  FOR  EACH  TIME.T - 

sm*z(2)-zm 

S2NCC 1>*Z(3)-2.«Z<2>*Z< 1) 

K*N— I 

00  200  I  *2  .  K 
S(  I >*Zt  I ♦ 1 >-Z<  I  I 

200  S2N01  I)=Z<  I*1)-2.*Z<  I  »«-Z(  I  — 1  > 

WRITE <6. 91?)  I  I tS( I ) ,S2NO< I >. I*1.N) 

C 

190  REAOI 5*920 )  FREO 1 .FREQ2 .OELTAF 

IF( ABS(  FREQ  1 -FREQ 2 )  — • IE-20)  SOI .601.191 

191  CONTINUE 

IF  I IP3) 192. 195. 192 

192  TEMP*FREQ2-FREQ 1 

IF! AMODI TEMP tOELTAF )-OEL TAF/2 . ) 193 . 1 9* . 1 94 

193  NCFREO*TEMP/OELT  AF 
GO  TO  195 

194  N0FRE0*TEMP/0ELTAF*1.0001 

195  CONTINUE 

C - DAMPEO  OR  UNDAMPED.  IF  UNDAMPED  FIND  NUMBER  OF  RESIDUAL  POINTS 

00  S00  IX*1.NALPHA 
N  IN*  1 
INI  *0 

FR£Q*FREQ 1 

IFCAHSC ALPHA ( IX) )-. 16-09)201.201.202 

201  I  ALPHA* 1 
GO  TO  222 

202  1 ALPHA*2 

222  GO  TO  (203.207) . I  ALPHA 

203  TT*TZERO-T( 1) 

1F( AMOOI TT . H )-H/2» >209.204.204 

204  KH-TT/H41.000001 
GO  TO  207 

205  KH*TT/H+. 000001 
207  CONTINUE 

C 

C — - — — COMPUTE  CONSTANTS- - - — 

0MEGA*FRE0*2.«PI 
RACCAL*SORT( l.-ALPHA(IX)»*2) 

OH»OMEGA*H 
P*OMEGA»RADCAL 
EX*EXP( -ALPHA( IX)*OH) 

COSPH*COS(P*H) 

S INPH*S IN ( P*H) 

A20H*2.» ALPHA ( IX )/OH 
RECOH*l./OH 
ALPHA 2* ALPHA ( I X ) ««2 
XI 2*1 .-2.* ALPHA 2 
FACT 1*EX*S INPH/HAOC AL 
FACT2*EX»C0SPH4FACT1»ALPhA( IX) 

IF(ABS( FREO)-. IE-20  >220.220. 221 

220  1FREQ*1 
GO  TO  223 

221  IFPEQ*2 

223  1F( 1SZ ) 300.206. 30 0 
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200  GO  TO  ( 400 , 209 ) , I FREQ 
209  CONTINUE 


. — COMPUTEXOMEGA(N) .XOOT  FOR  VE LOCI  TV  INPUT - 

GO  TO!224,225>.NlN 

224  XOMEGAI  t)*OMEGA‘XO 
FACT3*RECOH‘! I # -FACT 2) 

FACT4»RECOH‘! .5-2. ‘ALPHA! IX)‘RECOH‘EX‘| ( . 542. ‘ALPHA ( I  X ) *RECOH ) 
1*C0SPH-|XI2*R£C0H- ALPHA! I  X ) /2 . ) *SI NPH/RAOC AL > > 
FACT5*EX»COSPH-ALPHAI I  X ) ‘FACT  t 
F  ACT6*RECOH‘F ACT  1 

FACT7*!RECOH“2>‘! I ,-EX‘I C0SPH‘ I  ALPHA!  I  X > ‘OH/2. ) *SI NPH/RADCAL) ) 

XOOT!l)*XOOTO 

00  210  I»2.N 

XOMEGAI  I IsXOMEGAI  I- I ) *F ACT2‘X00T I  1-1)  *F ACTl-S I  I- 1 > ‘FACT3-S2N0! 1-1  I 
1*FACT4 

210  XDOTI I ) a-xOMEGA ( 1-1 )‘FACTt‘XOOT| I-I)‘FACTS-SI I-i> •FACT6-S2NOI 1-1  )* 
1FACTT 

GO  TO  756 

225  CONTINUE 
KNO*KMX‘KOOK 
KRAsKMN+K ICK 
00  235  1*2.5 
VA«I-I 
TA*VA‘.2‘H 
COSPT*COS!P‘TA> 

S INPT*S IN ! P*TA ) 

EXT* EXP! —ALPHA!  I  X ) ‘OMEGA* T A ) 

FACTRl*EXT‘!COSpT‘ALPHA!  IX)*SINPT/RADCAL) 

FACTR2»EXT*SINPT/RADCAL 
FACTR3-U  .-FACTH  1  I*RECOH 

FACTR4*!TA/H-.5-A20H‘| ALPHA! IX) /2,-X12‘RECOH ) «F ACTR2‘ I ,5‘ A20H ) 

I *EXT*COSPT 1‘RECOH 

XMAXO  *XOMEGA ! KRO I ‘FACTR I ‘XOOT ! KRO )‘FACTR2-Si KRO ) *F ACTR3 

1-S2N0! KRO I ‘FACTR4 
IFiXMAX-XMAXO  >229,230.230 

229  XMAX*XMAXO 

230  XMINO  a X OMEGA I KRA ) ‘FACTR I ‘XOOT( XRA > *F ACTR2-S! XRA )‘FACTR3 

1-S2NC(KRA>‘FACTR4 

IF1XMIN-XMINO  >235,235,232 
232  XMIN>XMINO 
23S  CONTINUE 
N  IN*  1 

211  IF! IP1 >500,212,500 

212  IF! IP2)550, 213,550 

213  IF! IP3I600, 214,600 

214  IF! IP4)650.21S,6S0 

215  IF! IPS) 700, 765, 700 


- COMPUTE  XOMEGAI N ) •  XOOTIN),  FOR  ACCELRATION  INPUT 

300  CONTINUE 

GO  TO 1350, 305) , I FREQ 
305  CONTINUE 

ECCS1*I«-EX‘C0SPH 
0FACT3-! l.-FACT2)/OMEGA 

CFACT4* ! 1 , - A20H*ECOS 1-X  1 2‘FAC  T I ‘RECOH) /OMEGA 
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CO  TO(308»325)»NIN 
306  XOttEGA(l)*OMEGA»XO 

CFACT5*(-4.*ALPhA(IX)-<2.*< 1 • -4 «* ALPHA 2 I /OH-2 •• ALPHA I IX)) *EC0S14 
11X1242. * ALPHA (  IX)»(3.-4.*ALPHA2)/0H)*FACTI >«REC0H/(2.«0MEGA I 
OFACT 6*< l.-ECOSl ) -F ACT 1 »ALPHA ( I  X > 

CFACTTaFACT l/OMEGA 

0FACT8*( I •— FACT2 I *RECOH/QMEGA 

CFACT9«(2*-( 1. 44. • ALPHA ( I  X  I «RECOH >«EC0Sl-( 2. •X124REC0H-ALPHAI I X> ) 
1 *FACT 1 l*RECOH/< 2.40MEGA I 
XDOT ( 1 ) «XDOTO 
OO  310  l>2*M 

XOMEGA( t )aXOMEGA(  I  —  1 ) *F ACT24XD0T {  I-I )*FACT1-Z< I - 1 I *DF ACT3-S ( I - 1) 
l«OFACT4-S2NO(  I-l)*OFACT5 

310  XOOT ( I ) a-XOMEGA ( 1-1 )«FACT14XOOTI  I- 1 ) *DFACT6-Z ( I - 1  I *DF ACT7 
1  -S( I-I I *0 FACT 8- S 2ND < I-I >*DFACT9 

GO  TO  756 
325  CONTINUE 

KR0*KMX4KC0K 

KRA*KMN4KICK 

H2«H»«2 

00  335  I >2. 5 

VA«I-I 

TA=VA*.2*H 

COSPT«COS(P«TA» 

SINPT«SINCP»TA) 

EXT  *EXPI — ALPHA ( IX) *OMEGA*TA I 

FACTRI  =  EXT*( COSPT  4 ALPHA!  I X ) *S INPT/RAOC AL ) 

FACTR2»EXT*SINPT/RADCAL 
FACTR3*( l.-FACTRt l/OMEGA 

F  ACTR4* ( T  A/H— A20H4 ( I . -E XT*COSPT  )-X 1 2»F AC TR2 4REC0H ) /OMEGA 
FACTR5*ITA*»2/H2-TA/H-(2.«< 1 .-4 .4ALPHA2 I •REC0H«*2-A20H ) * 

K  l.-EXT*C0SPT)4(  X124REC0H42.»ALPHAI|X)*(3.-4.*ALPHA2) *REC0H«*2) 
2«FACTR2 1/(2. •OMEGA) 

XMAXOSXOMEGA  (  KRO ) *FACTR 1 4X00T ( KRO ) «F AC  TR2-Z ( KRO ) *F ACTR3 
1-SCKR0)»FACTR4-S2ND(KR0)*FACTR5 
IF( XMAX-XMAXO  >329,330.330 

329  XMAXaXMAXO 

330  XMINO«XOMEGA(KRA)«FACTRI4XOOT(KRA)«FACTR2-Z(KRA)*FACTR3 
l— S( KRAI *FACTR4— S2NOIKRA ) •FACTR5 

IFIXMIN-XMINO  1335.335,332 
332  XMIN>XMINO 
335  CONTINUE 
NIN«t 
GO  TO  211 


— — FREQI-ZERO  FOR  FOUNDATION  ACCELERATION - - 

350  GO  T0(352.36S).NIN 
352  XQMEGAI 1 >»X0 
XOOT 1 1 laXOOTO 
H2»M**2 
00  360  I«2.N 

XOMEGAI DaXOMEGAI  I-I >4H«XDOT( I-ll-ZI 1-1 )*H2/2.-SI I-I )*H2/6.4 
1S2N0I 1-1 )*H2/24« 

360  XOOT I  1 1 Z (  I-I >*H-S< I- 1 )»H/2. 4S2NDI 1-1 )*H/12. 

GO  TO  756 
365  CONTINUE 
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N  IN*  1 
00  TO  211 

- COMPUTE  X.  XOOT  FOR  OMEGA*ZERO — - — VELOCITY  INPUT 

400  GO  TO( 402*425 1  »N  IN 
402  XOMEGAI I )*X0 
XOOT I 11*X00T0 
00  410  1*2. N 

XOMEGA! I|*XOMEGAt I- 1 ) ♦XOOT ( I- l  > «H-S ( I- 1 1 *H*. 5+S2N0 I 1-1 }*H/12. 
410  XOOT ( I  )  =  XDOT ( I— 11-S( I— 11 
GO  TO  756 
42S  CONTINUE 
NIN*1 
GO  TO  211 

- PLOT  X  VS.  T  ANO  XOCT  VS.  T- - - - 

500  GO  TO(SOS.SIO). IFREQ 

505  00  S06  I  * 1 . N 

506  X(II*XOMEGA{II 
GO  TO  S18 

S10  00  SIS  I *1 . N 

518  X(I)*XOMEGA( Il/OMEGA 

516  XMX*Xtll 

XOMX-XOOTI 1 1 
XMN*Xt 1 I 
XOMN*XOOTI 1 ) 

00  530  1*2. N 

I  FI  X  ( I ) —XMX  1521.521.520 

520  XMX*X(II 

521  IFIXIII-XMN  1522.523.523 

522  XMN*X (II 

523  IFIXOOTI 1 1  — XOMX  1525. 525. 524 

524  XCMX*XOOT(II 

525  IFIXOOTI | 1-XDMN  1526. 530.530 

526  XCMN*XOOT(II 

530  CONTINUE 

0 Y* I XMX  -XMN  |/ZN 

CALL  GR I01V (1.TI1I.TINI. XMN.  XMX.  0X.0Y.4 .4 .4 .4 .6 .6 1 
CALL  APLOTVIN.T.X.I, 1.1.42, 1ERI 
IF! I ER 1531.532.531 

531  MRITE 16.9131  ICR. FREQ 

532  NN*N-1 

00  535  1*1. NN 

CALL  LINEVINXVITI I  I l.NYVIXI II  1 ,  NXV I T 1 1 ♦ 1 } 1 , NYV I X  I  I ♦ 1 1  1 1 
535  CONTINUE 

CALL  PRINTVI-15. 15HT1ME  IN  SEC0N0S.4S2 .61 

CALL  PR INTVl -23. 23HFREQ*  CYCLES/SEC. .20 .31 

CALL  LABLVI FREQ. 60 ,3,6* 1.4| 

CALL  PR  I NT V 1-8 »  8H ALPHA*  ,20.171 
C^LL  LABLVI ALPHA! 1X1,72.17,6.  1.  II 
CALL  APRNTVI0.-14.-11.11HX  RESPONSE .4 .600 1 
OY* I XOMX  — XOMN  l/ZN 

CALL  GRI01VI1.TI1I.TIN1.X0MN,  XOMX,  OX ,0Y,4 .4 .4 ,4 *6.61 
CALL  APLOTV  IN.T.XOOT*  1 . 1 . 1 . 42 , IERA 1 
IFIIERA1  538.539.538 
538  MRITEI6.9I4I  IERA. FREQ 
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539  00  540  1*1. NN 

CALL  LINEVINXVITI I) ).NYVIXOOTII  I  1 .NXV( T< I ♦! ) ) ,NVV < XDOT < I ♦ 1 )) ) 

540  CONTINUE 

CALL  PR INTVl— 15. 15HT I ME  IN  SEC0N0S.452 .6) 

CALL  PRINTV<-23.23HFRE0*  CYCLES/SEC. .20.3) 

CALL  LABLV <  FREQ. 60. 3.6.1.4) 

CALL  PRINTV<-8.8HALPHA*  .20.17) 

CALL  LABLVI ALPHA < IX ) . 72 . 1 7. 6. 1 • 1 ) 

CALL  APRNTVI0.-14.-22.22HVEL0CITV  RESPONSE  XDOT. 4. 688) 

00  TO  212 

- PLOT  Y< ACCELERATION)  VS.  TIME - 

550  CONTINUE 

GO  T0I213.SS1). IFREO 

551  CONTINUE 
560  00  562  1*1 .N 

562  X( I )*— XOMEGAt I ) •OMEGA 

XMX-XI 1 ) 

XMN*X( 1 ) 

00  566  1*2. N 

IFIXMX  -X< I  )  >563,564.564 

563  XMX*X( I  ) 

564  IFIXMN  -XII) >566.566.565 

565  XMN*X< 1  ) 

566  CONTINUE 

0 Y*  <  XMX  -XMN  )/ZN 

CALL  GRI01V(1.T(1).T(N).XMN,  XMX.  OX.OV.4.4,4.4.6.6) 

CALL  APLOTVIN. T .X, 1,1. 1.42. IERB ) 

IF  I IERB >570.571.570 

570  MRITEI6.919)  IERB. FREQ 

571  NN*N-1 

DO  575  1*1. NN 

CALL  LINE VI NXV (Till) »NYV! XII)) .NX VI T{ I +1 ) ) , NVV <  X ( I  + 1 ) ) ) 

575  CONTINUE 

CALL  PR  INTV I ”15. 1SHT I ME  IN  SEC0N0S.452 .6) 

CALL  PRINTVC— 23.23HFREQ*  CYCLES/SEC. .25 .3 > 

CALL  LABLVI FREQ. 60. 3,6. 1.4) 

CALL  PR  I NT VI -8. eH ALPHA.  ,20,17) 

CALL  LABLVIALPHAI IX). 72. 17, 6. 1.1) 

CALL  APRNTVIO .“ 14 ,”26,26HABS0LU TE  OR  Y  ACCELERAT ION, 4 , 720 ) 

GO  TO  214 


- - - - - FOUR  CCOROINATE  GRIO - 

600  CONTINUE 
IN1*1«INI 
OSCILI INI)* FREQ 

VELt INI )*AMAX| I ABSI XMAX ) , ABSI XM IN ) ) 

GO  TO  1601.602), l ALPHA 

601  RESIOI INI )*SORT^XOMEGA|KH)442«XOOT|KH)442) 
IFlREStOI INI )“. IE” 19)604, 604, 602 

602  IFIFREQ”.1E”19)604. 604.603 

603  IFIVELI  INI)-. 15-19)604,604, 605 

604  NOFREQ*NOFREQ- 1 
INI. INI-1 

605  IF<  ABSI FRE0-FREG2)” • IE- 09  >606,606 ,214 
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606  CONTINUE 

CAU.  SMXVVI 1*1) 

FRQSML«OSClL( 1 > 

VELSML«VEL( I) 

RE SS ML* RES  1 0( 1 ) 

00  616  I-2.N0FREQ 

IFIOSCILI I ) -FRQSML >611.612.612 

611  FR0SML-0SC1L( I ) 

612  IF( VELt I )— VELSMD613.614.614 

613  VELSML-VELI 1 ) 

614  IFIHESIOI I  >-RESSML)61S.615.616 
613  RESSML«R£SIO( I ) 

616  CONTINUE 

ORDSML-AM IN1( VELSML  .RESSML ) 


TO  FIND  THE  LIMITS  FOR  THE  (4X3)  CYCLES - 

IF( AL0010 (OROSML ) 1620*621 *621 

620  LOGORO-ALOG 10 (OROSML >— 1. 

60  TO  622 

621  LOGORO-ALOG 10( OROSML) 

622  IF( ALOGtO(FROSML) >623*624.624 

623  LOGFRO*ALOG10(FRQSML)-1* 

GO  TO  623 

624  L0GFRQ-AL0G10( FRQSML) 

623  OROSML- 10 •••LOGORO 

OROLG-OROSMLMO<««4 
FRQSML- 10 • ••LOGFRQ 
FR0LG«FRQSML«10.*O 

C — ...PLOT  LOG-LOG  GRIO  FOR  VELOCITY  VS  FREQUENCY - — - 

CALL  GR 101V( 1. FRQSML. FRQLG. OROSML .ORDLG. 1 *0*1*0*1*1*1*1.— 2*-2) 

C - TO  FIND  LARGEST  1X1G  To  THE  PTH  POM£R  LINE  FOR  X— - 

0MSTRT-2.«PI«FR0SML 
TEMP- ALOG 10 ( OROLG/OMSTRT ) 

IF( TEMP >627.626. 628 

627  LOP-TEMP- 1. 

GO  TO  629 

628  LOP-TEMP 

629  01- 10  »**LOP 
V1-0MSTRT401 
MIG- 10*  *OROSML 

XMARGN-NX V ( FRQLG )-NX V( FRQSML ) 

YMARGN»NYV( OROLG )— N VV( MIG) 

SLOPE- YMARGN/XMARGN 
I XC0R-4+ NX V ( FRQSML ) 

CALL  PRINTV(  — 1. 1H0.  I XCOR »NVV( VI ) ) 

CALL  PR INTV(— 2. 2H0- .900. I 7 ) 

CALL  LABLVI 01 .924. 1 7.-2 .1.3) 

C - ORAM  LINES  UP  FROM  VI - - - - — 

VIS-0. 

00  632  1-1*9 
V1S-VIS4V1 

IFIOROLG— VIS >647.630.630 

630  xTCH-NVVIOROlS)-nVV(VIS) 

ITCH-XTCH/SLOPE 

I TCH-NXV I FRQSML ) 4  ITCH 
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K I TCH-NYV ( ORDLG ) 

CALL  LINEV  ( NXV ( FROSML ) t NVV  (  V  IS ) • ITCH.KITCH) 

632  CONTINUE 

C - - CRAW  LINES  FROM  VI  DOWN— - 

647  00  638  J-1.20 
LAL-LOP-J 

L AN-LAL ♦ 1 
OLAN- 10 . **L AN 
CL AL=lO • **L  AL 
00  638  1  =  1.9 
EFT1 *  I 

VIS=OMSTRT*DLAN-OMSTRT»OLAL*EFT  1 
IF< VI S-OROSML >639.633.633 

633  KITCH*NYVCOROLG>-NVV(VIS) 

I YKARG-YMARCN 

IF ( I YMARG-K ITCH >635*634. 634 

634  ZITCH*KITCH 

I TCH*Z I TCH/SLOPE 
GO  TO  637 

635  ITCHsXMARGN 

K ITCH*VMARGN 

637  I TCH- I TCH +N XV  <  FRQSML ) 

K I TCH-K I TCH4NVV ( V  IS  > 

CALL  LINEV  (NXV(FROSML>.NYV(  V  I S 1 » I TCH. K I TCH > 

638  CONTINUE 

C— - DRAW  FROM  THE  RIGHT  OROINATE  THE  REMAINING  X-LINES- - - 

639  0MLAST«2.«PI*FRQLG 
LEFT  2=J 

LEFT  1-EFT 1 
00  641  J-LEFT2.30 

lal-lop-j 

lan*lal+i 

CLAN* 10 •••LAN 
0LAL=10.««LAL 
00  640  I-LEFTI.9 
EFT1-I 

V  I S»OMLAST*DLAN-OMLAST*CL AL*EFT  1 
1FIVIS-0R0SML >642.636. 636 

636  K ITCH-NYVC  V  I S  >-NYV( OROSML ) 

CAPT-KtTCH 

I TCH-C APT/SLOPE 
ITCHaNXVI FRQLG)— ITCH 
KITCH-NYV(OROSML) 

CALL  LINEV  I NXV I FRQLG >  »NYV<  VIS).  ITCH.  KITCH) 

640  CONTINUE 
LEFT1-1 

641  CONTINUE 

642  CONTINUE 

C - -PLOT  POINTS - 

CALL  APLOTVINOFREQ.OSCIL.VCL. 1.1*1 .44.LIT) 

IF (L IT  >643.644.643 

643  WRITEI6.922)  L I T, ALPHA! IX ) 

644  GO  TO  1 648.646). I ALPHA 

648  CALL  APL0TV(N0FREQ*0SCIL.RESI0.1.1.1.38.KIT> 

CALL  PR INTV I -28 . 28H000  -RESIDUAL  SHOCK  SPECTRUM. NXVIFROSML >» 1006 > 
IF (KIT >649.646. 645 
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649  WRITE(6.923)  K 1 T .ALPHA! IX ) 

646  CONTINUE 

CALL  PRINTVI-15. 15HFRE0UENCY  (CPS). 452. 6) 

CALL  APRNTV(0»- 14.-8. 8HVELOC I TY.9.576) 

CALL  PR  I NTV ( — 8. 8HALPHA*  .20.17) 

CALL  LABLV ( ALPHA( IX). 72, 17. 6. 1.1) 

CALL  PRINTV<-19.19H»»»  «SHOCK  SPECTRUM, NXV(FROSML) .1015) 


C - -GENERATION  OF  ACCELERATION  GRID 

C — - TO  FINO  A 1 — - - — 


TEMP*ALOG10(ORDSML*OMSTRT> 

IF( TEMP  >690,851 .851 
850  LGPsTEMP 
GO  TO  8S 3 
e51  L0P*TEMP4l. 

853  A 1= 10 «**LOP 
V1*A1/0MSTRT 

CALL  PR  INTV! -1. 1HA. IXCOR.NYVI Y1 ) ) 

CALL  PR  INTV (-2. 2HA* .900,4) 

CALL  LABL V ( Al, 924. 4*— 2.1*3) 

C - -CRAW  LINES  DOWN  FROM  Al - 

K ITCH=NY V ( OROSML ) 

00  896  1*1.9 
EFT1* l— 1 

VIS*V1»! 1.-.14EFT1) 

IF! OROSML-VIS) 859.895. 857 

855  XTCH*NYV (VI S)—NY VI OROSML) 

1TCH*XTCH/SL0PE 

ITCH* ITCH ♦NXV(FROSML) 

CALL  LINEV  (NXV(FRQSML) «NYV( VIS) .ITCH.  KITCH) 

856  CONTINUE 

C - DRAW  LINES  FROM  Al  UP - - - - - - 

857  DU  870  J*l, 20 
LAN*L0P4J-I 
ALAN*10. ••LAN 
00  870  1*2.10 
EFT1-I 

V I S*AL AN/OMSTRT  «EFT  1 
IF! VIS-OROLG >858.858. 871 

858  KITCH*NYV(VIS)-NVV( OROSML) 

I YMARG*YMARGN 

IF! I YMARG-K ITCH  >862.86 1,861 

861  Z ITCH*K ITCH 

I TCH*Z ITCH/SLOPE 
K I TCH*NYV! OROSML ) 

GO  TO  868 

862  ITCH*XMARGN 

K I TCH«NYV ( V I S )- I VMARG 
868  I TCH* I TCH4NXV I FRQSML ) 

CALL  LINEV  (NXVIFRQSML ) .NYV < V IS ) , I TCH,  KITCH) 

870  CONTINUE 

C - CHAW  FROM  THE  RIGHT  ORDINATE  .REMAINING  A-LINES- 

871  LEFT2*J 
LEFTUEFT1 
K!TCH*NVV(0HDLG) 

00  875  J*LEFT2,30 
LAN*L0P4J— 1 


30 


u  V  U  U  uu 


PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  1  (Continued) 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFNIS) 


06/10/ 


ALAN*10.«*LAN 
00  874  l*LEFT1.10 
EFT  1“ I 

V lS*ALAN/OMLAST»EFT 1 
IF< VIS-0R0LG)872.872.876 
872  C APT*NY V ( ORDLC ) —NYVlVIS) 

I TCH»C APT/SLOPE 
ITCH*NXV(FRQuG)-I TCH 

CALL  LtNEV  INXVIFROLG) ,NVV|  VI  S)  .  ITCH. KI TCH) 
874  CONTINUE 
LEFT  1*2 
879  CONTINUE 
876  CONTINUE 
00  TO  214 


- SHORT  FORM - 

650  CONTINUE 

IF( ABSIFREQ-FREQ I  I—  •  IE* 04 ) 651 .651.652 

651  MRITEI8.912)  ALPHA!  IX) 

652  00  T0(658.655). I  ALPHA 

655  X0MAXSAMAX1 ( ABS( XMAX ) . ABS!  XM  IN)  ) 

GO  TOI668.669). IFREO 
658  RMAX* SORT ( XOMEGA ( KH )**2*XD0T I KH )**2 ) 

X0MAX>AMAX1 ( ABS ( X MAX  )  . ABS! XMIN) .ABS(RMAX) ) 

GO  T0(681. 678). IFREO 

668  HR  I TE ! 8.906 )  FREO.XOMAX 
GO  TO  215 

669  MRITE! 8.90S )  FREO.XOMAX 
GO  TO  215 

678  MR! TE ( 8.907 )  FREO.XOMAX. RMAX 
GO  TO  215 

681  MRITE! 8.908  )  FREQ .XOMAX.RMAX 
GO  TO  215 

- - LONG  FORM- - 

700  CONTINUE 

I F I  ABS ( FREO-FREO 1 )- . IE-04 ) 70 1 . 70 1  • 702 

701  MR ITE(6 .9 IS )  ALPHA! IX) 

702  GO  TO! 704.705). I ALPHA 

704  RM AX*SQRT I XOMEGA ! KH  I  ••2»X00T| KH ) **2 ) 

GO  T0(743.742). IFREO 

705  GO  T0I726. 727). IFREO 

726  MRITEI6.916  )  FREQ • XMAX. XMIN 
GO  TO  785 

727  MRITEI6.910)  FREQ  .  XMAX. XM  IN 
GO  TO  785 

742  MR  I TE 1 6.910 )  FREQ. XMAX. XMIN. RMAX 
GO  TO  785 

743  MRITEI6.91 1 )  FREQ, XMAX. XMIN. RMAX 
GO  TO  785 

>— — —PRELIMINARY  CALCULATIONS  FOR  MIN.  AND  MAX.  RESPONSE 
750  GO  TO  1758.757). I ALPHA 

757  NX>N 

GO  TO  759 

758  NX>KH 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  1  (Continued) 


SNERD  -  EFN  SOURCE  STATEMENT  -  IFN(S) 


06/1 


759  XMAX*XOMEGA(  1) 

KMX*  1 

XM!N*XOMEGA< 1 ) 

KMN*1 

00  760  1*2 .NX 

IF(XMAX-XOMEGA(  11)761.762.762 

761  XMAX*XOMEGA( I  ) 

KMX*  I 

762  IF(XMIN-XOMEGA( I ) >760.765.  765 
765  XMIN*XOMEGA( I ) 

KMN*  I 

760  CONTINUE 

IFCKMX-l I77I.77C .771 

770  KMX*KMX4l 

771  IF ( KMX-NX  >773.772.773 

772  KMX*KMX-1 

773  IF(KMN-1 >775.774,775 
77*  KMN-KMN41 

775  1FCKMN-NX1777.776.777 

776  KMN*KMN-1 

777  CONTINUE 

IF <  XOOT ( KMX— 1  ) *XDOT (KMX) )77fl.77E,779 
776  KOOK*— 1 
GO  TO  781 
779  KOOK*0 

761  CONTINUE 

I  F<  XOOT  (KMN-  1  > 11  XOOT  (KMN  )  >783. 783. 764 

783  KICK*— 1 
GO  TO  950 

784  KICK-0 
950  N IN— 2 

GO  TO  2v7 
C 

C— — — EITHER  RESTART  WITH  NEW  FREQ.  OR  RESTART  WITH  NEW  ALPHA  OR  STOP 
765  CONTINUE 

IF( ABSC  FRE0-FRE02 )- . IE-03  > 600 ,755. 755 
755  FRE0*FREQ40ELTAf 
GO  TO  207 

600  CONTINUE 
GO  TO  190 

601  WRITE ( 6.909 ) 

CALL  FRAMEV 
STOP 

C 

900  FORMAT! 1H1 • 49X. 22HSHOCK  SPECTRUM  PROGRAM//50X.21H«PPLIED  MATH  LAB 
10TMB. /// ) 

901  FORMAT ( 40X, 36HF0UN0ATI0N  ACCELERATION  VS  TIME!  INPUT)//18X.23HF0UN0 
1 AT  I  ON  ACCELERATION. 36X. I5HTIME  IN  SECONDS) 

902  FORMAT ( 42X • 3SHF0UN0ATI0N  VELOCITV  VS  TIME  C INPUT) //20X • 19HF0UN0AT I 
ION  VELOC ITV.3SX.1SHTIME  IN  SECONDS) 

903  FORMAT! IX. F3S.9.  44X.f20.91 

904  FORMAT! 47H  PLOTTING  ERROR  NUMBER  OF  POINTS  OUT  OF  RANGE*.I2//> 

905  FORMAT! 1SX.F10.S.29X.F12. 7) 

906  FORMAT! 10X. I 1HFIRST  FREO-.F 10 .5 .7X. 1 3HMAX  RESPONSE* ,F12. 7/ ) 

907  FORMAT ( 1 SX.F 10.5.2! 29X.F1 2.7 1  ) 

908  FORMAT! 10X* 11HFIRST  FREO-.F 10.S.7X.I 3HMAX  RESPONSE- .F 1 2.7. 1 2X . 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  1  (Continued) 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFN(S) 


06/10/ 


117HMAX  RES  RESPONSE-.F 12.7 ) 

909  FORMAT ( ///38X. 25HPROGRAM  RAN  TO  COMPLETION) 

910  FORMAT ! 10X.F10.S.3!  18X.FI2.7)  ) 

911  FORMAT ( SX • 1 1HF IRST  FREO-.F 10. St 9X .9HMAX  RESP- .F 12. 7 . 10X .9HMI N  RESP 
1-.F12.7.7X. 1SHMAX  RES10  RESP-.F12.7) 

912  FORMAT!  1H1.SSX. 10HSH0RT  FORM//4SX . 20 HD AMP I NG  COEFF IC IENT* .F10.  8/// 
1/9X • 22HFREQUENCY! CYCLE  S/SEC • 1 . 1 9X . 20HMAX  RESPONSE  X  OMEGA. 18X.26HM 
2AX  RESIO  RESPONSE  X  OMEGA) 

913  FORMAT (S9H  PLOTTING  ERROR  FOR  X  VS  T.  NUMBER  OF  POINTS  OUT  OF  RAN 
1GE-.I2.6H  FREQ*.F10.7//I 

914  FORMAT ( IX* I  2 . 43H  POINTS  OUT  OF  RANGE  FOR  XOOT  VS  T I  ME .FREQ*. FI 0.  7 
1//) 

915  FORMAT!  1H1 .S5X.9HLONG  FORM//* 5X, 2 OHO AMP I NG  COEFF IC  IENT- .F 10. B//// 
17X  * 16HFREQ! CVCLES/SEC 1  * 12X  * 20HMAX  RESPONSE  X  OMEGA. 10X.20HMIN  RESP 
20NSE  X  OMEGA. BX. 22HMAX  RESIO  RESP  X  OMEGA) 

916  FORMAT  1 4X . 1 1HF IRST  FREQ- .F 1 0. 5. 10X  .  9HMAX  RESP- .FI  2. 7. 10X.9HMI  N  RES 
IP- .F 12. 7 ) 

917  FORMAT!  1H1 . 16X. 1HI.3BX.4HSI I ) .34X.7HS2ND! 1 ) // 

1! 17X. 14 . 34 X .F12.7.28X.F12.7)) 

918  FORMAT! 1H1 .SOX. 19H INTERMEDIATE  VALUES/ 19X. 1HI .3SX.9HX0MEGA! I )  . 

1 3 1 X . 7HXDOT ! I )// 

2! 17X. I3.34X.F12.7.28X.F12.7)) 

919  FORMAT! IX. 12. 43H  POINTS  OUT  OF  RANGE  FOR  V  VS  T IME .FREQ-.F 10. 7 
I//) 

920  FORMAT! 3F20. 8) 

921  FORMAT ( 6F 1 S.6 ) 

922  FORMAT !  i  X • 1 2.49H  POINTS  OUT  OF  RANGE  FOR  VEL  VS.  FREQ  WITH  ALPHA-. 
1F10.8) 

923  FORMAT!  IX. 12. SOh  POINTS  OUT  OF  RANGE  FOR  RESIDUAL  PLOT  NITH  ALPHA- 
1.  F 10.8 ) 

ENO 
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ATION  ACCELERATION 


129.09 


Figure  3  -  Input  Shock  Function  (Sample  Problem  1) 
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Figure  4  -  Four-Coordinate  Grid  for  Sample  Problem  1,  a  =  0 

FREQUENCY  RANGE  =  10  TO  1000  CPS 
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Figure  5  —  Four-Coordinate  Grid  for  Sample  Problem  1,  a  >.  0.1 

FREQUENCY  RANOE  =  10  TO  1000  CPS 
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Figure  6  -  Four-Coordinate  Grid  for  Sample  Problem  1,  a  c  0.5 

FREQUENCY  RANOE  =  10  TO  1000  CPS 


TABULATION  OF  SHOCK  FUNCTION  FOR  SAMPLE  PROBLEM  1 

SHOCK  SPECTRUM  PROCRAM 
APPLIED  MATH  LAB  OTMB 


FOUNDATION  ACCELERATION  VS  TIME!  INPUT) 


FOUNDATION  ACCELERATION 
0. 

7.494259298 

14.609627389 

22.196642709 

29.291196346 

36.272976239 

43.009920267 

49.719232649 

96.147649769 

62.369493907 

66.366116266 

74.131931309 

79.690019646 

64.912377397 

69.909680636 

94.634310722 

99.070369326 

103.239664366 

107.100766069 

110.669174194 

113.937322617 

116.902992699 

119.963299179 

121.916684999 

123.960916693 

129.719063099 

127.199060909 

126.303602490 

129.192914047 

129.710926096 

129.963192390 

129.979661309 

129.699392994 

129.149709702 

126.346902996 

127.296049303 

126.006376266 

124.486022604 

122.791962119 

120.606099692 

116.669028202 

116.346379260 

113.692384967 

111.199636499 

100.400992393 

109.469936761 

102.416932372 

99.261381149 

96.011973791 

92.682640699 

89.266141291 


TIME  !N  SECONDS 

0. 

0.000^29996 
0.000060000 
O.OOOC69996 
0 .000 1 19999 
0. *00149999 
0.000179999 
0.0002C9999 
0.000239998 
0.1*00269994 
0.000299996 
0.000329994 
0.000399997 
0.000369993 
0.000419997 
O.C00449993 
0.000479996 
0.000910000 
0.000539996 
0.000969999 
0.000599999 
0.000629999 
0.900659995 
0.000689998 
0.000719994 
0.000749998 
0.000779994 
0. 900 Ou 9997 
0.900839993 
0.900869997 
0.000899993 
0.900929996 
0.900960000 
0.000989996 
0.00 1C  19999 
0.001049999 
0 . 'C 10  79999 
0.001109999 
0 .90 1 139998 
0.001169994 
0.001199996 
0.301229994 
0.301299997 
0.301269994 
0.. '01319997 
0.901349993 
0.001379997 
0.C01409993 
0. 301439996 
0.001470000 
0.001499996 
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PRINTOUT  OF  FIRST  AND  SECOND  FORWARD  DIFFERENCE  FOR  SAMPLE  PROBLEM  1 

I  5(11,  S2N0(I>  s2 


1 

7.4942993 

-0.1028912 

2 

7.3913681 

-C. 1028912 

3 

7.2712193 

-0.1201528 

4 

7.1343936 

-C. 1368617 

S 

6.9813819 

-0.1529717 

6 

6.8129420 

-0.1684399 

7 

6.6297126 

-C. 1832294 

a 

6.4324169 

-0.1972957 

9 

6.2218041 

-0.2106128 

10 

9.9986644 

-C. 2231398 

11 

9.7638130 

-0.2348514 

12 

9.9180883 

-0.2457247 

13 

9.2623977 

-0.2957306 

14 

4.9975033 

-0.2648544 

19 

4.7244301 

-C.2  73C732 

1« 

4.4440546 

-0.2803755 

17 

4.1972990 

-0.2867596 

ia 

3.8651037 

-0.2921953 

19 

3.5684061 

—0.2966975 

20 

3.2681484 

-0.3002977 

21 

2.9652700 

-0.3028784 

22 

2.6607065 

-0.3045635 

23 

2.3593858 

-0.3053207 

24 

2.0902319 

-0.3051536 

29 

1.7461462 

-0.3040857 

26 

1.4440174 

-0.3021288 

27 

1.1447220 

-0.2992954 

2a 

0.8491116 

-0.2956104 

29 

0.9980120 

-0.2910995 

30 

0.2722263 

-0.2897857 

31 

-0.00  74711 

-0.2796974 

32 

-0.2803288 

-0.2728577 

33 

-0.9456429 

-3.2653141 

34 

-0.0027267 

-0.2970839 

39 

-1.0509377 

-3.2482109 

36 

-1.2896690 

-0.2387314 

37 

-1.9183535 

*0.2286844 

3a 

-1.7364607 

-y. 2181372 

39 

-1.9435024 

-0.2073417 

40 

-2.1390314 

-0.1995290 

41 

-2.3226490 

-0.1836176 

42 

•2.4939947 

-0.1713457 

43 

-2.6927481 

-0.1987534 

44 

-2.7986441 

-0.1458960 

49 

-2.9314956 

-0.1328119 

46 

-3.0510044 

-0.1195488 

47 

-3.1571912 

-0.1061468 

48 

-3.2498074 

-0.0926561 

49 

-3.3289291 

-0.0791178 

90 

-3.3945074 

-0.0655823 

91 

-3.4465914 

-0.0920840 

92 

-3.4892648 

-0.0386734 

93 

-3.5106916 

-0.0253868 

94 

-3.9229206 

-9.0122690 

99 

-3.5222907 

0.0006399 

96 

-3.5089779 

0.0133028 

97 

-3.4832994 

0.0256824 
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SAMPLE  PROBLEM  1 


shout  form 

DAMPING  COEFFICIENT 


:1 


TNWOKT I  CYCLES/ SCC  .  I 

IS.SSOOS 

to  •M000 

90.00000 
40.00000 
80.00000 
80.00000 
TO .00000 
00.00000 
00.00000 
100.00000 
110.00000 
110.00000 
130.00000 
140.00000 
180.00000 
140.00000 
170.00000 
100.00000 
140.00000 
100.00000 
910.00000 

too. 00000 

130.00400 
940.00040 
980.00040 
940.00000 
970.00000 
940.00004 
940.00000 
300.00040 
310.00000 
330.00440 
330.00004 
340.04400 
380.00000 
340.00000 
370.00000 
380.00000 
340.00040 
400.00000 
410.00400 
490.00040 
430.00000 
440.00400 
480.00000 
440.00000 
474.00000 
400.40000 
440.00000 
800 .00000 
810.00000 
890.00400 


max  Response  x  omcoa 
0.9079498 
0.2139544 
0.93*3971 
4.9809*00 
0.98T4S04 
4.4483*40 
4.3944093 
0.34*4747 
0.40 94*38 
0.344*411 
0.328*349 
0.9942964 
0.2041444 
0.1740678 
0.16*3494 
0.1612798 
0.19616*4 
0.1991390 
0.1310937 
0.1*06604 
0.1*57131 
0.1*29994 
0.139*123 
0.193*681 
0.1719*7* 

0. 1S9S0S* 
0.1961798 
0.4170794 
0.43*1441 
0.4*39210 
0.4399940 
0.41701*4 
0.14446*1 
0.1044044 
0.1641084 
0.1976931 
0.1*89047 
0.1330379 
0.1443907 
0.1444743 
0.11*0373 
0.1077670 
0.1019402 
0.0493746 
0.0993973 
0.993908* 
0.0706344 
0.47*9374 
0.0780430 
0.073*003 
0.0717896 
0.0701900 


max  Res to 


espoNse  x  omega 

.2079929 

.4139996 

.4198**8 

.2491704 

.491800* 

.2703299 

.30273*2 

.3994694 

.*029*39 

.392*702 

.343137* 

.4471617 

.1407629 

. 1046990 

.07*9*00 

.0904*7* 

.4*426*7 

.4*07670 

.0*3461 1 

.0846901 

.0*34114 

.0739*99 

. 10*906* 
.1434144 
.146199* 
.1739992 
.40*3399 
.4303319 

.4497889 

.10088*9 

.1999079 

.1383144 

.11737*4 

.143384* 

.0449031 

.07**132 

.0*0*170 

.0*47487 

.0873943 

.0932774 

.0*41904 

.0*3061* 

.0*07940 

.03*4173 

.039749* 

03*0900 


V.  0193618 


40 


930*00000 

0.0304324 

940.00000 

0.0471014 

0*0290216 

990.00000 

0.0494341 

0.0870279 

900.00000 

0.0443174 

0.0263309 

9T0 .00000 

0.0434349 

0.0290704 

900.00000 

0o0014966 

0.0242022 

990.00000 

0.0401929 

0.0230929 

000  .09000 

0.0509277 

0*0219662 

010.00400 

0.0976940 

0.0212421 

030.00000 

0o0505023 

0.0204936 

030.00000 

0.0953399 

0*0196742 

040.00000 

0.0542099 

0.0144002 

090.00000 

0.0931119 

0.0102437 

090.00000 

0.0930434 

0*0174431 

670.00000 

0.0910093 

0.0167792 

090.00000 

0*0490952 

0.0163960 

090.00000 

0.0440134 

0.0197409 

700.00000 

0.0440942 

0*0151047 

710.00000 

0.0471344 

0.0147314 

730.00000 

0.0442192 

0.0142906 

730.00000 

0*0453304 

0.0137106 

740.00000 

0.0444794 

0*0133310 

790.00000 

0*0436436 

0.0130249 

700.00000 

0.0424304 

0.0129393 

770.00000 

0.0420341 

0.0121319 

700.00000 

0.0412493 

0.0119004 

790.00000 

0.0409114 

0.0119239 

900. 00000 

0.0347744 

0.0111103 

910.00000 

0.0390437 

0.0106992 

930.00000 

0.0303440 

0.0106309 

030.00400 

0.0374494 

0.0102409 

040.00000 

0.0370394 

0.0100090 

090.00000 

0.0343404 

040.00000 

0.0397903 

0.0094932 

070.00000 

0.0391344 

0.0092249 

940.00000 

0.0349340 

0.0090942 

090.00000 

0.0339910 

0.0094363 

900.00000 

0.0333744 

0.0064906 

910.00000 

0.0324210 

0.0044174 

430.00000 

0.0332779 

930.00000 

0.0317440 

940.00000 

0.0311171 

0*0076119 

490.00000 

0.0307201 

0.0476974 

940.00000 

0.0302247 

970.00000 

0.0207409 

0.0072724 

400.00000 

0.0392479 

0.0071602 

440.00000 

0.0244091 

0.0076000 

1000.00009 

0.0843933 

0.0064009 
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SAMPLE  PROBLEM  1 


SHOUT  FORM 

OAMPING  COEPF 10  tENTM. 00900000 


FREQUENCY  1CVCLES/SEC. I 

M. OOOOO 
SO.  MOtt 
40.00000 
00.00000 
00.00000 
ro.ooooo 
00.00000 
00.00000 
too .00000 
ito.eooto 

I SO .00000 
I SO .00000 
140.00000 
100.00000 
100.00000 
ITO.OOOOO 
100.00000 
100.00000 
too. 00000 

t 10. 00000 

tao .ooooo 

ISO .00000 
too. ooooo 
tOO. MOOO 
too .00000 
ITO.OOOOO 
tM.00000 
too. ooooo 
300.00000 
S10.M000 
sto.ooooo 
330. MOOO 
340.00000 
3M.M0M 
300.00000 
3TO.OOMO 
300. MOOO 
3M.MM0 
400 .00000 
410.00000 
4M.OOOOO 
430 .MOOO 
440.00MO 
400.00000 
440.00000 
4TO.MOOO 
4M.OOMO 
4M.M000 
MO  .OOOOO 
010.00000 
OM.MOOO 


MAX  RESPONSE  X  OMEGA 
O.tTTOtTA 
O.ITOOSSS 
0.3041074 
0.S169432 
O.ttTtSAO 
0.2S47007 
0.2I3S1O1 
0. 1313120 
0.1230210 
Q.20SS132 
0.1494174 
0.1444007 
0.1470100 
0.1440704 
0. 1412331 
0.13041S0 
0.1340210 
0.1342103 
0.1314370 
0.1104007 
0.1100170 
O.ltatOlt 
0.1103031 
0.1174404 
0.1 140040 
0.1113314 
0.1004743 
0.1100744 
0.1110343 
0.1104107 
0.1003103 
0.1047713 
0.1M1M1 
0.0040008 
0.0004002 
O.MOMOt 
0.0044070 
0.M44M3 
0.0010143 
0.0M4M7 
0.0707004 
0.0740740 
0.0701101 
0.0730343 
0. 0710774 
0.0701001 
0.04071 00 
0.047IOM 
0.0407107 
0.0041040 
0.0410001 
0.0010077 


MAX  MES10  RESPONSE  X  ONES* 
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0902279 

940 *#0000 

0909709 

990*00000 

0977129 

900*00000 

0999292 

970*00000 

0 

0993420 

900.00000 

0942143 

990.00000 

0931090 

900.00000 

9920244 

910*00000 

0909099 

020.00000 

0499490 

OJO.OOOOO 

0409034 

940*00000 

0400091 

990*00000 

0470742 

990.00000 

0491730 

970.00000 

0492927 

900.00000 

0444191 

990.00000 

0439003 

700.00000 

•407374 

710.00000 

0419693 

720.00000 

0411092 

730.00000 

0404099 

740.00000 

0396997 

790.00000 

0399910 

790.00000 

0302937 

770.00000 

0379924 

700.00000 

0369379 

790.00000 

0392904 

000.90000 

0399379 

910.00000 

0390391 

020.00000 

0344396 

030.00000 

0339399 

040.00000 

0332977 

090.00000 

0327227 

090.00000 

0321773 

070.00000 

0319322 

900.00000 

0311124 

too.  00k  00 

0309104 

900.00000 

0301249 

910.00000 

0299319 

020.00000 

0291900 

930.00000 

0207000 

040.00000 

0202999 

990.00000 

0279112 

990.00000 

0273904 

070.00000 

0299470 

900.00000 

0299499 

990.00000 

0291499 

1900.00000 

0297409 
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SAMPLE  PROBLEM  1 


SHOUT  PORN 

DAMPING  COEPPIClENTaO'SOOOOOOO 


p*eouENCviCTCt.es/  sec.  i 
IO.SOOOO 
SO.OOOOO 
30*00000 
40.00000 
SO.OOOOO 
SO.OOOOO 

ro.ooooo 

SO.OOOOO 
•0.00000 
100.00000 
110.00000 
ISO  .00000 
130.00000 
140.00000 
180.00000 
|40.00o00 
170.00000 
100.00000 
>•0.00000 
ISO .00000 
310.00000 
330.00000 
330.00000 
340.00000 
3SO.OOOOO 
360.00000 
370.00000 
300.00000 
3*0.00000 
300.00000 
310.00000 
330.00000 
330.00000 
340.00000 
3S0. 00000 
360.00000 
370.00000 
300.00000 
3*0.00000 
400.00000 
410.00000 
430.00000 
430.00000 
440.00000 
4S0 .00000 
460.00000 
470.00000 
400.00000 
4«0 .00000 
800. 00000 
810.00000 


MAS  RESPONSE  X  OMEGA 
6. 6*86034 
0.1743736 
0.137*774 
0.1738866 
0.1741084 
0.1*001*8 
••144**70 
0.13*6800 
0.1187130 
0. 1081648 
0.0*72676 
0.0*832*4 
0.0*74387 
•.0*8*716 
0.0*40406 
0.0*1*184 
0.08*83** 
0.0877301 
0.0086670 
0.0836308 
0.0016244 
0.07*7440 
0.0770310 
0.0760846 
0.0743176 
•••73624* 

0. 070*7*6 
0.06*3832 
0.0670387 
0.0663363 
0.0640032 
0.0634730 
0.0621083 
0.0607742 
0.08*4771 
0.0802032 
0.0871008 
0.088*0*0 
0.08480*1 
0.0830074 
0.0S377O* 
0.0817840 
0.0807732 
0.04*0104 
0. 0400764 
0.047*060 
0.0470020 
0.0403100 
0.0483041 
0. 04*8610 
0.0437786 
0.0*2**82 


MAX  RESID  RESPONSE  X  OMEGA 
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0.0422499 

940.00060 

0.04,19147 

990.00000 

0.0*07990 

990.00000 

0*0401094 

970.00000 

0.0394272 

990.00000 

0.0397692 

990.00000 

0.0391319 

600.00000 

0.037900* 

910.00000 

0.0349902 

920.00000 

0.0393021 

930.00000 

0.0397199 

9*0.00000 

0.0391494 

990.00000 

0.0349066 

990.00000 

0*0340700 

970.00000 

0.0339417 

600  oOOOOO 

0.0330300 

690.00000 

0.0329393 

700.00000 

0.0320929 

710.00000 

0*0319734 

720  oOOOOO 

0.0311007 

730.00000 

0.0309929 

740.00000 

0*0302233 

790.00000 

0.0297909 

790.00000 

0.0293949 

770.00000 

0.0299960 

790.00000 

0.0299996 

790.00000 

0.0291691 

000.00000 

0.0277949 

410.00000 

0.027*071 

920.00000 

0.0270399 

930.00000 

0*0266072 

040.00000 

0.0263409 

990.00000 

0.0299990 

990.00000 

0.0296691 

970.00000 

0.0293364 

990.00000 

0.0290130 

0.02*70*9 

900.00000 

0.02**022 

910.00000 

0.0241094 

920.00000 

0.0239140 

930 .00009 

0.0239290 

9*0.00000 

0.0232479 

990.00000 

0.0229723 

990.00000 

0.0227020 

470.00000 

0.0224*02 

990.00000 

0.0221062 

990.00000 

0.0219369 

0.0216923 

45 


SAMPLE  PROBLEM  1 


Mtatcvci.es/sec  > 
10.00000 
20.00000 
30.00000 
*0.00000 
90.00000 
60.00000 
70.00000 
00.00000 
90.00000 
100.00000 
110.00000 
120.00000 
130.00000 
1*0.00000 
ISO. 00000 
160.00000 
170.00000 
100.00000 
190.00000 
200.00000 
210.00000 
220.00000 
230.00000 
240.00000 
290.00000 
260.00000 
270.00000 
200.00000 
290.00000 
300.00000 
310.00000 
320.00000 
330.00000 
3*0.00000 
390.60000 
360.00000 
370.00000 
300.00000 
390.00000 
400.00000 
410.00000 
420.00000 
430.00000 
440.00000 
*90.00000 
*60.00000 
*70.00000 
400.00000 
490.00000 
900.00000 
910.00000 
920.00000 


LONO  porn 

OANPING  COEMlCItMTaO. 


MAX  RESPONSE  X  ONES* 
0.0000029 
0.1116299 
0.2197229 
0.2299621 
0.2901399 
0.2993460 
0.3222017 
0.36*2797 
0.3901919 
0.3924911 
0.329*392 
0.2929902 
0.2091999 
0.1790679 
0. 1*70611 
0.1139742 
0.0913069 
0.0997927 
0.0992999 
0.06*999* 
0.0912001 
0.10  7*940 
0.1319333 
0.1934991 
0.1719*74 
0.1999094 
0.1961799 
0.2139926 
0.23*1621 
0.2*39210 
0.2390194 
0.21701*9 
0*1949299 
0.1909992 
0. 1691 996 
0.1976931 
0.1*99097 
0.1330372 
0.120*794 
0.10904*2 
0.0999097 
0.0942173 
0.0730973 
0.0626136 
0.0616163 
0.060*020 
0.0992929 
0.0993163 
0.0919079 
0.0919449 
0.091*743 
0.0911100 


MIN  RESPONSE  X  OMECA  MAX 
-0.2069933 
-0.2091239 
-0.2343271 
-0.2S02409 
-0.2970906 
-0.2669993 
-0.3229092 
-0.3611291 
-0.3990324 
-0.3992212 
-0.3292009 
-0.2962969 
-0.2029997 
-0.167621* 

-0.1643999 
-0.1612799 
-0.1991949 
-0.1991090 
-0.1919937 
-0.1499609 
-0.1497131 
-0.1429999 
-0.139*123 
-0.1362791 
-0.1377002 
-0.1699999 
-0.1961367 
-0.2170799 
-0.2336992 
-0.2*33912 
-C. 2399290 
-0.2169204 
-0.199*6*1 
-0.1791679 
-0.19*9310 
-0.1*31226 
-0.13793*3 
-0.1323199 
-0.1263997 
-0.12027*3 
-0.11*0373 
-0.1077670 
-0.1019292 
-0.0993766 
-0*0993973 
-0.003909* 

-0.070632* 

-0.0769374 

-0.0790939 

-0.073*003 

-0.0717996 

-0.0701900 


resio  *esp  x  omeca 

0.2079929 

0.2139996 

0.21994*9 

0.229170* 

0*2919094 

0.2703299 

0.3027342 

0.399299* 

0. *029439 

0.3924702 

0.323137* 

0.2271617 

0.1497629 

0.1049990 

0.07*9*00 

0.0909*74 

0.04226*7 

0.0*07670 

0.0*32611 

0.0929901 

0.0632112 

0.0739499 

0.0992626 

0.10*900* 

0.123*190 

0.1*6199* 

0.1739992 

0.2043399 

0.2303319 

0.2422160 

0.23*0*09 

0.2097999 

0.19099*9 

0.1999979 

0.1393160 

0.11737*2 

0.1033906 

0.092616* 

0.0929031 

0.07*6132 

0.069*170 

0.0627997 

0.0973923 

0.093277* 

0.0499002 

0.0461209 

0.0*30616 

0.0*07900 

0.0392173 

0.039799* 

0.03*0900 

0.0323619 
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930*00000 

0.0901319 

940.00000 

0.0407903 

990.00000 

0.0469905 

900.00000 

0.0449161 

570. 00000 

0.0429542 

900.00000 

0.0399636 

990.00000 

0.0371932 

600.00000 

0.0343024 

610.00000 

0.0313290 

620.00000 

0.0292234 

630.00000 

0.0202349 

640.00000 

0.0207046 

690.00000 

0.0290902 

660.00000 

0.0293749 

670.00000 

0.0295037 

660.00000 

0.0294607 

690.00000 

0.0292600 

700.00000 

0.0200793 

710.00000 

0.0203209 

720.00000 

0.0276102 

730.00000 

0.0267590 

740.00000 

0.0297604 

790.00000 

0.0246569 

760.00000 

0.0234433 

770.00000 

0.0221469 

700.00000 

0.0207001 

790.00000 

0.0193360 

000.00000 

0.0196544 

010.00060 

.0.0199203 

020.00000 

0.0201350 

030.00000 

0.0202630 

040.00000 

0.0203002 

090.00000 

0.0202632 

060.00000 

0.0201299 

070.00000 

0*0199066 

060.00000 

0.01 95903 

090.00000 

0.0192060 

900.00000 

0.0107363 

910.00000 

0.0101950 

920.00000 

0.0175060 

930.00000 

0.0169109 

940.00000 

0.0161979 

990.00000 

0.0194303 

960.00000 

0.0140003 

970.00000 

0.0190021 

900.00000 

0.0191531 

990.00000 

0.0!5i'592 

1000. 00000 

0.0153194 

*.0304324 

-0.0671016 

0.0290*16 

-0.0696301 

0.0270279 

*0*0642174 

0.0263300 

-0.0620369 

0.0230700 

-0.0614966 

0.02*202* 

-0.0601925 

0.02309*3 

-0.0569277 

0.021966* 

-0.0976900 

0.021*4*1 

-0.0969023 

0.020*336 

-0*0553399 

5*0 194741 

*0*0542099 

0.0100002 

-0.0931115 

*.0102437 

-0.0320436 

0.017**31 

-0.0510052 

0.016779* 

-0.0499952 

0.0163960 

-0.0490126 

0.0197*09 

-0.0400562 

0.0131047 

*0*0471246 

0.01*7314 

-0.0462192 

0.0142906 

-0. 049330* 

0.0127106 

-0.0444796 

0.0133310 

-0.0436436 

0.0130209 

*0*0426304 

0.0129393 

-0.0420361 

0.0121319 

-0.0412653 

0.011900* 

-0.040511* 

0.0113235 

-0.0397796 

0.0111103 

-0.0390637 

0.010099* 

-0.0303600 

0.0106309 

-0. 0376006 

*.010*409 

-0.0370256 

0.0100090 

-0.0363006 

0. 0090240 

-0.0337503 

*.009493* 

-0.03313*6 

0.0092256 

-0.03*3360 

*.009009* 

-0.0339510 

0.0000363 

-0.033379* 

0.0009900 

-0.0320210 

0.0004170 

-0.Q322770 

0.000*439 

-0.0317460 

*.0079604 

-0.0312271 

0.0070119 

-0.0307201 

0.0076976 

-0.03022*7 

*.0074609 

-0.0297405 

0.00727*0 

-0.0292675 

0.007100* 

-0.0200051 

*.0070000 

-0.0203933 

0.« 0000*9 
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SAMPLE  PROBLEM  1 


LONG  FORM 

DAMPING  COEFFICIENTS.  09999999 


FREOICVCLES/SEC)  MAX  RESPONSE  X  OMEGA 


10.00000 

0.1390010 

20.00000 

0.1304083 

30*00900 

Qo 1385235 

40.00000 

0.  1403955 

50.00000 

0.1639449 

60,00000 

0.2009101 

70.00000 

0.2235191 

60.00000 

0.2313129 

90.00000 

0.2230216 

100.00000 

0.2055132 

110.00000 

0.1056170 

120.00000 

0.1664997 

130.00000 

0.1479100 

140.00000 

0.1207061 

150.00000 

0.1006053 

160.00000 

0'0ta33S7 

170.00000 

0.0693412 

100.00000 

0.0543430 

190.00000 

0.0470772 

200.00000 

0.0513437 

210.00000 

0.0603020 

220.00000 

0.0700200 

230.00000 

0.0011709 

240.00000 

0.0904091 

250.00000 

0.0902499 

200.00000 

0.1043040 

270.00000 

0.1005026 

200.00000 

0.1100744 

290.00000 

0.11 15363 

300.00000 

0.11 06 1 07 

310.00000 

0.1003103 

320.00000 

0.1047713 

330.00000 

0.1002321 

340.00000 

0 .0940908 

350.00000 

0.0000957 

340.00000 

0.0024535 

370.00000 

0.0737169 

300.00000 

0.0600C93 

390.00000 

0.0610631 

400.00000 

0.05S002S 

410.00000 

0.0403075 

420.00000 

0.0410500 

430.00000 

0.0357064 

0.0331957 

400.00000 

0.031 3SSI 

400.00000 

0.0295201 

470.00000 

0.0276970 

400.00000 

0.0250992 

490.00000 

0.0241330 

800.00000 

0.0234211 

S10. 00000 

0.024  0  307 

520.00000 

0.0242630 

M|N  RESPONSE  X  OMEGA  MAX  RESIO  RESP 

-0.1770270 
-0.I700SSS 
-0.2001070 
-0.2109032 
—0.2272509 
-0.2207007 
-0.2139103 
-0.1900390 
-0. 2057001 
-0.1939235 
-0.1717000 
— 0.1SC0029 
-O.I0700SO 
-0.1000700 
-0.1022331 
-0.1390120 
-0.1309220 
-0.1302103 
-0.1310370 
-0.1200957 
-0.1259270 
-0. 1231010 
-0.1203031 
-0.1 170050 
-0.1109900 
-0. 1123320 
-0.1090703 
-0. 1071005 
-0.1000330 
-0.1021093 
-0.0997299 
-0.0973770 
-0.0950990 
-0.0920506 
-0.0900052 
-0.0005052 
-0. 0000970 
-0.0000003 
-0.0025103 
-0.0000307 
-0.0707050 
-0. 0709709 
-0.0752101 
-0.0735303 
-0.0710770 
-0.0702592 
-0.0007190 
-0.0071000 
•0.0057257 
-0.0  002945 
-0.0020902 
-0.0015577 


OMEGA 
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530.00000 
S«0. 00000 
100.00000 
100.00000 
1T0. 00000 
110.00000 
590.00000 
600.00000 
010.00000 
010.00000 
030.00000 
040.00000 
050.00000 
060.00000 
670.00000 
610.00000 
090.00000 
TOO. 00000 
TIO. 00009 
710.00000 
730.00000 
740.00000 
750.00000 
700.00000 
770.00000 
700.00000 
790.00000 
500.00000 
010.00000 
010.00000 
130.00000 
040.00000 
030.00000 
000. 00000 
070.00000 
000.00000 
090.00000 
900.00000 
910.00000 
910.00000 
930.00000 
940.00000 
910.00000 
900.00000 
970.00000 
900.00000 
990.00000 
1000.00000 


1.0141590 
1.0137599 
1.1131051 
1.1111400 
1. ItllttO 
1. 1100000 
1.1100111 
1.1175411 
1.0101001 
0.1150000 
1. 0141939 
1.1141070 
0.1144499 
1.1147097 
1.C140099 
1.1149793 
1.0149733 
1.0140731 
1.0I400SO 
1.0144170 
1.0140031 
0.1136934 
1.0131011 
0.0117905 
1.0113317 
0.1110040 
0.01 15031 
1.1111014 
1.1111100 
0.11 11 000 
0.0111900 
1.0111959 
1.0111004 
0.0111090 
0.0110167 
1.0100940 
1.1107456 
1.1105003 
1.0103735 
1.0101610 
1.0199454 
1.0197307 
0.0193410 
1.0193017 
1.0191071 
1.1091911 
0.1191414 
1.1190997 


•0.0001175 
-0.0509705 
-0.0577116 
-0.0505161 
-0.0553410 
-0.0541143 
-0.0531030 
-0.0510144 
-0.0509195 
-0.0499450 
-0.0409034 
-0.0400091 
-0.0470741 
-0.0461730 
-0.0451017 
-0.0444191 
-0.0435003 
-0.0417374 
-0.0419053 
-0.0411091 
-0.0404009 
-0.0390957 
-0.0309011 
-0.0301037 
-0.0375914 
-0.0309375 
-0.0301014 
-0.0350375 
-0.0350391 
-0.0344396 
-0.0330395 
-0.0331077 
-0.0317117 
-0.0311773 
-0.0310311 
-0.0311114 
-0.0300104 
-0. 0301140 
-0.0190315 
-0.0191500 
-0.0107000 
-0.0101555 
-0*0170111 
-0.0173054 
-0.0109470 
-0.0100456 
-0*0101405 
-0.0157409 
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SAMPLE  PROBLEM  1 

LONG  PORN 

DAMPING  COEPP I CIENTaO. 50000000 

PREOICTCLES/SECI  MAI  RESPONSE  X  OMEGA  MIN  RESPONSE  X  OMEGA  MAX  RESIO  RESP  X  OMEGA 


10 .00000 

0.49S0034 

-0.1007304 

<0.00000 

0.01T3T30 

-0.1743734 

30.00000 

0.0314043 

-0.1379774 

AO. 00000 

0.01*1291 

-0.1730044 

SO. 00000 

0.03I04T4 

-0.1741004 

00.00000 

O.OSI64TO 

-0.1400190 

TO. 00000 

0.0409949 

-0.1449970 

00.00000 

O.OTS4T44 

-Ooimsae 

90.00000 

0.0746953 

-0.1157130 

100. 00000 

0.0T0I0S9 

-0.1051440 

1 10.00000 

0. 0490794 

1 to. 00000 

0.0590004 

-0.0943394 

130.00000 

0.0S496T4 

-0.0974307 

140.00000 

0.0513543 

-0.0900714 

ISO. 00000 

0.0490T10 

-0.0940404 

100.00000 

0.0476330 

-0o091«IM 

1 TO .00000 

0.0441553 

-0.0090399 

100.00000 

0.0447773 

-0.0077301 

190.00000 

0.0435200 

-0.0004470 

<00.00000 

0. 0434415 

-0.0034300 

<10.00000 

0.0414701 

MO. 00000 

0.0404375 

<30.00000 

0.0390413 

-0.0770310 

<40.00000 

0. 0391340 

-0.0740044 

<S0. 00000 

0.0304004 

-0.0743174 

<40.00000 

0.0374797 

-0.0734349 

*70.00000 

0.0349330 

-0.0709794 

<00.00000 

0.0341071 

-0.0493033 

<90.00000 

0.0303439 

-0.0470307 

300.00000 

0. 0340194 

-0.0443343 

310.00000 

0.0334019 

-0.0440033 

330.00000 

0.0347744 

-0.0434734 

330.00000 

0.0310033 

-0.0431003 

340.00000 

0.030990t 

-0.0407743 

3S0.00000 

0.0301349 

— 0.0S94771 

300.00000 

0.0393403 

-0.0043033 

3T0 .00000 

-0.0071000 

300.00000 

-0.0009090 

390.00000 

0.0347034 

-0.0040091 

400.00000 

0.0309419 

-0.0030074 

410.00000 

0.0301440 

-0.0037744 

430.00000 

0.0344413 

-0.0017040 

430.00000 

0.0337044 

-0.0007733 

440.00000 

0.0330031 

-0.0490104 

4S0. 00000 

0.0334443 

-0.0400744 

440.00000 

0.0310319 

-0.0479440 

4T0 .00000 

0.0313043 

-0.0470480 

400.00000 

0.0307000 

-0.0448190 

490.00000 

0.0301730 

000.00000 

0.0194700 

910.00000 

0.0191970 

•0.0437744 

S<0. 00000 

0.0107434 

-0.0489933 

50 
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sjp.oooeo 

0.0109129 

-9.0422442 

940*00000 

0.0170999 

-0*0413147 

950.00000 

0  oOl 75021 

-0*0407060 

500.00000 

0*01 71275 

-0.0491094 

970.00000 

0.0107601 

-0.0294272 

900.00000 

0.0164229 

-0*0307602 

900.00000 

0*01 69897 

-0;020t9l5 

000. 00000 

9.C1S7717 

-0.0275004 

OlO. 00000 

0.0194609 

-0*0360902 

020.00000 

9*0131769 

-0.0262021 

090.00000 

0.0140959 

-0.0257199 

040. 00000 

0.0146297 

-0.0251496 

090.00000 

0*010613 

-0*0346066 

000.00000 

OoOl 41 1 09 

-0.0940706 

070.00000 

0.0198701 

-0.0225417 

000.00000 

0. (196972 

-0.0920906 

090. 00000 

0.0194119 

-0.02259*9 

700.00000 

0*0131939 

-0.0220929 

710.00000 

0.0129027 

-0.0915794 

720.00000 

0*0127702 

-0.0911007 

790.00000 

0.0125801 

-0.0900620 

740.00000 

0.0129009 

-0.090  2299 

790.00000 

0.01 22C69 

-0.0297909 

700.00000 

0.0120279 

-0.0292649 

770.00000 

0.0110944 

700.00000 

0.01 I6B99 

700.00000 

0.0119221 

-0.0261691 

000.00000 

0.0119629 

-0.0277*4* 

010.00000 

0.0112000 

-4.0274071 

020.00000 

0.0110579 

-0.0270244 

090.00000 

0.0100109 

—0.4264472 

040.00000 

0.0107620 

-4.9269609 

090.00000 

0.0106912 

-4.029499* 

000.00000 

0.0104971 

-0.0296651 

070.00000 

0.0109664 

-4.429296* 

000.00000 

4.0102900 

-0.0290124 

090.00000 

0.010114* 

-0.92*70*4 

900.00000 

0.00*9991 

-0.0264022 

910.00000 

4. 0090746 

920.00000 

-4.0290160 

990.00000 

0.0096660 

-0.4294200 

940.00000 

0.0090940 

-4. 029*479 

990.00000 

0.0*94277 

940.00000 

4.0099222 

-4.4*170*9 

970.00000 

-4.4*2*60* 

400.00000 

0.009I194 

-4.0*2100* 

490.00000 

2.0294109 

-0.0219264 

1000. 90000 

9. 4944244 

-4.421*429 

22062*9  9 AN  70  COMMIT  ION 
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APPENDIX  B 

SAMPLE  PROBLEM  2 


Arrangement  of  Input  Data 
and 

All  the  Various  Forms  of  Output 
That  Can  Be  Generated 
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SAMPLE  PROBLEM  2 


The  data  Z  for  this  problem  was  furnished  by  George  O’Hara  of  the  Naval  Research  Laboratory. 
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ENO  0F  FILE 


SRESTORE 


JIBSYS 
*  • 


ENO  OF  FILE 


(SYSTEM  CARO)  \ 


(SYSTEM  CARO) 


(SYSTEM  CARD)  \ 


(SYSTEM  CARO)  \ 


(BLANK  CARD) 


10.  1000.  10.  (FREQUENCY  CARD) 


ISO.  1.  (FREQUENCY  CARD) 


SZAZA4  ALPHA(l)  =  0,  .1,  NALPHA  =  25  (4TH  NAMELIST) 


,  IP2  «=  0,  IP3  =  1,  IP4  =  1,  IP5  -  IS  (3RD  NAMELIST) 


,  TZERO  =  .145,  H  -  .12464165E  -  03.  XO  -  0,  N  - 1164,  XDOT  0-  0,  ISZ  -  OS  (2ND  NAMELIST) 


-0072,  -00571  -0057,  -0050S 


-0050,  -0036, 


SZAZA1  Z 


SDATA 


(1ST  NAMELIST) 


(SYSTEM  CARD) 


i  i  i  •  j  i  !  i  1 1* »  » ii *  »*  ii  « k iai< nnxnn nnif » n hum n  im 
II,  1,  1  III  II  III  II  III  1 II II  M 1 II  III  1  II  II  1  III  1 1 1  III  1  III  II  1  1 II II 1  111  II  III  1 II  III  III  1 II 

7  7  7  7 

77777777777 

7  7  7  7  7 

7  7  7  7  7 

7777777777 

7777777777 

7777777777777777777777177 

7  7  7  7  7  7  7  m 

(ill,. 

1117)31)}]] 

lllll 

))))} 

1)1)111))) 

13)3)13))) 

))))))))))))))))))))))))) 

)))))))))) 

44444444444 

4  4  4  4  4 

4  4  4  4  4 

4444444114 

4444444444 

4444444444444444444444444 

4444444444 

5  5  5  5 

55555555555 

ssm 

ssm 

SSSSSSSSS5 

sssssmss 

smms5smsmmsssss 

sssssmss 

UK 

lllllllllll 

Kill 

IKK 

KIIIKIII 

KKIIIKKIIIKIIKIKK 

KIIIKIII 

7  7  7  7 

77777777777 

7  7  7  7  7 

7  7  7  7  7 

7777777777 

7777777777 

7777777777777777777777777 

III 

lllllllllll 

Kill 

IKK 

KIIIKIII 

unit  mi 

lllllllllllllllltllllllll 

llllllllll 

ms 

'in 

9999999999999999999999999999999999999999999999999999999999999999999999999999 

SIM  t  IM17  0H  •  •  1?  IUU1  NTS 2177* *»JI  JM)  )«  »  N  17  W  »«« IM|  «)«  «t«  «?«#•»  >7  UU  >4  »  M  «l»  »M  It  MU  M  It  N  »  M  M  ■  >1  1}  7)  «  II  *  H  71  VI  •» 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2 


SNERO  -  EFN  SOURCE  STATEMENT  -  IPN(S) 


06/10/66 


0 (MENS ION  Z  < 2000  > . T  <  2000  > . S 1 2000  > . S2nO ( 2000  > . XOMEGA ( 2000  >  * X ( 2000 ) • 
IXOOTt 2000 1 • ALPHA! 23 1 .OSCILt 1000)  .VeH 1000 ) .RESIOI 1000 1 
EQUIVALENCE  (RESIO.Xl 1 1  ) . < VEL #X < 1001 > > 

NAMEL I ST/Z AZ A 1/Z/Z A2 A2/T #  TZERG •  H*  XO.N, 

1 XOOTO • I SZ/Z AZ  A3/ IP1#IP2#IPS#I P4 • I P5/ZA  ZA4/ ALPHA •  N ALPHA 
MR  I TE ( 6#900 ) 

PI-3. 1418927 
REA0C9.ZAZA1) 

REAO( 9*  ZAZA2 ) 

REAOI 9#ZAZA3) 

REAO( StZAZAA  I 
IF( IPS ) 50 • 31 • 50 
SO  IP1-0 
IP2-0 

81  CONTINUE 
00  75  1-1. N 

78  Zm-tZ<I)*80.)M.0139S9 

c 

C - - - — WRITE  ANO  PLOT  INPUT  OATA - - 

00  100  1>2#N 
100  TI I >-TI !-!>♦* 

IPC ISZ 1 104# 105# 104 
104  WRITE<6.901> 

GO  TO  106 
108  WR1TEI6.902) 

106  CONTINUE 

MR1TEI6. 903 )  (Z( I ) #TC I ) . I-I.N) 

ZMIN-Ztl) 

ZNAX-Z(I) 

00  110  I-2.N 

IP I ZMAX-ZI I ) )100tl09tl09 

ioo  zmax-z(i) 

60  TO  110 

100  IPIZNIN-ZdMUO.UO.lOT 

107  ZNIN-ZII) 

110  CONTINUE 

CALL  CAMRAVC  38 ) 

ZN-N 

0X-HMZN/10# 

0V-(ZNAX-ZM1N)/10. 

CALL  GNI01V(l#T(l)tT(N) tZMIN# ZMAXtOXtOVtl tltltlt6#6) 

CALL  APLOTVIN.T.Z. 1 1 1 • 1 #42 1 1 ERR ) 

1P1 IERR) 118*116# 118 
110  «RITE(6#904)  IERR 
110  CONTINUE 
NNaM-1 

00  120  I-l.NN 

CALL  L1NEVINXVITI  m.NVVIZII)  )#NXVITII41))»NVV(Z(I4t  )>> 

120  CONTINUE 

CALL  PRINTVI-IS.ISHTINE  IN  SEC0N0S#482 #6) 

IP! ISZ >121.122 #121 

121  CALL  APRNTVIO#*14#-23»23HPOUNOATION  ACCELERATION. 4. 696) 

60  TO  123 

122  CALL  APRNTV(0#”14#”19# 19HPOUNOATION  VELOClTV.4#664) 

123  CONTINUE 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 

SNERO  -  EFN  SOURCE  STATEMENT  -  JEMS) 


06/10/66 


C- — - COMPUTE  S(  M  AND  ScND(N)  FOR  EACH  TIME.T - 

S(l)*2(2)-Z(l> 

S2NO (  l)*Z(3)-2.*Z(2|+Z( 1  I 
K«N-1 

00  200  1*2. K 
S( I )*Z( !♦! *  — Z  < I 1 

200  S2NO(tl*2( !♦! )-2.*Z( I )*Z(  1-1 1 
WRITE! 6. 917)  (  l.S( I ) . S2ND (  I  1 . 1*1 .N 1 

C 

190  READ( 5.920 )  FREQ1.FRE02 . OELTAF 

IF! AB5IFR6Q1-FREQ21-. IE- 20)  801.001.191 

191  CONTINUE 

IF( IP3I192.19S.192 

192  TEMP*FRE02-FRE01 

IF I AM00I TEMP  .OELTAF  I— DEL TAF/2. 1193.194.194 

193  N0FRE0*TEMP /OELTAF 
60  TO  199 

194  N0FRE0«TEMP/0ELTAF4 1.0001 
199  CONTINUE 

C— — OANPEO  OR  UN0AMPE0.1F  UNOAMPEO  FINO  NUMBER  OF  RESIOUAL-  POINTS— 
00  800  IX*1 .NALPHA 
N  IN*  1 
IN  1*0 

FREQ*FREO 1 

IFIABSI ALPHA! IX 1 1-. IE-09 1 201 • 20 1 • 202 

201  1 ALPHA* 1 
60  TO  222 

202  1ALPHA-2 

222  60  TO  (203.2071. IALPHA 

203  TT*T2ER0-T( 1 1 

IF( AMOOI TT .HI-H/2. 1209.204.204 

204  KH*TT/H41. 000001 
60  TO  207 

209  KH*TT/H4. 000001 
207  CONTINUE 
C 

C - COMPUTE  CONSTANTS—— 

0ME6A*FRE0*2.*P I 

R ACC AL* SORT ( | .—ALPHA ( IX)«*2I 

OH*OMEGA*H 

P*OME6A*HAOCAL 

E  X»EXP  C  -  ALPHA  (  IXMOHl 

COSPH*COS(P«H) 

S 1NPH*S IN(P*H| 

A20H*2. * ALPHA ( IXl/OH 
REC0H*1./0H 
ALPHA2* ALPHA (  IX)«*2 
X 1 2* 1 a~2.*ALPHA2 
F ACT 1*EX*S INPH/R AOCAL 
F ACT2*EX*C0SPH4FACT l*ALPHA( IX 1 
IF( ABS( FREO )-« IE-20  1220 .220.221 

220  IFREQ*! 

60  TO  223 

221  IPREQ*2 

223  IFf 192)300.208.300 
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PROGRAM  LISTING  FOR  SAMPLE  PR03LEM  2  (Continued) 


06/10/66 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 

200  00  TO  (400*209 ) • IFREO 

209  CONTINUE 
C 

C— — — COMPUTE XONEGM N )  (XDOT  FOR  VELOCITY  INPUT - - - - — 

FACT3*REC0H«( 1.-FACT2) 

F ACT 4*REC0H« I .5-2.4 ALPHA! I  X )«REC0H4EX4 I ( . S42.«ALPhA! IX) 4RECOH) 

1 4COSPH-1 X 1 24REC0H-ALPHA (IX)/2*)*SI NPH/RAOCAL ) > 

FACT5*EX*C0SPH- ALPHA! IX)4FACT1 
FACT6*REC0H*FACT1 

FACT 7*1 RECOH**2  )•( l.—EX4(C0SPH4( ALPHA! IX140H/2. t*S I NPH/RAOCAL) > 

GO  TO (224 *2 25) *NIN 
224  XOMEGA! 1 )*OMEGA4XO 
XDOT! 1 )*XOOTO 
00  210  1*2. K 

XOMEGA! t l*XONEGA( 1-1 1*FACT2 4X00 T ( I— 1 I*FACT1-S( I- 1 > *FACT3-S2nO ( 1-1 ) 
14FACT4 

210  XOOT ( 1 l«-X0ME6A( 1-1 )*FACT14XOOT ( 1-1 )*FACT5-S( I— 1 1*FACT6— S2N0! 1-1)4 
1FACTT 

GO  TO  756 

229  CONTINUE 
KR0*KMX4k00K 
KRA*KMN4KICK 
00  235  1*2*5 
VA»I— 1 
TA«yA4«2*H 
COSPT«COS(P»TA) 

S1NPT*S 1N(P*TA ) 

EXT* EXP (  —  ALPHA (  IX  )40MEGA*TA ) 

FACTR1*EXT*(C0SPT4 ALPHA!  I  X ) *S INPT/RADC AL ) 

FACTR2*EXT  *S INPT/RAOCAL 
FACTR3-! 1 •— FACTR1 )*RECOH 

FACTR4* ( T A/H- . 5— A20H4  (  ALPHA  (lX)/2.-X12«  PECOH )*FACTR24| . 54A20H) 

1 *EXT  *COSPT ) 4REC0H 

XMAXO  —XOMEGA (KRO )*FACTR 14X00  T ( KRO ) *FAC TR2- S ( KRO ) *FACTR3 

1— S2N0(KR0 )«FACTR4 
IF ( XM AX— XMAXO  1229.230.230 

220  XMAX*XMAXO 

230  XMINO  “XOMEGA  CKR A )»FACTR 1 4X00T! KRA)4FACTR2-S(KRA )«FACTR3 

1-S2N0(KRA)«FACTR4 

IF (XMIN— XMINO  1235.235.232 

232  XM|N*XMINO 
239  CONTINUE 
NIN*1 

211  IF! IP1 1900.212 .500 

212  IF! IP2IS90.213.fiSO 

213  IF! IP3I600.214.600 

214  IF! IP4)6S0. 215.650 
219  IF! 1P9I7OO.705.7OO 

C 

C - COMPUTE  XOMEGAIN).  XOOT(N).  FOR  ACCELRATION  INPUT - 

300  CONTINUE 

GO  TO! 350 .305). IFREO 
309  CONTINUE 

EC0Sl*l.-eX4C0SPH 
0FACT3* 1 1 .—FACT  2 ) /OMEGA 

0FACT4*! 1 .— A20H«EC0S1-X124FACT1 4REC0H) /OMEGA 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


SNERC  -  EFN  SOURCE  STATEMENT  -  IFMS1 


06/10/66 


0FACT5*(-4.*ALPHA<  IXI-(2.»(  1 . -4 . • ALPHA 2 1 /OM-2. • ALPHA ( (X) )*EC0S1+ 

1  (XI2+2.* ALPHA' IX>M3.-4.«ALPHA2  )/0H)«FACTI  ) • REC OM/ ( 2 . *OMEG A ) 
DFACT6M  1.-EC0S1)-FACT1»ALPHA<  1  X) 

0FACT7>FACT1 /OMEGA 

OF ACTS* ( 1 (—FACT 2 ) • RE COH /OMEGA 

OF  ACT9*  (  2. ALPHA  (  I  X  )  ARECOH  »  »ECO  S 1  -  <  2.  •  X 1 2»RFC0H- ALPHA!  I  X  1  1 
t«FACTl )*REC0H/(2.*0MEGA| 

GO  T0C3C8.32S)  .MN 
308  XOMfcGA( 1 )*OM£GA*XC 
XCOT (  n=XOOTO 
CO  3 1C  1*2. K 

XOMfcGAJ  I »*XOMEGA(  |~ 1 > •»  AC  1 2.XCC ! 1  I - 1 ) *F AC T 1 - Z (I- 11  «0F,  CT3-SI 1-l» 
MCFACT4-S2NCI  l-l  )»UHACTS 

310  XOO  T  (  I  1  =  —  XOMEO  A  (  I-I)«FACT1«XOOT(  1  -  1  )  »r>F  AC  T6- 2  (I  - 1  » «0FACT7 
1  -SC  l-l  )»OACTM-S2NO<  I- I  1»CFACT<J 

GC  TO  756 
325  CONTINUE 

KRO=KMX.KOOK 

KRA*KMN»K ICK 

H2*H«»2 

00  335  1*2.5 

VA*I-l 

TA«VA*.2«H 

COSPT»COSlP*TA  > 

S1NPT*SIN(P*TA) 

EXT* EXP ( —ALPHA (  IX  >*OMEGA*T A ) 

FACTRlxEXT * I COSPT4 ALPHA{  I X I *S INPT/RAOC AL 1 

FACTR2*EXT«SINPT/RA0CAL 

FACTR3*! 1.-FACTR1 l/OMEGA 

F  ACTR4* I T  A/H— A2GM* ( 1 .-EX T*COSPT 1-X1 2*F ACTR2OREC0H > /OMEGA 
FACT RS* ( T A« • 2/H2-T A/H- 1  2 . • ( 1 . - 4 . • ALPHA  2 1 *REC0H**2- A20HJ • 

II 1 •— EXT*COSPT )  ♦ ( X 1 2»REC OH+2 • • ALPHA  I IX) ■  1 3 .-4 . * ALPHA2 1 *RECOH**2 ) 
24FACTR2I/I2.40MEGA) 

XMAXO*X OMEGA I KRO ) *F  ACT R 1  AX DOT  I KRO I *F AC  TR2-Z ( KRO) *FACTR3 
1  — St KRO ) *F ACTR4— S2N0 ( KRO )*F ACTRS 
IFIXMAX-XMAXO  1329.330.330 

329  XMAX«XMAXO 

330  XMIN0*X0MEGAIKRA1*FACTR14X00TIIIRA14FACTR2-ZIKRA)*FACTR3 
l-SIKRA)*FACTR4-S2N0(KRA)*FACTR5 

IF  <  XM IN— XM | NO  1335.335.332 

332  XMIN«XMINO 
335  CONTINUE 
NIN*1 
GO  TO  211 
C 

C - FREOl'ZERO  FOR  FOUNDATION  ACCELERATION - 

350  GO  TO  1 352  *  365 1  .'N  IN 
352  XOMEGAI 1 >*X0 
XOOT ( 1 l*XOOTO 
H2«H**2 
00  360  I*2.N 

XOMEGAI I )*XOMEGA| I-1)*H»X00T( I-t )-Z( I- 1 »*H2/2.-S< 1-1 >»H2/6«4 
1S2N0I 1-1 14H2/24. 

360  XOOT 1 1 t*-ZI I-l1*H-SI  l-l  l*M/2. 4S2NDI 1-1 l»H/12. 

GC  TO  756 
365  CONTINUE 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


06/10/66 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFNI  SI 


N  IN*  1 
GO  TO  211 
C 

C- - —COMPUTE  X.  XOOT  FOR  OMEGA*ZERO - VELOCITY  INPUT 

400  GO  TO ( 40 2 • 425 )  »NIN 
402  XOMEGAI 1 l*X0 
XOOT I  1 )*X0QT0 
00  410  1*2. N 

XOMEGAI  I  )*XOMEGA<  I-D4X00TI  I- 1 »  »H-SI  I- 1 ) «H*. S4S2NDI 1-1  ) «H/12. 
410  XOOTII)*XOOTtI-l>-StI-l) 

GO  TO  756 
425  CONTINUE 
N  IN*  1 
GO  TO  211 
C 

C— — - -PLOT  X  vs.  T  ANO  XOOT  VS.  T - 

500  GO  T0I50S.S10) . IFREQ 

505  00  506  1*1. N 

506  XI I )*XOMEGAI 1 1 
GO  TO  SIS 

510  00  515  1*1. N 

SIS  XI 1>*X0MEGA( I l/OMEGA 

SIS  XNX*Xtl) 

XOMX-XDOT ( 1) 

XMN*Xt 1 > 

XOMN*XOOT 1 1 1 
00  530  1*2. N 
1PIXI11-XMX  >521.521.520 

520  XMXaXI I > 

521  IPIXIII-XMN  >522.523.523 

522  XNN-XIH 

523  IPIXOOTI  D-XOMX  >525.525.524 

524  XONX-XOOTU) 

525  IPIXOOTI I l-XOMN  >526.530.530 

526  XCMN«XOOTI I > 

530  CONTINUE 
OVatXMX  — XMN  )/ZN 

CALL  GRIOIVU.TIU.TINI.XMN.  XMX.  0X.0V.4.4.4.4.6.6I 
CALL  APLOTVIN. T *X. 1 . 1 • 1.42.  IER> 

IPI 1ER I S3 1.532.531 

531  MRITEI6.913)  IER.FREO 

532  NN*N-1 

00  535  1*1. NN 

CALL  LINEVINXVI T( l)).NVVIX(  I) > .NXV1TI I +1 >  >.NVV(X( 141) >> 

535  CONTINUE 

CALL  PRINTVI  — 15.  15HTIME  IN  SECONOS.482  .6> 

CALL  PRINTVI-23.23HPRE0*  CVCLES/SEC. .20.3) 

CALL  LAOLVIPREO. 60.3.6. 1.4) 

CALL  PRINTVI-8.0HALPHA*  .20.17) 

CALL  LABLWIALPHAI IX>.72. 17.6. 1.1) 

CALL  APRNTVIO.— 14.-1 1 • 1 1HX  RESPONSE .4 .600) 

OVbIXOMX  -XOMN  )/ZN 

CALL  GRI01VI 1. Til). TIN). XOMN.  XOMX.  0X.0Y.4>4.4.4.6.4) 

CALL  APLOTV  IN. T, XOOT.  1 • 1  *  1 .42 . IERA) 

1PI1ERA)  535.539.536 
530  MRITEI6.914)  IERA. FREQ 
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Pfl  /GRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


06/10/6C 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFNIS)  - 

839  00  840  I  * 1 , NN 

CALL  LINEVINXVlTI 1 ) ) , NYV(  XOOT ( I ) ) ,NXV<  T I 1 ♦ l ) ) , N YV I XOOT ( J ♦ l ) ) ) 

840  CONTINUE 

CALL  PRINTVt-lS.lSHTIME  IN  SECOND St 452. 61 
CALL  PR  INTV I  ”23  •  2 3HFREQ*  C VCLES/SEC . *20 .3) 

CALL  LABLVt FREQ.60.3.6,  1.4) 

CALL  PRINTV(-8.8HALPHAs  .20.17) 

call  lablvialphai ix). 72.17.6.1.1) 

CALL  APRNTV 10,-14  «-22, 22HVEL0CI TV  RESPONSE  XOOT. 4.688) 

CO  TO  212 

- - PLOT  VI ACCELERATION)  VS.  TIME - 

880  CONTINUE 
GO  T0I213.S81). IFREO 

881  CONTINUE 
860  00  862  1*1.  N 

862  X  C I )a—X OMEGA ( I ) *ONEGA 
XMXaXI 1 ) 

XMNaXI 1 ) 

00  866  1*2. N 
1PIXMX  -XI 1 )  1863. 864. 864 

863  XNX*XII> 

864  1PIXMN  -XI  I )  1866.866.868 
868  XMNsXI  I  ) 

866  CONTINUE 

OVbIXMX  -XMN  )/2N 

CALL  GRI01VI  l.T(l).TIN). XMN,  XMX.  OX .D V .4 .4  .4  .4 .6 .6) 

CALL  APLOTVCN.T.X. 1.1. 1,42. IERB I 
IF! IER8 >870.871 .870 

870  HRITE<6.91«)  IERB. FREQ 

871  NNBN-1 
00  878  1*1 .NN 

CALL  LINEVINXVITI  I) >.NVV(X( I) I.NXVITf I  41 ) ) .NVV I X ( I *1 > >  ) 

878  CONTINUE 

CALL  PR  INTV I  — 18. I SHT IME  IN  SECONOS .482 .6) 

CALL  PRINTVI— 23.23HFREQ*  CYCLES/ SEC. .20.3) 

CALL  LABLV! FRE0.60, 3.6. 1.4) 

CALL  PR INTV ! —8. 6H ALPHA*  .20.17) 

CALL  LABLVIALPHAI IX), 72. 17, 6.1.1) 

CALL  APRNTV 1 0 .-14 .-26* 26HABS0LUTE  OR  Y  ACCELERATION. 4. 720) 

GO  TO  214 
C 
C 

C- — - - -FOUR  CC0R01NATE  GRIO— - 

600  CONTINUE 
INI»14INI 

OSC ILl INI )«FREO 

VELI INI ) a AM AX 1 I A8S I XMAX ) . ABSI XM IN ) ) 

GO  TO  1601. 602). IALPHA 

601  RESIOt INI ) a  SORT I XOMEGAIKH) **24XDOTIKH>  *42) 

IF  IRES  101  INI)-. IE-19 >604,604, 602 

602  IFIFREQ-. IE— 10)604,604,603 

603  IFIVELI INI  )-. IE- 19)604.604.603 

604  NOFREQaNOFREO-1 
INIaINl-1 

608  IF! ABSI FREQ— FREQ 2 )-• IE-091606.606.214 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFN(S) 


06/10/66 


606  CONTINUE 

CALL  SMXVVtl.l) 

FRQSML*OSCIL(  1  I 

VELSML*VEL< 1  I 

RESSML*RES tO( 1  I 

00  616  I*2.N0FRE0 

1FI OSC IL ( I  l-FROSML  161 1.612*61  2 

611  FROSML*OSCIL( 1  ) 

612  IF1VELC  I  >-VELSML >613.614,614 

613  VELSML*VEL( 1 > 

614  IFIRESIOI  I  l-RESSMD61S.61S.616 

615  RESSML*RESIO( I) 

616  CONTINUE 

OROSML*AMIN1 ( VELSML • RE S SML  ) 

c 

C - - TO  FIND  THE  LIMITS  FOR  THE  <4X31  CYCLES - 

C 

IF <  AL0G10 ( OROSML ) 1620.621  .621 

620  LOGURO* ALOG 10 ( OROSML 1-1 . 

GO  TO  622 

621  LOGORO* ALOG 10 <  OROSML 1 

622  ,!F<  ALOG 10 < FRQSML  1 1623.624.624 

623  LOGFRO*ALOG 10 ( FRQSML  I- 1 . 

GC  TO  629 

624  L0GFR0>AL0G10(FR0SML) 

629  OROSML* 10  <  * •LOGORO 

OROLG«OROSML* 10 .**4 
FRQSML* 10. ••LOGFRO 
FR0LG*FRQSNLM0.««3 

C- — — PLOT  LOG-LOG  GRIO  FOR  VELOCITY  VS  FREQUENCY — - - 

CALL  GR 101 V< 1. FRQSML. FRQLG • OROSML .OROLG. 1.0. 1.0.1. 1.1.1. -2.-2) 

C - TO  FIND  LARGEST  1X10  TO  THE  PTH  POMER  LINE  FOR  X— ----- 

0MSTRT*2.«PI«FR0SML 
TEMP*  ALOG 10 ( ORDLG/OMSTRT 1 
IF<  TEMP  1627 .628.628 
627  LOPaTEMP-1. 

GO  TO  629 
626  LOP«TEMP 

629  01*10. ••LOP 
V1*0MSTRT«01 

• I G* 10 ••OROSML 

XNARGN-NXVI FRQLG l-NXVI FRQSML) 

VMARGN>NVV( OROLG >-NYV<aIG) 

SLOPE* YMARGN/XMARGN 
I XC0R*4 ANX V <  FRQSML  > 

CALL  PRINTVI-l. 1H0. IXCOR.NYVI  VI)) 

CALL  PRINTVI -2»2HO*.900.  IT  I 
CALL  LABLVIO 1.924. 17,-2.  1.3) 

C - CRAM  LINES  UP  FROM  VI - — - 

VIS*0. 

DO  632  1*1*9 
V1S-VIS*V1 

IFIOROLG-V IS >647.630, 630 

630  XTCH-NVVI OROLG )-NVV< VIS) 

1TCH*XTCH/ SLOPE 

I TCH*NXV< FRQSML  )4ITCH 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


06/10S6I 


SNERD  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 

KIT  CHVNYV  (  OROLG ) 

CALL  LINEV  (NKV«FRQSML).NVV(VIS).ITCM,KITCH» 

632  CONTINUE 

C - .-.-ORAM  LINES  FROM  VI  OOMN - 

647  DO  636  J*1.20 
LAL*LOP-J 
LAN-LAL41 
CLAN«10.**LAN 
0LAL*10.««LAL 
00  638  1*1.9 
EFTI* I 

VIS*0MSTRT*CLAN-0MSTRT»0LAL»EFT1 
IF! V IS-ORDSML 1639 *633 *633 

633  KITCH»NYV(OROLG>-NVV(  VIS) 

I VMARGaYMARGN 

IFI IYMARG-K I TCH >635 • 634 . 634 

634  2 ITCH*K ITCH 
ITCH*Z ITCH/SLOPE 
GO  TO  637 

639  ITCH-XMARGN 
KITCH*VMARGN 

637  ITCH* I T CH+NXV I FROSML  ) 

K  ITCH*K  ITCH4NYV! V  IS) 

CALL  LINEV  (NXVIFROSMLI.NYVIVISI.ITCH.KITCH) 

636  CONTINUE 

C— — ~ ORAM  FROM  THE  RIGHT  ORDINATE  THE  REMAINING  X-LINES - — 

639  0NLAST*2.«PI«FR0LG 
LEFT2*J 
LEFTI-EFTi 

00  641  J-LEFT2.30 

LAL*LOP-J 

LAN*LAL4l 

0LAN*10.««LAN 

DLAL*10.**L AL 

00  640  l*LEFTt.9 

EFTI*I 

VIS«OMLAST*OLAN-OMLAST«OLAL»EFTI 
IFI VIS-OROSML  >642.636.636 
636  K I TCH»NYVI V I  SI— NY V I OROSML ) 

CAPT-KITCH 
ITCHaCAPT/SLOPE 
ITCHaNXVt  FROLC )- ITCH 
K  |  TCHaNVV I OROSML  > 

CALL  LINEV  | NXV I FRQLG ) . NYV C  V IS).  ITCH.  KITCH) 

640  CONTINUE 
LEFTlal 

641  CONTINUE 

642  CONTINUE 

C - - - PLOT  POINTS- — - 

CALL  APL0TV|N0FRE0.0SC|L.VEL.l.l.t.44.LIT> 

IFILITI643.644.643 

643  MRITEI6.922)  L IT . ALPHA! IX ) 

644  60  TO  (640.646). IALPHA 

640  CALL  APLOTVINOFREO.OSC IL.RESID. 1.1.1.38.KIT) 

CALL  PR INTV I —26. 20HC00  -RESIDUAL  SHOCK  SPECTRUM. NXVIFROSML) .1006) 
IF! K IT  1640.646.648 
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SNERO  -  EFN  SOURCE  STATEMENT  -  IFMS) 
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645  NRITE16.923)  KIT. ALPHA! IX 1 

646  CONTINUE 

CALL  PRINfV(-19« 15HFRE0UENCV  (CPS). 452.61 
CALL  APRNTV10.— 14.— 0.6HVELOCI TV. 9. 576) 

CALL  PR INTW ( —6. 6H ALPHA*  .20.1T) 

CALL  LABLV(ALPHA( IX ) . 72 . 1 7.6. 1 . 1 ) 

CALL  PR INTV(  — 19  » 1 9H***  «SMOCK  SPECTRUM. NXV(FROSML) .1015) 


C — - — - GENERATION  OF  ACCELERATION  GRID 

C-— — — TO  PINO  AI - 


TEMP*AL0G10(0R0SML«0MSTRT) 

if (Temp) 650 .651 .esi 

650  LOPaTEMP 
GO  TO  653 

651  LOP*TEMP* 1 . 

653  A|aiO.**LOP 

VlaAl/OMSTRT 

CALL  PR INTV(-I. IH A » I XCGR .NVV( VI ) ) 

CALL  PRINTV(-2.2hAa, 900.4) 

CALL  LABLV( Al. 924. 4. —2.1.3) 

C— — — OR A«  LINES  OOMN  FROM  Al - 

K I TCN*NVV ( OROSML ) 

00  656  I >1 .9 
EFT la  I— l 

V1S«V1«( l.-.t«EFTl) 

IF( OROSML— VIS >655.655.657 

655  XTCHaNYV (VI S )— NVV ( OROSML ) 

ITCHaXTCH/SLOPE 

I TCHa ITCH+NXV (FROSML) 

CALL  LINEV  ( NXV ( FROSML ) *NYV( V  IS ) • ITCH.  KITCH) 

656  CONTINUE 

C - ORAM  LINES  FROM  Al  UP - 

657  00  670  J»1 *20 
LANaLOP* J— 1 
ALANa  10  •  >*L  AN 
00  670  IaZ.10 

EFTlal 

VISaALAN/OMSTRT«EFTl 
IF(VlS-0R0LG)6Se, 656.671 
656  KITCH-NVV(VIS)-NVV(OROSML) 

I YM ARG« VM ARGN 

IF( IVMARG-K ITCH  >662.661 <  661 

661  ZITCHaKITCH 

I TCM»Z I TCH/SLOPE 
K I TCHaNVV ( OROSML ) 

GO  TO  066 

662  ITCHaXMARGN 

K | TCHaNVV (VIS)— IVMARG 
660  ITCH* ITCH*NXV(FRQSML ) 

CALL  LINEV  ( NXV( FROSML ) .NYV( V IS ) • I TCH.  KITCH) 

670  CONTINUE 

C— — — ORAM  FROM  THE  RIGHT  0R01NATE  .REMAINING  A-lInES< 

671  LEPTEaJ 
LEFTlaEFTl 

K I TCHaNVV (OROLG ) 

00  675  J«LEFT2. 30 
L ANaLOP* J— 1 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


SNERO  -  fcFN  SOURCE  STATEMENT  -  IFN(S) 


06/10/6< 


ALAN*10.««LAN 
00  874  1 *LCFT 1 • 10 
EFT  1*1 

vis*alan/omlast*efti 

IF(VIS-0RDLG)872,B72,876 
872  CAPT*NYV(OROLG)-NVV(VIS> 

I tch»c apt/slope 

I  TCH*NX V ( frolg > - I tch 

CALL  L INEV  ( NXV( FROLC > ,NYV ( VI S > ,  ITCH.KITCH) 
874  CONTINUE 
LEFT  1*2 
879  CONTINUE 
876  CONTINUE 
GO  TO  214 


SHORT  FORM - 

690  CONTINUE 

IF! ABSIFREO-FREOl 1E-04I6S1  ,651,692 

691  NRITE(8,912)  ALPHA! IX) 

692  GO  T0!6S8,699),  I  ALPHA 

699  XOM AX* AM AX  1 ! A6S ( XMAX ),ABS(XMIN)  ) 

GO  TO (668 ,669  )  *  IFREO 

698  RM AX* SORT  I XOMEGA! KH  >**24X00  T ( KH ) **2 ) 

XOM  AX*  AM  AX  1  (  A0S  (  XMAX  )  ,A8S!  XM1N)  ,  ABS  (  RM  AX  )  ) 

GO  T0(681,678).1FREQ 

668  NRITE(8.906>  FREO.XOMAX 
GO  TO  219 

669  NR  I TE( 8,909  >  FREO.XOMAX 
GO  TO  215 

678  MRITE(8,907>  FREO.XOMAX ,RMAx 
GO  TO  219 

681  MR  I TE (8,908  )  FREO.XOMAX.RMAX 
GO  TO  215 
C 

C - — — — LONG  FORM - 

700  CONTINUE 

IF( ABS(FREQ-FREQl)-. IE-04 ) 70 1,701, 702 

701  MR ITE(  6,915)  ALPHA!  IX) 

702  GO  70(704,705), IALPHA 

704  RMAX*SQRT  (  XOMEGA  ( KH  )  »*2  AXDOT  (  KH  )  •  *2  > 

GO  70(743.742), IFREO 

705  GO  TO(726. 727), IFREO 

726  NRITE (6,916)  FREO • XMAX, XM IN 
GC  TO  789 

727  MRITE(6.910>  FREQ, XMAX, XMIN 
GO  TO  785 

742  MRI TE(6«910  )  FREO, XMAX, XMIN, RMAX 
GO  TO  785 

743  MRITE(6,91 1  )  FREO, XMAX, XMIN. RMAX 
GO  TO  785 

C 

C— — — PRELIMINARY  CALCULATIONS  FOR  MIN.  AND  MAX.  RESPONSE 

796  GO  TO  (758. 757). I ALPHA 

797  NX*N 

GO  TO  759 
788  NX*KH 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


06/10/66 


SNFRD  -  EFN  SOURCE  STATEMENT  -  IFN(S>  - 

759  XMAX=XOMEGA< 1 1 
K*X  =  1 

XKlN=XOMECA(  1  I 

KMN=  I 

Cd  764'  1=2. NX 

IF( XMAX  —  XOMEGA (  I  )  1761.762. 762 

761  XMAX=XOMfcGA<  I  I 

KMX=  I 

762  IFIXMIN— XQMEGA1  11)764.765.765 
765  XMlN=XOMEGAI  I  ) 

KMN=  I 

760  CONTINUE 

IF1KMX-1 1771, 77., 771 

770  KMXaKMX4l 

771  IF1KMX-NXI773.772.773 

772  KMX«KMX-1 

773  IF1KMN-1 1775.774,775 

774  KMN=KMN4l 

775  IF1KMN-NX1777.776.777 

776  KMN=KMN-1 

777  CONTINUE 

IF1X00TIKMX-I 14X00T1KMX1 1776.T78.779 
776  KCOKa-l 

CC  TO  781 
779  K00K=0 
701  CONTINUE 

IFIX00T1KMN-1 )«XOOT(KMNI 1763.703.764 

763  KICKa-l 
GO  TO  950 

764  KICKaO 

950  N IN«2 

00  TO  227 

C 

c - —EITHER  RESTART  WITH  NEW  FREO.  OR  RESTART  WITH  NEW  ALPHA  OR  STOP 

70S  CONTINUE 

IF!  ABSl  FREQ-FRE021-.  lE-G3)6«-0.7!5,755 
755  FREOaFREO  »OELT  AF 
00  TO  207 

600  CONTINUE 
00  TO  190 

601  WR1TE16.9091 
CALL  PRAMEV 
STOP 

C 

900  FORMAT! 1H1.49X.22HSH0CK  SPECTRUM  PROOR AM//S0 X . 21H APPLI ED  MATH  LAB 
I0TM6.///1 

901  FORMAT 140X.36HFCUNDAT ION  ACCELERATION  VS  TIME! INPUT  1 // 1 BX.23HF0UN0 
1  AT  ION  ACCELERAT  I0N.36X, 15HTIME  IN  SECONOS1 

962  FORMAT!  42X.35HFCUN0ATI0N  VELOCITY  VS  TIME  1  INPUT 1//20X  « I9HF0UN0ATI 
ION  VELOCITY, 3SX.1SHTIME  IN  SECONDS! 

903  FORMAT! IX, F35.9.  44X.F20.9I 

904  FORMAT! 47H  PLOTTING  ERROR  NUMBER  OF  POINTS  OUT  OF  RANGES. 12//! 

905  FORMAT! 15X.F 10.5, 29X.F 12. 7 1 

906  FORMAT! 10X.11HFIRST  FREO« .F 10.S.7X.1 3HMAX  RESP0NSE-.F12.7/I 

907  FORMAT!  15X.F  10 .5, 21  29X.F  1  2. 7}  1 

906  FORMAT!  10X.11HFIRST  FREO».F 10 .5 .7X.1  3HMAX  R£ SP0NSEa,F|2. 7. 12X. 
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PROGRAM  LISTING  FOR  SAMPLE  PROBLEM  2  (Continued) 


i 


SNERO  -  EFN  SOURCE  STATEMENT  -  IFN(S) 


06/10/6 


117HMAX  RES  RESPONSE*. F12. 7) 

900  FORMAT I ///36X,  2SHPR0GR Am  RAN  TO  COMPLETION) 

910  FORMAT ! 10X.F 10 . S, 31 taX.F12.7l) 

911  FORMATISX.llHFIRST  FREO*.F10.S.9X.9HMAX  RESP* .F 1 2. 7 , 10X.9HM IN  RESP 
1*»F12»7,?X» 1SHMAX  RES  10  RESP*.F12.7> 

912  FORMAT ( IH1.S5X.  10HSHORT  F0RM//4SX . 20H0AMPING  COEFFICIENT* .F10.6/// 
1/9X . 22HFRE0UENC V  (  CYCLES /SEC . ) . 1 9X . 20HM AX  RESPONSE  X  OMEGA, 1 6X . 26HM 
2 AX  RESIO  RESPONSE  X  OMEGA) 

913  FORMAT ( 59H  PLOTTING  ERROR  FOR  X  VS  T.  NUMBER  OF  POINTS  OUT  OF  RAN 
1GE*. 12. 6H  FREQ*. F 10.7// ) 

914  FORMAT ilX.l2.43H  POINTS  OUT  OF  RANGE  FOR  XOOT  VS  T IME .FREQ* »F 10. 7 
1//) 

915  FORMAT!  IHI.SSX.9HL0NG  F0RM//4SX  .20HDAMP I NG  COEFF I C  1ENT* .F  10.8//// 

1 7X* 16HFREQ ( CYCLES /SEC ) • 12X • 20HM AX  RESPONSE  X  OMEGA • 10 X .20HM IN  RESP 
2CNSE  X  OMEGA, BX.22HMAX  RESIO  RESP  X  OMEGA) 

916  FORMAT (4X«  1 1HF1RST  FREO* .F 10.  5 , 10X  , 9HMAX  RESP*  .F 1 2. 7 •  I0X.9HMI N  RES 
1P-.F12.7) 

917  FORMAT! lHl.iax. INI .38X.4HS!  1 )  .34X.7HS2NOI I )// 

1 ! 17X* 14.34X.F12.7. 28X .F12.7) ) 

918  FORMAT! 1H1.S0X. 19H INTERNED  I ATE  VALUES/ 19X,  1H I  ,3SX*9hX0MEGA( I) . 
131X.7HXOOT! I )// 

2 ( 17X. 1 3 . 34X ,F1 2 .7 • 28X.F  12.7 )  ) 

919  FORMAT! IX. 12, 43H  POINTS  OUT  OF  RANGE  FOR  V  VS  TIME ,FRE0*.F10. 7 
1//) 

920  FORMAT! 3F20. a) 

921  FORMAT ! 6F1S.6  > 

922  FORMAT!  IX.  12, 49H  POINTS  OUT  OF  RANGE  FOR  VEL  VS.  FREQ  HITH  ALPHA*. 
1F10.8) 

923  FORMAT!  IX.  12, SOH  POINTS  OUT  OF  RANGE  FOR  RESIDUAL  PLOT  HITH  ALPHA* 

1.  F10.8I 

END 


66 


X  RESPONSE 


Figure  9  -  X  versus  Time  (/PI),  o  =  0,  f  =  98  Cycles  per  Second 
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Figure  10  -  X  versus  Time  (/PI),  a  =  0,  /  =  98  Cycles  per  Second 
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ABSOLUTS  0«  Y  ACCE  LIBATION 


Figure  11  -  versus  Time  (IP2),  a  ■  0,  /  -  98  Cycles  per  Second 
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X  RESPONSE 


VELOCITY  RESPONSE  XDOT 
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Figure  13  -  X  versus  Time  (/Pi),  o  =  0.1,  f  ■  98  Cycles  per  Second 


ABSOLUTE  OK  V  ACCELERATION 


FREQ  =  98.0  CYCLE8/8EC.  TIME  IN  SECONDS 


Figure  14  — <n2X  versus  Time  ( IP2 ),  a  —  0.1,  /-  98  Cycles  per  Second 
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Figure  15  -  Four-Coordinate  Grid  for  Sample  Problem  2,  a  -  0 
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Figure  16  -  Four-Coordinate  Grid  for  Sample  Problem  2,  a  -  0.1 
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Figure  17  —  Four-Coordinate  Grid  for  Sample  Problem  2,  a  -  0 
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Figure  18  —  Four-Coordinate  Grid  for  Sample  Problem  2,  a  -  0.1 
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TABULATION  OF  SHOCK  FUNCTION  FOR  SAMPLE  PROBLEM  2 


SHOCK  SPECTRUM  PROGRAM 
APPLIEO  MATH  LAB  OTmB 


FOUNDATION  VELOCITY 


FOUNDATION  VELOCITY 
-0. 

0*069794998 

-0. 

•0* 

0.167907999 
1. 828629987 
1.628628987 
4.460836994 
4.341246969 
4.341248989 
1.963407967 
1.312149993 
0.670031999 
0.628194993 
0.628194993 
1.926341996 
3.070979983 
3.796847999 
4.131863992 
-0.410769993 
-0.418769993 
-2.289279974 
-0.963170998 
0.321096999 
0.626194993 
1.423017992 
1.670909989 
1.912302990 
1.912302990 
•0.949212000 
-1.679079986 
-2.121767990 
-2.024004974 
“ 1 • 091 311911 
-1.046924993 
-0.432720990 
-0.432720990 
-0.079416999 
-0.339010997 
-0.339019997 
0.809621997 
0.799662999 
0.799662999 
-1.921930906 
-1.921930906 
-1.204227997 
-0.097712994 
-0.097712994 
1.200473994 
2.107008977 
2.047439904 


VS  TIME  I  INPUT! 

TIME  IN  SECONDS 

0. 

0.000124641 

0.000249282 

0.000373922 

0.000496963 

0.000623204 

0.000747849 

0.000672409 

0.000997126 

0.001121774 

0*001246419 

0.001371096 

0.001499697 

0.001620337 

0.001744978 

0.001669619 

0.001994260 

0.002118900 

0.002243049 

0.002360109 

0.002492030 

0.002617471 

0.002742112 

0.002866792 

0.002991393 

0.003116034 

0.003240462 

0.003369323 

0.003409964 

0.003614404 

0.003739249 

0.003063006 

0.003900927 

0.004113160 

0.004237016 

0.004362496 

0.004407097 

0.004611730 

0.004736379 

0.004961020 

0.004909460 

0.009110301 

0.009234942 

0.009399990 

0.009404231 

0.008600072 

0.008733912 

0.008090193 

0.009902794 

0.006107439 

0.004232078 
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I 


tin 


stNom 


1 

0.0607080 

-0.1308000 

8 

-0.0607080 

-0.1308000 

3 

•0# 

0.0607080 

4 

0.1678000 

0.1678000 

S 

l.aautio 

1.4036130 

• 

0. 

-1.6611810 

7 

t. 6811100 

1.6888100 

a 

-0.1308000 

-t. 701 7000 

• 

0. 

0.1308000 

10 

-t. 7770410 

-8.7770410 

11 

-o.tsitoto 

8.8868700 

11 

-0.6411140 

-0.3000880 

13 

-0.0418770 

0.6008370 

14 

0. 

0.0410770 

10 

i.toaia70 

1.8001070 

10 

1. 14443a0 

-0.1038400 

17 

0.7180600 

-0.4107700 

ta 

0.3380160 

-0.3000881 

10 

-4.8806330 

-4.0086400 

to 

-0. 

4.0806330 

tl 

-1.0708060 

-1.0708060 

at 

1.3161080 

3.1066100 

ta 

1.8041100 

-0.0410770 

t4 

0.3070000 

-0.0771300 

ta 

0.7086630 

0.4000600 

ta 

0.4466000 

-0.3400700 

t7 

0.0410770 

-0.4040110 

ta 

0. 

-0.0410770 

to 

-8.0618080 

-8.8618040 

30 

-0.7880600 

8.1307870 

31 

-0.4466000 

0.8701000 

3t 

0.0077130 

0.0444010 

33 

0.4387800 

0. 3300160 

34 

0.8444010 

0.1116780 

sa 

0.6141060 

0.0607000 

30 

-0.6141060 

37 

-0.4466000 

-0.4466000 

30 

0.0016000 

30 

•0# 

-0.8444010 

40 

1. 1446300 

41 

-0.0130800 

-1.1000070 

4t 

0. 

0.0130000 

43 

-8.3171040 

-8.3171040 

44 

•0# 

8.3171040 

48 

0.8373030 

0.8373030 

40 

1.1068180 

0.0408180 

47 

•Po 

-1.1060100 

4a 

1.8001070 

1.8001070 

40 

0.0073380 

-0.3000880 

80 

0.7308870 

-0.1670000 

at 

-0.1308000 

-0.0704170 

at 

0. 

0.1300000 

83 

-1.0100070 

-1.0100070 

84 

-1.8001070 

-0.8701000 

88 

0. 

1.8001070 

86 

-1.0818310 

-1.8818310 

87 

0.9010000 

8.0186100 
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SAMPLE  PROBLEM  2 


SHORT  FORM 
DAMPING  COEFFIC IENT*0. 


FREQUENCY! CYCLES/SEC.)  MAX  RESPONSE  X  OMEGA  MAX 

FIRST  FREQ*  O*  MAX  RESPONSE*  0.097889 


1*00000 

t.00000 

3.00000 

4.00000 

S.00000 

4.00000 

7.00000 

0.00000 

0.00000 

10.00000 

11.00000 

1S.00000 

13.00000 

14.00000 

18.00000 

18.00000 

17.00000 

10.00000 

10.00000 

to.ooooo 

01.00000 

00.00000 

03.00000 

04.00000 

08.00000 

08.00000 

07.00000 

00.00000 

00.00000 

30.00000 

31.00000 

30.00000 

33.00000 

34.00000 

38.00000 

38.00000 

37.00000 

30.00000 

30.00000 

40.00000 

41.00000 

48.00000 

43.00000 

44.00000 

48.00000 

48.00000 

47.00000 

40.00000 

40.00000 

80.00000 

81.00000 


0.10VS09T 
0.1883880 
0.1480487 
0.1081878 
0.8007801 
0.4043840 
0.4088698 
0.8868813 
0.8588787 
0.8013889 
0.8340188 
0.8606716 
0.7348171 
0.8068347 
1.0188488 
1.8341474 
1.4833681 
1.6844868 
1.7734640 
8.84  380  60 
3.0140838 
3.6448861 
4.0876333 
4.8138840 
4.1131108 
3.0480118 
3.6078880 
3.1070430 
8.0843780 
8.8480430 
8.8308108 
1.0808474 
1.8840077 
1.4181010 
1.4080817 
1.8384036 
1.8476387 
1.6838837 
8.0801703 
8.3607807 
8.7181630 
3.0780144 
3.3107400 
3.8401844 
3.6183400 
3.7847884 
3.0378306 
3.0783816 
3.0383370 
3.0843171 
3.7184084 


RES 10  RESPONSE  X  OMEGA 
MAX  RES  RESPONSE* 
0.0978896 
0.1008007 
0.1883880 
0.1480487 
0.1961878 
0.8907801 
0.4043640 
0.49S869S 
0.8868813 
0.4768863 
0.3610668 
0.8814837 
0.3968718 
0.8903004 
0.7187693 
0.6914860 
0.4784660 
0.1808071 
0.6661818 
1.4840801 
8.8430080 
3.0149838 
3.6448861 
4.0876333 
4.8138840 
4.1130483 
3.7084818 
3.3436788 
8.8368744 
8.3836884 
1.0341080 
1.8700869 
1.8879770 
0.0880778 
0.8706067 
0.4035043 
0.4883808 
0.7043330 
1.0018493 
1.3801687 
1.8891307 
8.3879006 
8.0008731 
3.1480448 
3.3800761 
3.3884080 
3.3146841 
3.3007693 
3.3086401 
3.8811847 
3.6979987 
3.7184964 


SO 


3. 919*300 

3.919*899 

gj.aooo* 

3. 1031430 

3.103143* 

94*400## 

a. **at **3 

8.903#7#1 

99.00000 

3.43*4013 

1 . 0060*7* 

90.0000# 

3.1741##* 

1. 11610*1 

9T.#0#M 

1.43443«3 

0*94 384# 2 

9#. #000# 

1.71923*4 

0.833*13# 

94. #000# 

1.9013907 

0.3930*49 

0«.000#0 

1.9*19413 

0. *133413 

•l.#0#«# 

t.*3a*##a 

0.9#1489* 

#a.#o#«# 

l.##7#<#3 

1.4309049 

#9.#00«# 

3.0373M0 

1.0**933* 

64.000## 

3.3*9799# 

8.1*19393 

•o.oooo* 

3.9*#4#30 

3.3974114 

6#.«00#0 

a.#i#4i4* 

8.1943919 

67.000*0 

3.001M17 

1.949*##* 

*#.#00## 

3.1340*#* 

1.7*8174# 

*9.000## 

3.3147##* 

1.9134831 

70.0#0#0 

3.3339937 

8. 39# 3# #4 

71 .0006# 

3.3990993 

8. **47449 

73.0000# 

3.41M##1 

3.413*##1 

73.#00«# 

3. *3393#* 

3.433*309 

74.####« 

3*  9*40073 

3.994M73 

70.#00«# 

3.4379734 

3.389T#78 

7*.#0«M 

3.1*1779# 

8. *84*39* 

7T.#00#0 

t.*34#73a 

3.4M891* 

7#.#0#«# 

a.*#i43#a 

8.1137383 

7t.#0### 

a.3979*#* 

1. #341893 

•«.#0#M 

I. •***•• 7 

1. #13*333 

•l.#«0## 

l.#*#3«94 

1.418333# 

•a.#o««# 

1.9141 ••# 

t.8#l*47# 

#3.#00«# 

1.99347## 

1.874##*# 

•4 .#•••# 

3.9139934 

1.49*984* 

#9.##### 

••#731 t#4 

l.#38##8* 

a#.##### 

S.13l#37# 

8.13783#* 

#7.#00«# 

3.44*3#** 

3.333*4*0 

a#.##### 

B.*39#*a* 

8 .4 19*934 

•o.#oo#o 

3. 731 #3# 3 

3.3*194*9 

*#.#ooa# 

3.91*73## 

8. 39*3*47 

oi.aoo## 

B.4*34##7 

3.1#*7*48 

oa.ao### 

S.994M91 

8. #849339 

*3.000*0 

S.*#49##l 

1. #7*993* 

04.#0##« 

3.7# #• *44 

1 . 74*7993 

oo.aooa# 

3. *34*993 

1 • #3*4*41 

oo.aoo## 

3. #74133# 

1 . 9#**9*3 

07. #00## 

3.7331311 

1.9*7*130 

oa.#oa## 

8.7*#9747 

1.97*3139 

40.0000# 

3. #137073 

1«9*71#8* 

1##.#00#« 

8. #979# 9# 

1.343*914 

l#l.##0## 

8.9## 3*9# 

1.181*131 

i#a.#oo«# 

3. 9411037 

#.*89*30# 

1#3.#Q##0 

8.0##898* 

l.#*179«# 

l#*.##### 

S.#l#3#37 

1.3*1774* 

i#a.#oo## 

3.09497*1 

1. *834*3# 

l#*.#00## 

3.#9#*98* 

8.#43*89t 

1#7.#00M 

3. 1383#19 

8.943441* 

(••.•OO## 

3.9398719 

3. 14010*4 

1#0.#00## 

4.1S**##3 

3.M917M 

U#*#00«# 

4.99M31* 
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111.00000 
11 1.40000 
113.00000 
110.00000 
I 14.00000 
1 10.00000 
117.00000 
110.00000 
110.00000 
100.00000 
Itt. 00000 
IK. 00000 
103.00000 
104.00000 
100*04000 
100.00000 
107.00004 
100.00000 
100.00000 
134.00000 
131.00000 
130.00000 
133.00000 
134.00000 
130.00000 
130.04000 
137.00400 
130.00000 
130.00000 
140.00000 
141.00000 
140.00000 
143.00000 
144.00000 
140.00000 
140.00000 
147.00000 
140.00000 
140.00000 
100.00000 


4.01M3S3 
S.  1009774 
9. 1000000 
9.  3331073 
9.3264410 
9.1743302 
9.0932087 
9.0093000 
9.0400094 
4.9094300 
4.30900  07 
4.4377704 
4.0227702 
4.7900700 
4.0977090 
9.0041432 
9.1990739 
9.3307902 
0.9079200 
0.7734312 
0.1090104 
0.4919042 
0.7010001 
7.0730074 
7.0024002 
0.9371000 
7.3012240 
7.0401749 
7.9700109 
0.3794201 
0.7737249 
9.0009290 
9.1930417 
0.9099740 
0.0090907 
9.0400020 
9.2947002 
9.0440290 
9.9000109 
9.3120000 
SHOOT  POOH 


4.9100393 

9.1009774 

3.1000900 

4.7902090 

4.3203391 

3*9217709 

3.7702740 

3.0011407 

4.0900061 

4.0793949 

3.0303032 

3.3340023 

2.7690202 

2.4491791 

2.0009270 

3.0033700 

3.0292700 

3.9149902 

4.2901470 

4.7779020 

0.4244907 

0.1341902 

0.7340901 

7.0730074 

7.0024002 

0.0039070 

0.3090000 

0.0403009 

0.9201927 

0.0290410 

0.2370101 

0.4090099 

0.7100179 

7.1149190 

7.7130000 

0.4419030 

9.1100070 

9.0310430 

9.0000109 

9.1007001 


DAMPING  COEPPICIENT-O. 09999999 


PREOUENC V ( CYCLE O/SEC . I 
P1N9T  PRIG*  0.  MAX 

1.00000 
2.00000 
3.00000 
4.00000 
0.00000 
0.00000 
7.00000 
0.00000 
9.00000 
10.00000 
11.00000 


MAX  RCOMONS2  X  OMEGA 
RESPONSE*  0.0001190 

0.0379400 

0.0709043 

0.1121004 

0.1470073 

0.1902010 

0.2903012 

0.3192773 

0.3009440 

0.4300131 

0.4000043 

0.0300000 


MAX  RESID  RESPONSE  X  OMEGA 
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12.00000 
13.00000 
14.00000 
IS. 00000 
16.00000 
17.00000 
10.00000 
19.00000 
20.00000 
21.00000 
22.00000 
23.00000 
24.00000 
29.00000 
26.00000 
27.00000 
26.00060 
29.00000 
30.00000 
31.00000 
32.00000 
33.00000 
34.00000 
3S. 00000 
36.00000 
37.00000 
36.00000 
39.00000 
40.00000 
41.00000 
42.00000 
43.00000 
44.00000 
4S.00000 
46.00000 
47.00000 
46.00000 
49.00000 
SO. 00000 
SI. 00060 
62.00000 
83.00000 
84.00000 
SS.00000 
86.00000 
87.00000 
86.00000 
89.00000 
60.00000 
61.00000 
62.00000 
63.00000 
64.00000 
68.00000 
66.00000 
67.00000 
66.00000 
69.00000 
70.00000 


0.8686647 
0.6028072 
0.6711083 
0.8182684 
0.9723428 
1.1099622 
1.2376344 
1.3842498 
1.4817261 
1.8121902 
1.8718036 
1.6611360 
1.7121034 
1.7406730 
1.7244194 
1.7011797 
1.6498966 
1.9699608 
1.4861206 
1.3918268 
1.3046203 
1.2307312 
1.1768394 
1.1412443 
1.1278693 
1.1873616 
1.2236106 
1.3000136 
1.3783184 
1.4817876 
1.8203800 
1.9081121 
1.6169877 
1.6872430 
1.7248740 
1.7714903 
1.7698907 
1.7664930 
1 . 7808688 
1.6021202 
1.0000708 
1.7761844 
1.7419721 
1.7093011 
1.6671361 
1.6672338 
1.6498406 
1.6360376 
1.6334237 
1.6367144 
1.6470296 
1.6606974 
1.6741299 
1.6636090 
1.6947989 
1.7340482 
1.T647997 
1.7641267 
1.7903096 
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00009 

00**0 

7 

000** 

7 

00000 

0000* 

7 

0000* 

7 

00000 

7 

00000 

00000 

a 

0000* 

• 

00000 

000*0 

a 

00000 

a 

00000 

00000 

■ 

00000 

■ 

000*0 

• 

00000 

■ 

00000 

9 

00000 

9 

00000 

9 

0000* 

9 

00000 

9 

00000 

9 

00000 

00000 

9 

•0000 

9 

00000 

« 

00000 

10 

00000 

10 

00000 

10 

00000 

10 

0000* 

10 

00000 

10 

ooooo 

10 

00000 

10 

ooooo 

10 

ooooo 

10 

ooooo 

11 

ooooo 

11 

ooooo 

11 

ooooo 

11 

ooooo 

11 

ooooo 

11 

0000* 

11 

ooooo 

11 

ooooo 

11 

ObCCO 

11 

OOOOO 

12 

OUOOO 

1* 

ooooo 

12 

ooooo 

12 

ooooo 

1* 

ooooo 

1* 

ooooo 

12 

ooooo 

1* 

ocooo 

1* 

•0000 

12 

ooooo 

1.7M2394 
1.7600M2 
1.72949*7 
1 .**94*96 
1. *929079 
l.*ft**297 
I.***4*41 
I. *940*92 
1.7074174 
1.7204*99 
1.73329*1 
1.749*411 
1.79*113* 
1.770120* 
1.781**09 
1.79703*0 
1. *3*49*2 
1. *79227* 
1.9133*3* 
1.99047*4 
1.9*«914* 
2.02142** 
2.0991*** 
2.0**2*23 
2.1203047 
2.1910490 
2.1*04907 
2.2090*42 
2.2392*9* 
2.2«*79*9 
2.2970434 
2.32472*7 
2.3910329 
2.3799299 
2.3993**9 
2.4219209 
2.44334*9 
2.4*32*3* 
*.*01*3*7 
2.49*4304 
2.913**93 
2.92*1791 
2.9411911 
2.9929*43 
*. 9*2*74* 
2.9703*19 
2.9777414 
2.9*3910* 
2.9*797*0 
2.999*001 
2.027BI10 
*.*9*3*72 
*.*••**12 
2.71709*9 
2.744*299 
2.772*449 
2.79*90** 
*.**44790 
*.*409409 


84 


1M.IHH 

t.OTMTOt 

IIWIMH 

t.OOOtf 1) 

130.030*0 

1.01T4004 

I3MHH 

). 0)0410) 

II*. MOM 

I1I.IMM 

). 0040071 

1)0*00  IM 

).001)300 

1)7.0*00* 

3.1100044 

1)0.00000 

3.10M770 

1)0.00000 

3.0104103 

140.00000 

3.0000014 

141.00000 

).31004l) 

140.00000 

3. 3000000 

14). 00000 

3.4000047 

I44.0M00 

3.4431707 

140.00000 

3.4000741 

140.00000 

3.0170100 

147.00000 

3.0000001 

140.00000 

3.0003437 

140.00000 

3.0073000 

100.00000 

3.0334070 
SHOOT  7000 

04*7 INC  COEFFICIENT.*. 

FREQUENCVICrCt.ES/StC.  I 

HAH  RESPONSE  X  OMEGA 

10.00000 

0.M13SM 

0.3*10*00 

00.00000 

0.0430000 

0.0430000 

30.00000 

C.MS04  30 

1. *34 1*0* 

40.00000 

0.3007007 

1.00*1307 

00.00000 

3.0)43171 

3.0*7**07 

00.00000 

I.SOIS4IO 

0.0130413 

70.00000 

3.330*937 

00.00000 

1. 0004007 

1.013*033 

00.00000 

0.01*7300 

100.00000 

0.007*0*0 

1.3430*14 

110.00000 

4.0000310 

4.4310000 

100.00000 

4.  *0040*0 

130.00000 

0.77)4310 

4.777**2* 

140.00000 

0.37*4001 

0.009041) 

100.00000 

*.3100000 

9.1007001 

100.00000 

0.4*03300 

0.01*1*91 

170.00000 

0.10040*0 

100.00000 

7.07*0377 

7.3709419 

100.00000 

10.0070011 

*.007*4*2 

000.00000 

7.3007004 

9.7039270 

010.00000 

7.1770007 

0.1074044 

000*00000 

13.1304070 

I3.I3C407: 

030.00000 

0.07*1000 

3.0037943 

040.00000 

0.7*0373* 

4.9972454 

000.00000 

0.1301740 

4* 1900966 

000.00000 

0.741*034 

♦  .*10441 

070.00000 

0.137*1*0 

4.4807295 

000.00000 

4. 0004704 

1.3037090 

000.00000 

4.7)00)4* 

1.0047000 

300.00000 

4.  7*40*** 

3. 5465514 

310.00000 

4.M4070S 

2.2577650 
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4.8112682 
5.6847954 
9.9516458 
7.664261* 
8.3929031 
6.7073254 
5.3967382 
5.4726837 
7.9915000 
1C.0929618 
5.6722559 
4.9223868 
6.2064720 
16.6626980 
12.5990441 
14.0287548 
8.5671775 
7.6853912 
7.0414768 
6.8870090 
9.0612003 
9.5879017 
12.4685265 
11.7912351 
13.3713596 
15.4363683 
14.8852375 
14.0197783 
14.2594140 
15.1948453 
18.7635851 
12.7766875 
9.0435599 
10.9088503 
14.4142867 
I6.608C47C 
39.9168267 
41.0169787 
37.2481532 
33.5364890 
33.1622615 
33.0957823 
26.6072156 
36.0618491 
26.2506065 
24.2740149 
39.0256066 
21.6589115 
29.5010610 
57.1597290 
37.5797143 
40.2439532 
20.9892118 
21.1313868 
21.1676102 
21.1024454 
(C. 9411793 
*5.1282744 
*9.0582066 


3.6082111 
2.2581167 
9.5493848 
6.1251028 
7.8221811 
5.3521238 
0.9335095 
4.5590934 
5.5663314 
8.8656384 
2.4946889 
4.0816076 
5.8942370 
15.3588334 
8.5707424 
14.0240724 
5.8224086 
3.2154393 
4.0217979 
6.3562272 
6.5626441 
1.2931048 
11.7681925 
11.2566628 
11.7661412 
8.2025073 
11.5208479 
6.3820055 
8.930 64j9 
14.0344567 
18.5134978 
8.6071843 
5.8823307 
2.9627771 
13.1337078 
7.3286763 
39.4246798 
37.3638587 
27.6386724 
26.4047632 
27.3717346 
31.7338808 
6.9213288 
28.1872778 
11. 3558303 
1.6996583 
37.8539357 
9.0666267 
23. 7623100 
57.1374674 
28.64607S3 
37.9155636 
8.5*53859 
7.1144326 
13.960  575 
16.2059939 
13.9906071 
23.0315237 
25.530371C 
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910.00000 

31.6416600 

920.00000 

37.9951039 

930.00000 

30.0013939 

940.00000 

31.5943122 

950.00000 

31.3997571 

960.00000 

29.7507390 

970.00000 

33.4940093 

900.00000 

27.3976500 

990.00000 

26.5197105 

1000.00000 

23.5731676 

SHORT  FORM 

DAMPING  COEFFICIENT*!:. 09999999 

FREQUENCY (CYCLES /SEC. ) 

MAX  RESPONSE  X  OMEGA 

10.00000 

0.4065043 

20.00000 

1.4517261 

30.00000 

1.4001206 

40.00000 

1.3703154 

50.00000 

1.7005650 

60.00000 

1.6334237 

70.00000 

1.7903096 

00.00000 

1.7204095 

90.00000 

1.9504764 

100.00000 

2.2607969 

110.00000 

2.4904304 

120.00000 

2*5956001 

130.00000 

2.0722792 

140.06000 

3.2005014 

180.00000 

3.6334579 

160.00000 

3.7313607 

I70.0J000 

3.0233107 

100.00000 

4.1514060 

190.00000 

4.3373303 

200.00000 

4.3296462 

210.00000 

4.1441401 

220.00000 

3.0000703 

230.00000 

3.3555207 

240.00000 

3.0701179 

280.00000 

3.2349019 

200.00000 

3.3023309 

270.00000 

3.4044560 

200.00000 

3.5450266 

290.00000 

3.0094345 

300.00000 

3.0499451 

310.00000 

3.0600000 

320.00000 

3.0544473 

330.00000 

3.0123303 

340.00000 

3.5450590 

380.00000 

3.4529017 

360.00000 

3. 3350006 

370.00000 

3.2902506 

300.00000 

3.3141514 

390.00000 

3.3343120 

400.00000 

3.3850110 

410.00000 

3.3740913 

420.00000 

3.3914393 

27.061 1905 
35. 0295106 
37.6041909 
31.0100202 
24.0365034 
21.4417690 
31.9906242 
10.0603947 
24.4730940 
23.2734306 


MAX  RES  ID  RESPONSE  X  ORES* 
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430 

.00000 

440 

.00000 

450 

.00000 

460 

.00000 

470 

.00000 

460 

.ooooo 

*90 

.00000 

500 

.00000 

510 

.06000 

520 

.00000 

530 

.00000 

540 

.06000 

550 

.00000 

560 

•ooooo 

570 

.00000 

580 

.00000 

590 

.06000 

600 

.06000 

610 

.00000 

620 

.09000 

636 

.00000 

640 

*00000 

650 

.06000 

660 

.09000 

670 

.00600 

660 

.09000 

690 

.00006 

700 

.00000 

710 

.00000 

720 

.00600 

730 

.06000 

740 

.ooooo 

750 

.00000 

760 

.00000 

776 

.06000 

780 

.09000 

790 

.00000 

000 

.00000 

610 

•ooooo 

020 

.00000 

630 

.00000 

640 

•ooooo 

650 

.00000 

660 

.ooooo 

070 

•ooooo 

600 

.ooooo 

690 

.ooooo 

900 

.00000 

910 

.00000 

920 

.00000 

930 

.00000 

940 

•ooooo 

950 

.00000 

960 

.00000 

970 

.00000 

*60 

.00000 

990 

1666 

.00800 

3. *0*6858 
3.414663* 
3.4214000 
3.4240610 
3.5307016 
3.6127075 
3.7326111 
3.0307307 
3.0460405 
4.0432333 
4.1256106 
4.2130045 
4.2060452 
4.3506643 
4.4250700 
4.7210372 
4.0707236 
5.2252000 
5.4502060 
5.6026020 
5.0103302 
6.1351700 
6.3670600 
6.5542150 
6.7350400 
6.0506504 
6.0155040 
6.0157705 
6.0361212 
6.0703006 
6.0702772 
7.0420007 
7.3435100 
7.6167054 
7.0047350 
0.1104016 
0.3164675 
0.4006007 
0.6310056 
0. 7620024 
0.0406022 
0.0005405 
0.0416541 
0.0300007 
0.0130564 
0.0600607 
0.7760501 
0.6021544 
0.5376770 
0.4046476 
0.2147046 
0.1023047 
0.0  1  50021 
7.0002602 
7.7010631 
7.7011210 
7.6215640 
7.5117304 


320.00000 

330.00000 

340.0COOO 

350.00000 

360.00000 

370.00000 

380.00000 

300.00000 

400.00000 

410.00000 

420.0C000 

*30.09000 

440.00000 

450.00000 

460.09000 

470.00000 

460.00000 

490.00000 

500.00000 

510.00000 

520.00000 

530.00000 

540.00000 

sso.ocooo 

560.00000 

570.00000 

560.00000 

590.00000 

600.00000 

610.00000 

620.00000 

630.00000 

6*0.00000 

650.00000 

660.00000 

670.00000 

660.00000 

69O.OCC0O 

700.00000 

710.00060 

720.09000 

730.00600 

740.00000 

750.00000 

760.00000 

770.00000 

780.00000 

790.00000 

600.00000 

610.00000 

620.00000 

630.00000 

640.00000 

650.00000 

660.00000 

670.06000 

660.00000 

690.00000 

900.00000 
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SAMPLE  PROBLEM  2 


LCNG  FORM 

DAMP  I  NO  COEFFICIENT*!). 


FREOICVCLES/SEC ) 
FIRST  FREO*  0. 

1.00000 

2.00000 

3.00000 

4.00000 

S.  00000 
0.00000 

T.  00000 
0.00000 
9.00000 

10.00000 
11.00000 
12.00000 
13.00000 
14.00000 
IS. 00000 
16.00000 
17.00000 
10.00000 
19.00000 
20.0000C 
21.0000'' 
22.00000 
23.00000 
24.00000 
29.00000 
26.00000 
2T.C0000 
26.00000 
29.00000 
30.00000 
31.00000 
32.00000 
33.00000 
34.00000 
39.00000 
36.00000 
3T. 00000 
30.00000 
39.00000 
40.00000 
41.00000 
42.00000 
43.00000 
44.00000 
49.00000 
46.00000 
47.00000 
40.00000 
49.00000 
90.00000 
91.00000 


MAX  RESPONSE  X  OMEGA 
MAX  RESP*  0.0040391 
0.0260267 
0.CS6309I 
0.0996439 
0. 1427370 
0. 1964645 
0.2666994 
0.3427713 
0.4227999 
0.9014  796 
0.9723199 
0.6340126 
0.6896716 
0.7342171 
0.7821219 
C. 8111041 
0.8466261 
1.0913280 
1.3739737 
I. 7041911 
2.2223926 
2.7704621 
3.2003321 
3.9199902 
3.722663*1 
4.1131192 
3.9420119 
3.6970990 
3.1070430 
2.0243790 
2.41 441 1  0 
2.0047071 
1.6124741 
1.3241403 
1.2929934 
1.300692  2 
1.3039326 
1*2996910 
1.3961394 
1.6670373 
2.12724  90 
2.9020031 
2.7071109 
2.9391120 
3.2419096 
3.9292770 
3.7247994 
3.0372306 
3.0723216 
3.6039960 
3.0243171 
3.6442407 


PIN  RESPONSE  X  OMEGA 
MIN  RESP*  -0.0061196 
-0.0384237 
-0.0769012 
-0.1144403 
—0. 1519441 
-0.222230C 
-0.2530273 
-3.3190194 
-3.44967CC 
-0.9599767 
-0.6013559 
-0.9619901 
-C. 4292086 
-0.5924194 
-0.0060347 
-1.0185422 
-1.2341474 
-1.4233621 
-1.6244966 
-1.7734649 
-2.1979271 
-2.6311766 
-3.3646207 
-3.0109224 
-4.0739144 
-4.0671614 
-3.0646404 
-3.9441637 
-3.1640061 
-2.7722511 
-2.0494430 
-2.2300102 

-1.0006474 

-1.6940077 

-1.4191010 

-1.4924217 

-1.9394036 

-1.9476327 

-1.6932937 

-2.0291703 

-2.3607207 

-2.7121639 

-3.0720144 

-3.3107409 

-3.9401544 

-3.6193490 

•3.6449209 

-3.7997294 

-3.0210937 

-3.0393370 

-3.7214329 

-3.6339170 


MAX  RESIO  RESP  X  OMEGA 
MAX  RESIO  RESP*  3. :976896 
0.1045997 
0.1263260 
0.1480467 
0.1961672 
0.2907201 
0.4043640 
0.4992699 
0.9262613 
0.4760863 
0.3610662 
0.2614637 
0.3969710 
0.5903394 
0.7157693 
0.6914869 
0.4754660 
0.1609071 
0.6661912 
1.4249041 
2.2430060 
3.0149630 
3.6440261 
4.0970333 
4.2130940 
5.1130403 
3.7904012 
3.3436722 
2.0369744 


1.9341904 

1.5700969 

1.2279779 

0.0000779 

0.9700007 

0.4039043 

0.4023209 

0.7043330 

1.0019493 

1.3091627 

1.0991307 

2.3979096 

2.0092731 

3.1424442 

3.3200761 

3.3994050 

3.3140041 

3.3007693 

3.3096401 

3.5911247 

3.6979927 

3.7124464 
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so. ooooo 

3.3427003 

-3.2684262 

3.5196200 

S3. 06000 

2.0765460 

-2.0322416 

3.1031439 

54.00000 

2.6021543 

-2.3074032 

2.5030701 

SS .00000 

2.4364013 

-2.0400121 

1.0040076 

SO. 00000 

2.1741006 

-1.7467531 

1.1161041 

sr.ooooo 

1.0244303 

-1.4666013 

0.8432402 

SO. 00000 

1.7102324 

-1.4474320 

0.2330130 

so.ooooo 

1. SO  13007 

-1.4004230 

0.3530640 

so. 00000 

1.5615412 

-1.8006125 

0.6132413 

01.00000 

1.6326002 

-1.5303022 

0.0014209 

00.00000 

1.0070203 

-1.880 0450 

1.4300040 

03.00000 

2.C372060 

-1.6746003 

1.0648234 

04. 00000 

2.2057000 

-1.0176130 

2.1018383 

os.ooooo 

2.5604030 

-2.1100600 

2.2074114 

00.00000 

2.0104144 

-2.3320462 

2.1042018 

67.00000 

3.6010617 

-2.4018858 

1.0480000 

00.00000 

3.1240505 

-2.6166604 

1.7621740 

00.00000 

3.21470S4 

-2.7408800 

1.0134231 

TO. 00000 

3.3320537 

-2.0008788 

2.3063404 

71.00000 

3.3500552 

-3.1033287 

2.0647440 

72.00000 

3.2030400 

-3.1726068 

3.4134001 

73.00000 

3.3SSS246 

-3.4848123 

3.6338300 

74.00000 

3.5041400 

-3.8003002 

3.5040073 

TS.00000 

3.437S724 

-3.4130071 

3.3207072 

70.00000 

3. 101 7750 

-3.1471448 

8.0240300 

77.00000 

2.0262S43 

-2.0240732 

2.4002010 

70.00000 

2.8341170 

-2.0014302 

2.1137323 

70.00000 

2.2240130 

-2.2070606 

1.0341203 

00.00000 

1.0407114 

•1.0066047 

1.6130233 

01.00000 

1.0307004 

-1.0003684 

1.4123330 

oa. ooooo 

1.0010481 

-1.0141068 

1.2010470 

03.00000 

1 .0824744 

-1.0207420 

1.2740000 

04.00000 

2.0128024 

-1.0480382 

1.4000240 

os.ooooo 

2.0721104 

-1.0000604 

1.0226020 

00.00000 

2.1310270 

-2.1123307 

8.1872308 

07.00006 

2.4402006 

-2.3000000 

8.3380400 

00.00000 

2.6200020 

-2.0002030 

2.4100024 

00.00000 

2.7310303 

•2.0017700 

2.3010400 

00.00000 

2.8008307 

-2.8147308 

2.2003667 

01.00006 

2.4024007 

-2.4001867 

2.1667062 

oa.ooooo 

2.8010772 

•2.8040801 

8.0260380 

03.00000 

2.0247053 

-2.6048001 

1.0700030 

04.00000 

2.6686100 

-2.7000044 

1.7607002 

08.00000 

2.6246803 

-2.8008787 

1.6366661 

00.00000 

2.4741230 

•2.4000838 

1.0000063 

07.00000 

2.7221211 

-8.3012780 

1.0070130 

00.00000 

2.7608747 

-2.1424810 

1.0703130 

00.00000 

2.8137072 

-8. 1412381 

1.0071080 

100.00000 

2.0070008 

-8.1820202 

1.3630016 

101.00000 

2.0003000 

-2.1028888 

1.1210131 

tea. ooooo 

2.0411037 

-2.1726008 

0.0000200 

103.00000 

2.0002024 

-2.1026203 

1.0017000 

104. OOOOO 

3.0103037 

-2.1026143 

1.2617760 

IOS. ooooo 

3.6848701 

-2.80I3I8T 

1.6836030 

100. ooooo 

3.0000820 

-8.0470473 

2.0430206 

107.00000 

3.1223018 

-3.1147000 

2.8484010 

100.00000 

3.8302710 

•3.8300304 

3.1401006 

100.00000 

4.1200603 

-4.0078887 

3.0001700 

110.00000 

4.0804310 

-4.8438200 

4.4312800 
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111.00000 
112.00000 
113.00000 
1 1 • • 00000 
1 19*00000 
1 10*00000 
117.00000 
110.00000 
119.00000 
120.00000 
121.00000 
122.00000 
123.00000 
124.00000 
120. 00000 
120.00000 
127.00000 
120.00000 
124.00000 
130.00000 
131.00000 
132.00000 
133.00000 
134.00000 
130.00000 
130.00000 
137.00000 
130.00000 
139.00000 
140.00000 
141.00000 
142.00000 
143.00000 
144.00000 
140.00000 
140.00000 
147.00000 
140.00000 
149.00000 
100.00000 


4.0124017 
S.CI43716 
9.C349724 
4.0429094 
9.0972102 
9.1743202 
9.0992927 
9.0093004 
9.0440994 
4  .  40  9  4  990 
4.9690097 
4.4016733 
4.6124903 
4.7602334 
4.0772433 
9.0020366 
9.1003609 
9.3907962 
9.9079209 
9.7734312 
0.1650104 
6.3412770 
0.7054290 
6.0442219 
7. 0006906 
0.4371600 
7.1246003 
7.5329942 
7.4470620 
0.379*201 
0.7737244 
9.C66S246 
4.1936417 
0.9059740 
0.6090907 
9.0400020 
9.2123531 
9.5446290 
4.5344342 
9.2207503 


-4.0370142 
-4.9063667 
-4.9931400 
-5.2221972 
-5.2264419 
-5.059004C 
-4.0921160 
-4.7221107 
-4.4204204 
-4.0070023 
-4.2516337 
-4.4377704 
-4.6227762 
-4.7990706 
-4.9577050 
-5.0041432 
-5.1950735 
-5.2777900 
-5.3340325 
-5. 6110770 
-6.1403020 
•6.4919042 
-0.7610561 
-6.9613044 
-6.9003497 
-6.9201233 
-7.3012240 
-7.6401749 
-7.9750109 
-0.2769*24 
-0.3973036 
-0.3927690 
-0.2740026 
-0.2579001 
-0.6720413 
-9.002*075 
-9.2947002 
•9.41001*3 
-9.5100416 
-9.3120*0* 


4.91*0303 
5.1605774 
5.10**900 
4.7902650 
4.3263391 
3.92177*0 
3.7702740 
3.0011407 
4.05*00*1 
4.0795549 
3.0303032 
3.3340623 
2.7090262 
2.44*1791 
2.6009270 
3.0533760 
3.5252700 
3.91*9*02 
4.2901470 
4.7  779020 
0.4**4907 
6.13*1*02 
0.73*0*01 
7.0736*7* 
7.002*0*2 
0.003*070 
0.3**005* 
0.0*03*0* 
0.9201*27 
0.02*541* 
0.23701*1 
0.4090S9* 
0.7100179 
7.114*100 
7.713*00* 
0.4*19020 
9.1150070 
9.5310*30 
9.59**10* 
9.10*70*1 
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SAMPLE  PROBLEM  2 


LCNG  FORM 

DAMPING  COEFFICIENTS. ;9«99999 


FREU(CVCLCS/SFC  I 
FIRST  FWEg  =  (.  . 

I  •  *  J03- 
it.swoc; 

3*  CUGOG 

A. roOPO 
S.  jICOi, 

6 • CCCC 

r.coooi. 

B. CCOjj 
v.caoit 

I0.CC0CC 

ii.oc.cc>-. 

I2.OS0GG 
13.C30G*. 
14.00000 
15. COO  CC 
16.00000 
ir.coooc 
16.00000 
19.00000 
20.00000 
2I.00CCC 
22.00000 
23.00000 
24.00000 
25.00000 
26.00000 
27.00000 
20.00000 
29.00000 
30.00000 
3I.0000C 
32.00000 
33.00000 
34.00000 
30.00000 
36.00000 
37.C0000 
30.00000 
39.00000 
40.00000 
41.00000 
42.00000 
43.00000 
44.00000 
40.00000 
40.00000 
47.00000 
40.00000 
49.00000 
50.00000 
51.00000 


MAX  MFSPONSE  X  OMEGA 
MAX  MtSMs  n. 0040351 
t.'  271717 
C.C577969 
C.C  959644 
C.  1422809 
C. 1962510 
C. 2563312 
U. 3192773 
0.380944C 
b. 4366131 
6.4565643 
0.5306656 
0.5656647 
0.6C25072 
0.6251661 
0.64  80966 
0.6692952 
0.7266459 
0.8622626 
1.04c  I  653 
1.2497623 
1.4257137 
I. 5715036 
1.6611300 
1.79  66653 
1.7240752 
1.7143936 
1.7011797 
1.6455966 
1.5659605 
1.4681206 
1.3915266 
1.2916945 
|. 20 35173 
1.1395501 
1.1086662 
1.1140342 
1.1530259 
1.2179627 
1.2976037 
1.3783154 
I. 4517576 
1.5241575 
1.5651121 
1.6165577 
I.6S7243C 
1.7248740 
1.7714903 
1.7690907 
1.7864930 
1.7605658 
1.6021202 


KIN  RESPONSE  X  OMEGA 
KIN  RESP«  -0.0061196 
-0.0379406 
-0.07S9543 
-C. 1121664 
-0.1470673 
-0.1607973 
-0.2130117 
-C. 2600  0  7  3 
-U. 3695401 
-0.4343356 
-0.4529236 
-0.4216341 
-C.  3830994 
-0.5247611 
-0.6711083 
-0.8182864 
-0.9723425 
-1.1059622 
-1.2376344 
-1.3542456 
-1.4517261 
-1.5121502 
-1.5467764 
-1.6475601 
-1.7121034 
-1.7400730 
-1.7244194 
-1.6394011 
-1.5326239 
-1.4993467 
—  I  .450  4  690 
-1.3873041 
-1.3046203 
-1.2307312 
-1.1766354 
-I. 1412443 
-1.1275653 
-I. 1573016 
-1.2230106 
-1.3000130 
-1.3720370 
-1.4453105 
-1.5203500 
-1.5747020 
-1.5797869 
-1.5463766 
-1.5363979 
-1.5747972 
-1.5607947 
-1.5943269 
-I. 5952716 
-1.5794198 


MAX  RES10  RESP  X  OMEGA 
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53.00000 

1.6000708 

-1.0444130 

83.00*00 

1.7701844 

-1.8004833 

84.00000 

1.7419731 

-1.4611470 

88.00000 

1.7083011 

-1.4138031 

80.00000 

1.6071301 

-1.3673101 

87.00000 

1.6073338 

-1.3886043 

80.00000 

1.6403486 

-1.374(170 

80.00000 

1.6330376 

-1.3928117 

00.00000 

1.6334337 

-1.4108703 

01.00000 

1.6367144 

-1.4364043 

03. 00000 

1.6470306 

-1.4440089 

03.00000 

1.6606874 

-1.4633819 

04.00000 

1.6741300 

-1.4004067 

08.00000 

1.4030090 

-1.4973849 

00.00000 

1.6947889 

-1.8139981 

07.00000 

1.7340483 

-1.8303004 

08.00000 

1.7847897 

—1.8408387 

00.00000 

1.7341367 

-1.8414309 

70.00000 

1.7903096 

-1.8700787 

71.00000 

1.7033394 

-1.8934000 

73.00000 

1.7600003 

-1.0000187 

73.00000 

1.7384907 

-1.0334901 

74.00000 

1.8004006 

-1.0301310 

78.00000 

1.6379006 

-1.0838078 

76.00000 

1.8713713 

-1.0408307 

77.00000 

1.8139166 

-1.0004041 

70.00000 

1.4893313 

-1.4940003 

79.00000 

1.4991477 

-1.7074174 

00.00000 

1.8434306 

-1. 7104090 

01.00000 

1.6074370 

-1.7331901 

03 .00006 

1.4300033 

-1.7480411 

03.00000 

1.4738007 

-1.7801136 

04.00000 

1.7184070 

-1.7701100 

es.ooooo 

1.7864038 

-1.7018008 

06.00003 

1.797038* 

-1.7933164 

07.00000 

1.8344983 

-1.0048338 

00.00000 

1.3783370 

-*.0184808 

09.0000U 

1.9133030 

-1.0101030 

9e.ooooo 

1.9804704 

-1.0304030 

9I.O00OC 

1.4068140 

-1.0460730 

92.CO000 

1.0314300 

-1 .0000398 

93.00000 

3.0881000 

-1.0607(30 

94.00003 

3.0003033 

-1.0703441 

98.00000 

3.130  30  47 

-1.0080978 

46  .1/0000 

1.1810480 

-1.0947040 

97.00000 

1. 1004807 

-1.9038000 

98.00000 

3.3090043 

-1.9131373 

99.C0003 

1.3391080 

-1.9104101 

100.00000 

3.3607909 

-1.9204000 

101 .00060 

3.3970434 

-1.9361360 

103. 1000 o 

1.33  473  0  7 

-1.9438074 

163.00000 

1.3810339 

-1.9807690 

104.30000 

3.3789380 

-1.9870747 

108.00000 

3.3943609 

-1.9443004 

106. COCOg 

3.4319309 

-1.9700000 

107.90000 

(.4433409 

-1.9707830 

I06.0000C 

3.463(630 

-1.9038033 

109.00090 

2.4810337 

-1.9000970 

110. 00090 

2.4904304 

-1.9933078 
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111.90000 

2.51 36293 

-1.9983436 

1 1 2.00000 

2. 52  61 75 1 

-2.0030557 

113.00000 

2.5411911 

-2.0074919 

114.00000 

2.5525643 

-2.0 1 1 6535 

US. 00000 

2.5622746 

-2.0155394 

116.00000 

2.5703015 

-2.0191509 

117.00000 

2.5777414 

-2.0224874 

no. ooooo 

2.5835108 

-2.0255508 

119.00000 

2.5875728 

-2.0283419 

120.00000 

2.5956801 

-2.0306568 

121.00000 

2.6275118 

-2.0549282 

122.00000 

2.6583672 

-2.1308325 

123.00000 

2.6882212 

-2.2056735 

124.00000 

2.71 70503 

-2.2793104 

125.00000 

2.7448299 

-2.3523175 

126.00000 

2.7722445 

-2.4243357 

1 27 .00000 

2.7989026 

-2.4946165 

128.00000 

2.8244750 

-2.5636270 

129.O0C0C 

2.8489405 

-2.6315377 

130.00000 

2.8722792 

-2.6980584 

131.00000 

2.8952513 

-2.7626036 

I32.C000C 

2.91  74094 

-2.8253742 

133.00000 

2.9384122 

-2.8871990 

134.00000 

2.9582444 

-2.9467676 

135.00000 

2.9776415 

-3.0C48671 

136.00000 

2.9962702 

-3.0613390 

137.00000 

3.0147992 

-3.1 156944 

138.00000 

3.0336479 

-3.1688772 

139.00000 

3.0523041 

-3.2194193 

140.00000 

3.07C7277 

-3.2685814 

141.00000 

3.0882752 

-3.3155413 

142.00000 

3. 1057992 

-3.3596568 

143.00000 

3. 1222547 

-3.4028847 

144.00000 

3.1375520 

-3.4431707 

145.00000 

3.1522428 

-3.4806741 

I46.COOOO 

3. 1664564 

-3.5170129 

147.00000 

3.1795024 

-3.5502561 

148.00000 

3. 1913713 

-3.5803437 

149.00000 

3.2021606 

-3.6075555 

150.00000 

3.21  29558 

-3.6334579 
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SAMPLE  PROBLEM  2 


CONG  FORM 

OAMPING  CQEFF I C I ENT*0 • 


FREGtOCLES/SEC  »  fax  RESPONSE  x  QFEGA  PIN  RESPONSE  X  QMEGA  MAX  RESlO  RESP  X  OMEGA 


IO.OOOOC 

0.5723199 

-0*6013959 

0*3610662 

20.  FQOOG- 

2.2223526 

-2*1979271 

2.2430000 

30. coco: 

2.4144116 

-2*9499436 

1.9341909 

40.SCOOO 

2. 1272496 

-2*3607207 

1*0591367 

5C  .COOOO 

3.6243171 

-3*7214325 

3*6979927 

6C.000CV, 

1.5619412 

-1*9396125 

0.6132413 

73.COOO: 

3.3329537 

-2*6699755 

2*3963664 

flC.COOCw 

1.6407114 

-1*6966067 

1*6139233 

90.CCC0C 

2.5065307 

-2*9147305 

2.2*636.7 

100**000:. 

2.8579095 

-2*1520202 

1*3436914 

UO.OOOCt 

4.5506316 

-4.9435206 

*..31230. 

I20.COOOC 

4.9054596 

-4*0670623 

4.C795549 

uo.cooc: 

5.7734312 

-5*6116770 

4.7779626 

1 40.0000 J 

6.3794201 

-6*2769624 

6.S295410 

ISO.COOCC 

9.2287503 

-9.31266C6 

9*1567661 

I60.C0CCC 

6.2057672 

-•..*0  3336 

9.C191691 

170. JOOOC 

6.0960002 

-6*1024095 

6*663U  763 

leo.ccocc 

7.6701775 

-7. .7*2377 

7.37.5*1* 

110.3000 

9.9505621 

-t  C  *097661 1 

9*6079402 

2cc.:cooc 

7.3057204 

-6*9496645 

9*7939270 

210*0003 C 

7. 1*5728® 

-7.1775007 

6.  127.6*. 

220.: 300 : 

12. .272.92 

-1 2. *5276. C 

13.1304.70 

230.0OCG0 

6.0791562 

-6*7040646 

3.6237543 

2AC. 00000 

6.7953739 

-6.63«6*97 

..*•726*4 

250.CO00C 

4.9294916 

-3.13517*6 

4*1900060 

2«3.3tiC8C 

9.7415234 

-3.726.6.C 

... 100*01 

270. COOOO 

9.0764190 

-3.13761*6 

4*4607299 

28C.OOOOG 

». 6«0*7g« 

-3. ..97976 

1.3.3769. 

29C.C00C0 

4.7303349 

-3*2965604 

1.62470.2 

30C  *00003 

4.7940996 

3*6469514 

310.00000 

4.9940769 

-4*6510790 

2*2977690 

320 .03000 

4.01 12662 

•3*6694795 

3. 60 .2111 

330.03000 

-9*0 J39369 

2.29.1167 

340.C3QCC 

9.9463107 

-9*9916453 

9*5493040 

3SC.33000 

7.2914316 

-7.66.26I2 

6*  1 291C20 

360 .03003 

8.3929031 

-0*3669973 

7*6221611 

370 . *0009 

6.7073294 

-•.6.796.7 

9.392123. 

36C.0000C 

9.3967302 

-4*7079370 

0.9339C99 

390.00000 

9.4726637 

-9.00691.2 

..95*3*3* 

400.00000 

7.3466106 

-7.991 900C 

9*5663314 

4 10 .30000 

9.0315334 

-W  *0929615 

...69630. 

420.00000 

9*6722999 

-9*6121054 

2*4946660 

430.00000 

4.9223666 

•4*9163995 

4*0016076 

440.00000 

6.2664720 

-6*0059653 

9*6942370 

490.00000 

-16*4112132 

19*3960334 

400.00000 

12*0612262 

-12*5990441 

0*5707424 

470.00000 

14.0207940 

•13*0333476 

14*0240724 

400.00000 

6*3241466 

-..9.71779 

9*0224? 06 

490.00000 

7*6093912 

-•.*•01313 

3.219*3*3 

SCO • 000CO 

6*9662094 

-7*0414766 

4*0217970 

910.0000C 

6*6649709 

-6*0670090 

6.39*2272 

920.00000 

9*0612003 

-9*0344149 

95 


530 

•  coooo 

540 

•  00000 

550 

•  OOOCO 

560 

.<*000** 

570 

•  coooo 

560 

.00000 

590 

•  coocc 

60C 

•  30CCC* 

610 

.eooco 

620 

•  coooo 

630 

•  OOOCO 

64C 

•  :;o  cv 

650 

•  OOOCO 

660 

.00000 

670 

.:ooor 

680 

.0000' 

690 

.00000 

700 

.00000 

710 

•  cooco 

721* 

.1000C 

7  3ft 

•  ftOOCO 

740 

•  OOOCO 

750 

.coooo 

760 

770 

•cococ 

780 

.cftoc: 

790 

.JftftOO 

6wC 

810 

.ooooc 

620 

•  COft'JJ 

630 

.30000 

640 

•  09009 

650 

•  OOOCft 

860 

•ooocc 

670 

.  ft  C  C-  v  ■» 

060 

.CCOOw 

690 

•OOOCo 

909 

•ccoco 

910 

.0350- 

920 

.  C  «*  J  P  * 

930 

•  003CC 

940 

•coooo 

950 

•  ?CuCl 

960 

•330CO 

970 

•ftoocc 

960 

•coo:c 

99C 

1000 

.30  0  00 

9.5879017 
12.2454479 
11.7912351 
13.3713596 
15.4383653 
14.6800256 
13. 9738264 
13.5483565 
15. 1948453 
18.7635851 
12. 7766675 
9. C  4  35599 
1C. 9068503 
14.4142667 
15.4063463 
39.9166267 
40.9539375 
37. 1573563 
33.5364690 
32.9295487 
33.C957623 
25.9561469 
3S.829C7C6 
25.2506065 

24.  C4  13256 
3  9.02 56  J 66 
21.6389115 
29.35C3990 
57. 1597290 
37.5707143 
39.77925C6 
19.8272700 
19. 7017763 
20.4062905 
29.S70318C 
2C.S93H89 

25 .  C 1 97773 
28.5140569 
31.3670921 
37.6792321 
38 .8544083 
31.5943122 
31.3997571 
29. 75  C  7398 
33.2609253 
27.3976500 
26.5197195 
23.5656412 


-9.4926915 
-12.4665265 
-11.7773446 
-13.1167963 
-1 5. 3C76662 
-14.6652375 
-14.0197783 
-14.2S9414C 
M 4.9064192 
-18.4727407 
-12.6451106 
-8.9427893 
-10.6612426 
-14.3406677 
-15.6C8047I. 
-39.6359294 
-41.0169787 
-37.2461532 
-33.5256643 
-33.1622615 
•32.6447542 
-26.6072156 
-36 • C  61 849 i 
-25.0946736 
-24.2740149 
-39.0116711 
-21.4430659 
-29.501061 C 
-56.9387345 
-37.2578034 
-40.2439532 
-29.9692118 
-21.1313868 
-21.I6761C2 
-21.1.24454 
-20.9411763 
-25.1282744 
•29.0562066 
-31.6416666 
-37.9951839 
-36.6013939 
-31.3670361 
-31.270514C 
•26.8960285 
-33.4948893 
-27.0516916 
-26. 1C  39162 
-23.5731676 


1 . 2931C48 
11.7681925 
11.2566626 
11.7661412 

8. 202SC73 
11.5206*79 
6. 362 v ' 55 
6.93C6409 
14. C  344567 
16.5134978 
6.6671643 
5.6823307 
2.9627771 
1  3. 1 337r78 
7.3266763 
39.4246790 
37.3636507 
27.6366724 
26. 4v 47632 
27.3717346 
31.7336608 
6.9213266 
29.1672776 
11.3556303 
1.6996563 
37.6539357 
9.  666267 
23. 76231C0 
57.1374674 
26.646:753 
37.9155636 
6.5253059 
7.1144326 
13.96C1S75 
16.2359939 
1 3. 9906071 

23.  *315237 
25. 53C  3710 
27.' 61 1965 
35.8295166 
37.6641969 
31 166262 

24. (365634 
21.4417696 
31.9906242 
10.8603947 
24.4736946 
23.2734306 


SAMPLE  PROBLEM  2 


LCNG  FORM 

DAMPING  C0EFF|CI€NT»0. 09999999 


FREC ( CYCLES  / SEC  I  max  RESPONSE  X  OMEGA  MIN  RESPONSE  X  OMEGA  MAX  RES  ID  RESP  X  OMEGA 


IO.OOOJO 

a. 4865843 

•0*4529236 

20 • C  UU  OC 

1. 2497623 

-1*4517261 

30 • 0000 v 

1  • 48  81 206 

-1.456 4690 

40 • JOOOc 

1. 3783IS4 

-1.3720375 

50.CCOOO 

1.7905658 

-1*5952718 

60.0CCOS 

1.6334237 

•1 *4105702 

70«9eeo« 

l.79C3w96 

-1.5780757 

ac.ooooc 

1.54  34386 

- 1 • 720489? 

9o.ecco.: 

1.9504764 

-1*8364835 

ioc .ooaoc 

2.2687969 

-1*928468'" 

i ic.ccgo: 

2. 49843C4 

-1*9933375 

I20.0C0CC 

2.595o801 

-2*0308568 

I3C  •  C  w  w  C  C 

2. 8722792 

-2*6980584 

I4C.OOOCC 

3. C 7C7277 

-3*2683614 

15C  .:C0O;. 

3.21 29558 

•3.6334579 

160. jccc; 

3.2414393 

-3*7313687 

I70.CC0C: 

3.3894143 

-3*8233167 

I8P.P000C 

3.4477171 

-4.1514867 

190.CCCU' 

3.8501138 

•*.3373303 

200  . ?COC . 

4.C0C7235 

-4.3298462 

2IO.Ci.OOC 

3*8944754 

-4*1441451 

220.:'.«5C 

3.6192376 

-3*6‘iS67C3 

23r  .CuO.f 

3* 26 1 0633 

-J.3SSS*;? 

240  *<i  3CC0 

3*0781179 

-2.8  390 1 *6 

2SC.JOCC0 

3*2349319 

-2*7411546 

26C  .  '? 0  000 

3.36233.4 

-2.6C  707^3 

272.J000: 

3.4644568 

”2 . .6.91 C  * 

260*  :OCCO 

3.5456268 

-2*9266416 

290.C2C0C 

3.60  94345 

-2*980238' 

soe.oaojc 

3*6499451 

-3.029.7.7 

312.CCSJC 

3*6660808 

-3*0749504 

J2C .C  3C0C 

3*6544473 

•3*1 160452 

ijo.ccco: 

3*61 233C3 

-3. IS696C5 

34C .OCCCC 

3*5458596 

-3*1959827 

JSO.POGOC 

3*4329017 

-3*2311901 

360.:coc: 

3. J3S6C66 

-3*2626109 

370 • 39000 

3*196376/ 

-3*2902566 

36C  .£  OOCw 

3.  .*379597 

-3*3141514 

39-?.  300C0 

2*8586736 

-3*3343126 

4CC. >0300 

2*6675712 

-3*35501 1C 

413.C00CC 

2*5990867 

•3*3748913 

420.79 oc: 

2* 7644484 

-3*3914393 

430 . JOCOC 

2*9166189 

440*  30000 

3.C730952 

-3*4146634 

450. COCO C 

3*2146007 

-3. *21 4090 

460. 00020 

3.33.19*0 

-3*4249616 

470*0900C 

3*4866225 

-3*5397916 

4BC.OOOCC 

3*6127973 

490.C30C: 

3*7326111 

-3*5850275 

SOO.COOOv 

3.  63.7307 

-3.5465646 

910.30CCC 

3*9468405 

-3*6634175 

520*00000 

4*0432333 

-3*4270800 
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930.00000 

4.1256196 

-3.6237302 

9*0.00000 

4.21 30045 

-3.6767179 

990.00000 

4.2669452 

-4.07773*6 

900.00000 

4.3506643 

-4.2474*42 

970.0C000 

4.4250706 

—4 • 3765202 

900.00000 

4.7216372 

-4.4899491 

990.00000 

4.9767236 

—4.5945634 

000.00000 

5.2252969 

-4.6959245 

oio.ooooc 

5.4592669 

-4.616*653 

020.00000 

5.6626926 

-4.9521055 

030.00000 

5.9193302 

-5.0996525 

040.00000 

6.1351799 

-5.2737662 

090.00000 

6.3670666 

-5.4526036 

000.00000 

6.5542150 

-5.6260 30C 

070.00000 

6.7356499 

-5.000*25* 

oao.ooooc 

6.6506564 

-6.1279792 

090.0000C 

6.9155640 

-6.4259195 

700.40000 

6.9157765 

-6.6547459 

710.00000 

6.6361212 

-6.6066060 

720.00000 

6.6605494 

-6.6793806 

730.0000C 

6.71 13206 

-6.6792772 

740.00000 

7.0420097 

-6.8001143 

790.00000 

7.3435109 

-6.6503145 

760.00000 

7.6167054 

-6.7043719 

770.00000 

7.6647350 

-6.9075*6* 

700.00000 

6.1104016 

-7.0920961 

790.00000 

6.3164675 

-7.2*70399 

000.00000 

6.4966007 

-7.612602* 

010.00000 

8.6319956 

-7.5*52302 

020.00000 

8.7626024 

-7.6113077 

030.00000 

8.6496922 

-7.6909260 

040.00000 

6.9005405 

-7.7316601 

090.00000 

8.9416541 

-7.9379763 

060.00000 

6.9306807 

-7.960*20* 

070.00000 

6.9139564 

-6.0390376 

000.00000 

6.6609697 

-6.0676199 

090.00000 

6.7759501 

-6.0630035 

900.00000 

6.6621544 

-0.090039* 

910.00000 

6.5376770 

-0.075*725 

920.00000 

6.4046476 

-0. 0601012 

930.00004 

6.2147046 

-0.021010* 

940.00000 

6.1023047 

-7.9620142 

990.00000 

0.0158021 

-7.92*0157 

960.00000 

7.9092662 

-7.059*635 

970.00000 

7.7910631 

-7.786076* 

900.00000 

7.6672070 

-7.7011219 

990.00000 

7.5651563 

-7.62156*0 

1000.00000 

7.4660490 

-7.511736* 

PROCRAM  RAN  TO  COMPLETION 
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APPENDIX  C 

LISTING  OF  THE  GENERAL  PROGRAM  AND  THE  PROGRAM 
WITHOUT  THE  S-C  4020  SUBROUTINES 
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Lilting  of  General  Program 


SRESTORE 

$*  CHARACTRON  OUTPUT  TA®E  ONB9 


SATTACH 

*AS 

SEXECUTE 
*IBJ0B 
SIBFTC  SNERD 
C 


B9 

SYSLB4 
IB  JOB 
MAP.F IOCS 
L I ST.NODECK 

PROGRAMMER - ANTHONY  MELODIA - 

*  APPLIED  MATHEMATICS  LABORATORY- 


MELODIA 


C- 

C- 

C 

r 


-DAVID  TAYLOR  MODEL  BASIN*  WASHINGTON  D.C. - 

- .TELEPHONE  NUMBER— AREA  CODE  301*  995-1514 - 


C 

C 


DIMENSION  Z  (  2000  )  *T  (  2000  )  »S  (  2000  )  »S2ND(  2000  >  .XOMEGA! 2000  )  »X  (  2000  )  • 
1XDOT ( 2000 ) .ALPHA ( 25 )  .OSC IL ( 1000)  *VEL ! 1000 ) *RESID(  1000 ) 

EQUIVALENCE  (RESID*X(1) ) * ( VEL *X ( 1001 )  ) 

NAMELIST/ZAZA1/Z/ZAZA2/T.TZERO.  H*  X0*N» 

1XD0T0. ISZ/ZAZA3/IP1.IP2. IP3.IP4.IP5/ZAZA4/ALPHA.NALPHA 
WRITE! 6*900  I 
P I *3*1415977 
READ ( 5  *ZAZA1 ) 

READ ! 5  *ZAZ  A2 ) 

READ ( 5  *7AZ A3 ) 

READ ( 5  »ZAZA4 ) 

I F ( I P3  )  50. 51 .50 

50  IP1*0 
IP2*0 

51  CONTINUE 

< - 


delete  If  Z's  are  generated  by  the  program 


insert  loop  here  if  Z's  are  to  be  calibrated  or 
generated  by  the  program,  see  listings  of 
sample  problems  1  ft  2 


-  WRITE  AND  PLOT  INPUT  DATA - 

DO  100  1*2. N 
100  T!I)*T!I-1)+H 

I F ( I SZ  )  104.105 .104 

104  WRITE! 6.901 ) 

GO  TO  106 

105  WR!TF(6.90?) 

106  CONTINMF 

WRITEI6.903)  !Z!I).T(I).I*1.N) 

ZMIN*Z 1 1 ) 

ZMAX-ZU) 

DO  110  1*2. N 
IF(ZMAX-Z! I ))108» 109.109 

108  ZMAX*Z  !  I ) 

GO  TO  110 

109  IF(ZMIN-Z! 1)1110. 110. 107 

107  ZMIN«Z!i) 

110  CONTINUE 

CALL  CAMRAV!15) 

?N«N 

DX«H#ZN/10. 

DY*(ZMAX-ZMIN)/10. 

CALL  GRIDlV(l.T(n.T(N).ZMlN.ZMAX.DX.DY.l.l.l.l#6#6) 
CALL  APL0TV1N.T .Z.1.1.1.42.IERR) 

IF( I  ERR) 115.1 16. 115 

115  WRITE(6.904)  I ERR 

116  CONTINUE 


NN«N-1 

DO  120  1*1. NN 

CALL  LINEV(NXV!T! I)).NYV!Z! I )  ).NXV!T!  I  +  1))*NYV(Z(  1  +  1))) 
120  CONTINUE 


100 


Cl  c» 


CALL  PRINTVI-15.15HTIME  IN  SECONDS.452.6 > 

IF(ISZ)121. 122*121 

121  CALL  APRNTVI0.-14.-23.23HFOUNDAT  ION  ACCELERATION.4.696) 

GO  TO  123 

122  CALL  APRNTVIO.-14.-19.19HFOUNDAT  ION  VELOCI TY. 4.664 ) 

123  CONTINUE 
C 

C - COMPUTE  SIN)  AND  S2NOIN)  FOR  EACH  TIME «  T - 

sm»2(2)-zm 

S2NDI1  )  =  Z(3)-2.*Z(2)+Z( 1 ) 

K*N-1 

00  200  I=2*K 
5111*2(1+11-2111 

200  S2NDI I  )*Z( I+1)-2.*Z(  I )+Z( I-l) 

WRITE (6.91 7)  (I*Sm»S2ND(l  ).I*1.N) 

C 

190  READ  I  5.920 )  FRE01 .FRF02.DELTAF 

I  FI ABS I FRE01-FRE02 )  -• IE-20 )  801 .801.191 

191  CONTINUE 

IF! IP3 )  192.195.192 
19?  TEMP*FRFQ2-FRFQ1 

I F ( AMODI TEMP .DEL TAF ) -DEL  TAF/2.)193.194.194 

193  NOFRFO*TEMP/OFLTAF 
GO  TO  195 

194  nofreo=temp/deltaf+i.oooi 

195  CONTINUE 

C— • - DAMPED  OR  UNDAMPED.  IF  UNDAMPED  FIND  NUMBER  OF  RESIDUAL  POINTS 

DO  800  IXM.NALPHA 

N  IN*  1 

INI-0 

FRF0-FRE01 

I F I  ABS  I  ALPHA I  IX) >-. 1F-09 )201 .201 .202 

201  IALPHAM 
GO  TO  2?.? 

?0?  I  alpha*? 

2?2  GO  TO  (703.207 ) .IALPHA 

203  TT-TZERO-TIl) 

IF! AMODI TT.H)-H/2. >205.204.204 

204  KH«TT/H+1. 000001 
GO  TO  207 

205  KH*TT/H+. 000001 
207  CONTINUE 

- COMPUTE  CONSTANTS - 

OMEGA»FREO*2.*PI 
RADCAL*SQRT( 1. -ALPHA I  IX ) **2 ) 

OH-OMEGA++H 
P*OMEGA«RAOCAL 
FX-EXPI-ALPHAI IX)«OH) 

COSPH*COS(P»H) 

SINPH*SIN(P++H) 

A20H*2.#ALPHA( IX) /OH 
RECOH-l./OH 
ALPHA2*ALPHA( I X ) **2 
X12=1.-?.*ALPHA2 
FACT1-FX++SINPH/RADCAI 
FACT2«FX*C0SPH+FACT1*ALPHA(  IX) 

IF(ABS(FRFO)“. 1F-20 ) ??0  »  2?0 .221 
220  IFREO*) 

GO  TO  2?3 
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\J  \JSJ 


221  IFREQ*2 

223  I F  ( ISZ  >300.208.300 

208  GO  TO  (400.209). IFREQ 

209  CONTINUE 

- COMPUTEXOME.GAlNI.XDOf  FOR  VELOCITY  INPUT - 

GO  TO( 224.225 ) .NlN 
724  XOMpGA ( 1 ) ®QMEGA*X0 

FACT3*RFCOH*( 1.-FACT2) 

FACT4*Rt‘C0H*  (  •  5-2 •* ALPHA  (IX)  #RFC0H+EX* (  (  . 5+2.*ALPHA  (  I X  ) *RECOH  ) 

1 *C0SPH- ( XI 2*RFC0H- ALPHA (  I X ) / 2 • ) *SI NPH/RADCAL  > ) 
FACT5*EX*C0SPH-ALPHA( IX)*FACT1 
FACT6=REC0H«FACT1 

FACT7*(RECOH*«2>»(  l.-FX*(C0SPH+(  ALPHA (  IXJ+OH/2.  > »S I NPH/RADCAL  )  ) 
XDOT ( 1 ) *XDOTO 
DO  210  1*2. N 

XOMEGA (  I  )=X0MEGA( I-l)*FACT2+XDOT( I-l >*FACT1-S< 1-1 )  *FACT3-S2ND ( 1-1) 
1*FACT4 

2 1C  XDOT ( I )*-XOMEGA( 1-1 )*FACT1+XD0T (  I  —  1 >*FACT5-S( 1-1 ) *FACT6-S2ND( 1-1  )* 
1FACT7 

GO  TO  756 
225  CONTINUE 

<RO*XMX+XOOK 
<RA*XMN+X I CK 
DO  235  1*2.5 
VA*I-1 
TA*VA#.2*H 
COSPT*COS(P*TA) 

SI NPT  =  S IN ( P*T  A ) 

EXT*EXP(-ALPHA( IX ) #OMEGA*TA ) 

FACTR1*EXT*(C0SPT+ALPHA(  IX  ) *S I NPT/RADCAL  ) 

FACTR2*EXT*SINPT/RADCAL 
FACTR3=( 1.-FACTR1 )»REC0H 

FACTR4*(TA/H-.5-A20H+( ALPHA (  IX >/2.-X12*RECOH)»FACTR2+( .5+A20H) 
1*FXT*COSPT)*RFCOH 

XMAXO  -XOMEGA ( KRO )*FACTR  1 +XDOT  (KRO  >  *FACTR2-S ( KRO » #FACTR3 

l-S2ND(KRO)*FACTR4 
I F ( XMAX-XMAXO  J229.23O.230 

229  XMAX*XMAXO 

230  XMINO  «XOMEGA(<RA)»FACTR1+XDOT(XRA)»FACTR2-S(KRA)*FACTR3 
1— S2ND( KRA )*FACTR4 

IF( XMIN-XMINO  >235.235.232 
23?  XMIN-XMINO 
235  CONTINUE 
NIN*1 

211  IF( IP1I500. 212.500 

212  IF( I P7 >550. 21^. 550 

213  IF( IP3  >600.214.600 

214  IF(IP4>650.215.650 

215  IF( IP5)700.785.700 

- COMPUTE  XOMEGA(N).  XDOT(N).  FOR  ACCELRAT ION  INPUT - 

300  CONTINUE 

GO  TO( 350.305). I FREQ 
305  CONTINUE 

GO  T0( 308.325 ) . N I N 
308  XOMEGA( 1 )»0MFGA«X0 
FCOSl«l.-FX*COSPH 
DFACT3=(1.-FACT2) /OMEGA 

DFACT4»(1.-A20H*EC0S1-X12*FACT1#REC0H)/0MEG« 
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DFACT5=!-4.#ALPHA<  I X )-? <  2.*  !  1 . -4.*ALPHA2  )  /0H-2.*ALPHA t  IX  )  >*EC0S1+ 

1 1X12+2. *ALPHA!  IX)*!3.-4.*ALPHA2)/OH)*FACT1)«RECOH/<2.*OMEGA) 
DFACT6* ( 1.-EC0S1 ) -FACT 1+ALPHA ! IX  ) 

DFAC T7= FAC T1 /OMEGA 
DFACT8=(1.— FACT2) +RECOH /OMEGA 

DFACT9*!2.-!1.+4.»ALPHA!  IX )  »REC0H)*EC0S1-(2.»X12*REC0H-ALPHA!  IX)  ) 
1 *FACT 1  >*RFC0H/(2.*0MFGA) 

XDOT ( 1  )  =XDOTO 
DO  310  1*2  »N 

XOMEGA! I )=X0MFGA< I-l )#FACT2+XD0T ! I-l )*FACT1-Z! I-l )*DFACT3-S! I-l ) 
1*DFACT4-S2ND< t-1 )#DFACT5 

310  XDOT  ( I  )=-XOMEGA( I-l  )*FACT  1 +  XDOT  (  I-l )  +DFACT6-Z  ( 1-1 )  +DFACT7 
1  -S!I-1)»DFACT8-S2ND!  I-l >*DFACT9 

GO  TO  756 
325  CONTINUE 

<RO*<MX+KOOK 
KRA*KMN+KICK 
H2*H**2 
DO  335  1*2*5 
VA* I —1 

ta*va#.?#h 
COSPT=COS ( P*TA ) 

SINPT*SIN(P*TA) 

FXT*EXP( -ALPHA! I X  )#OMEGA*TA > 

FACTR1*EXT*( CO SPT+ ALPHA!  IX  )  *S I NPT/RADCAL  ) 

FACTR2*EXT*SINPT/RADCAL 
FACTR3* ( 1 .-FACTR1 ) /OMEGA 

FACTR4* !  TA/H-A20H* !  1  .-EXT+COSPT )  -X 1 2»FACTR2*RECOH )  /OMEGA 
FACTR5*(TA**2/H2-TA/H-!2.*(  l.-4.*ALPHA2)»REC0H«»2-A20H)* 

1  (1.-EXT«C0SPT)  +  (X12*REC0H+2.*ALPHA(  IX)*(3.-4.*ALPHA2)*RECOH**2) 
2*FACTR2)/!P.*0MEGA) 

XMAXO*XOMEGA(KRO)»FACTR1+XDOT!KRO)*FACTR2-Z(KRO)*FACTR3 
l-«i!KRO)#FACTR4-S2ND!KRO)*FACTR5 
IF! XMAX-XMAXO  )329. 330*330 

329  XMAX*XMAXO 

330  XMINO*XOMEGA!<RA)»FACTR1+XDOT(KRA)*FACTR2-Z!KRA)*FACTR3 
1— S(KRA)#FACTR4-S2ND!KRA)*FACTR5 

IFJXMIN-XMINO  1335*335*332 
332  XMIN*XMIN0 
335  CONTINUE 
NIN»1 
GO  TO  211 

- FRE01=ZER0  FOR  FOUNDATION  ACCELERATION— - 

350  GO  TO(352*365) *NIN 
352  XOMEGA ( 1 ) *X0 
XDOT ! 1 ) *XDOTO 
H?*H«»2 
DO  360  I *2  *N 

XOMEGA!  l)«XOMEGA!  I-l  )+H*XDOT<  I-l  1-Z  C  I-l  )*H2/2.-S<  I-l)*H2/6.+ 
1S2ND! I-l)*H2/24. 

360  XDOT! I  )*-Z<  I-1)»H-S!l-1 1+H/2.+S2ND! I-1)*H/12. 

GO  TO  756 
365  CONTINUE 
NIN*1 
GO  TO  2)1 

- COMPUTE  X •  XDOT  FOR  OMEGA-ZERO - VELOCITY  INPUT 

400  GO  TO! 402*425) »NIN 
402  XOMEGA! 1 ) *X0 
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XD0TC1  )*XDOTO 
DO  410  I*2»N 

X0MEGA(I)*X0MEGA(I-1)+XD0T< I-1)*H-S( I-1)*H».5+S2ND< I-1)*H/12. 
410  XDOT ( 1  )*XDOT ( 1-1 >-S( l-l ) 

GO  TO  756 
425  CONT1NUF 
NlN«l 
GO  TO  211 

- PLOT  X  VS.  T  AND  XDOT  VS.  T - 

500  GO  TO ( 505*510 ) » I FREQ 

505  DO  506  1*1 »N 

506  X ( I ) “XOMEGA ( I ) 

GO  TO  518 

510  DO  515  1*1 *N 

515  X! II *XOMEGA»I ) /OMEGA 

518  XMX*Xm 

XDMX*XDOT ( 1 ) 

XMN*X(1) 

XDMN*XDOT ( 1 ) 

DO  550  1*2. N 
1F(X(I)-XMX  ) 521 *521 *520 

520  XMX*XUI 

521  IE(X( I )-XMN  )522. 523*523 

522  XMN*X(I) 

523  IF( XDOT ( 1 )— XDMX  >525.525*524 

524  XDMX*XDOT ( I > 

525  IF!XDOT(I)-XDMN  >526.530.530 

526  XDMN*XDOT<I) 

530  CONTINUE 

DY*  (  XMX  rXMN  ,>/10, 

CALL  GRTD1V(1*T(1>.T(N)*XMN*  XMX*  DX.DY* 1 *1  *  1* 1 *6*6 ) 

CALL  APLOTV(N*T*X *1*1* 1.42 »IER) 

IF { I ER >531*532*531 

531  WRITF(6.913>  IER.FREO 

532  NN*N-1 

DO  535  1*1 »NN 

CALL  LINEV(NXV(T( I >  >  »NYV(X( I ) >  »NXV(T( I  +  l>  >.NYV(X( I+l>  >> 

535  CONTINUE 

CALL  PRINTV(-15*15HTIME  IN  SECONDS. 452. 6 > 

CALL  PRINTV(-23*23HFREO»  CYCLES/SEC. *20.3 > 

CALL  LABLVIFRFQ* 60*3*6*1 *4  > 

CALL  PRINTV(-8»8HALPHA*  *20.17> 

CALL  LABLV(ALl>HAf  IX>  *72.17.6.1.1  > 

CALL  APRNTV(0«-14.-11.11HX  RESPONSE»4.600> 

DY* ( XDMX  -XDMN  >/10. 

CALL  GRID1V(1»T(1I»T(N)*XDMN*  XDMX*  DX*DY*1»1*1*1»6*6> 

CALL  APLOTV  (N.T.XDOT*  1 *1 *1 *42 » IERAI 
I F  C I  ERA)  538*539*538 

538  WRITE(6»914)  IERA.FREO 

539  DO  540  1*1 *NN 

CALL  LINEV(NXV(T( I >  >.NYVIXDOT( I ) ) *NXV(T( 1+1) ) .NYVIXDOTI I  +  l >  > ) 

540  CONTINUE 

CALL  PRINTV(-15*15HTIME  IN  SECONDS. 452. 6 > 

CALL  PRINTV(-23»23HFREO*  CYCLES/SEC. *20*3 > 

CALL  LARLV(FREQ*60*3*6»1»4> 

CALL  PRINTV(-8*8HALPHA*  .20*17) 

CALL  LABLV( ALPHA! IX) *72*17*6*1*1 ) 

CALL  APRNTV(0*-14»-22*22HVELOCITY  RESPONSE  XPOT»4.688> 

GO  TO  212 
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- PLOT  Y( ACCELERATION)  VS.  TIME - 

550  CONTINUE 

GO  TOI213.551 ) » IFREO 

551  CONTINUE 
560  00  56?  1=1 #N 

562  X( I )=-XOMFGA(I)#OMFGA 

XMX«X( 1 ) 

XMN=X(1) 

DO  566  I =? »N 
I F ( XMX  -X(  I  ))563. 564.564 

563  XMX=X(I) 

564  I F ( XMN  -X( 1)1566.566.565 

565  XMN=X(I) 

566  CONTINUE 
OY=(XMX  -XMN  )/10. 

CALL  GRID1V(1.T(1).T(N).XMN.  XMX.  DX . DY . 1 . 1 . 1 ♦ 1 .6 .6 ) 

CALL  APLOTVIN.T.X.1.1 .1.42.IERR) 

IF( IERBJ570.571.570 

570  WRITEI6.919)  IERB.FRFO 

571  NN=N-1 

00  575  1=1. NN 

CALL  LINEV(NXV(T< I ) ) .NYV(X< I ) ) »NXV( T( 1  +  1 ) ) .NYVI  XI 1  +  1) ) ) 

575  CONTINUE 

CALL  PRINTVI-15.15HTIME  IN  SECONDS. 452. 6 ) 

CALL  PR  INTV ( -23.23HFREO®  CYCLES/SEC. .20.3 ) 

CALL  LABLV( FRE0.60.3«6.d.4) 

CALL  PRINTVI-8.8HALPHA*  .20.17) 

CALL  LARLV ( ALPHA ( IX). 72 .17. 6. 1*1) 

CALL  APRNTV 1 0 .-14.-34.34HMINUS  PSEUDO  RESPONSE  ACCELERATION* 
14.844) 

GO  TO  214 


-  FOUR  COORDINATE  GRID - 

600  CONTINUE 
INI-1+INI 
OSCILI INI )=FREO 

VFL ( INI  )>AMAX1 (ABS(XMAX) .ARS(XMIN) ) 

GO  TO  (601 .602). I ALPHA 

601  RES  I  D(  INI  )=SQRT(XOMEGA(KH)*»2+XDOT(KH)*»2) 
I F ( RES  I D( INI )-. IE- 19 1604.604. 602 

60?  IF(FRE0-.1E-19)604. 604.603 

603  IF(VEL( INI )-.l £-19)604.6 04. 60 5 

604  NOFREO=NOFREQ-1 
IN  I* INI-1 

605  I F ( ABS ( FREQ-FREQ2 )-• IE-09 ) 606 .606.214 

606  CONTINUE 

CALL  SMXYV(l.l) 

FRQSML=OSC!L( 1 ) 

VFLSML=VELU) 

RESSML-RFSIDI1) 

DO  616  I  =2 .NOFREO 
IFCOSCILCD-FROSML  1611.612. 612 

611  FROSML=OSCIL(  I ) 

612  IF(VEL(  I  )-VELSMD613. 614.614 

613  VELSML*VEL(  I ) 

614  I F(RESIDU)-RESSML) 615.615 .616 

615  RESSML*RESID(  I ) 

616  CONTINUE 
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0RDSML=AMIN1 ( VELSML  .RFSSML ) 

- TO  FIND  THE  LIMITS  FOR  THE  14X3)  CYCLES - 

I F ( ALOGIO ( ORDSML  )  (620.621.621 

670  LOGORD=ALOG10(ORDSML)-1  . 

GO  TO  672 

671  LOGORD=ALOG10(ORDSML) 

622  I F ( ALOGIO ( FRQSML ) >623.624.624 

623  LOGFRQ=ALOG10 ( FRQSML ) -1 • 

GO  TO  625 

624  LOGFRQ=ALOG10 ( FRQSML ) 

625  ORDSML=10.**LOGORD 
ORDLG»ORDSML*10.«*4 
frqsml=io.**logfrq 

FRQLG*FRQSML*10.«*3 

c - plot  log-log  grid  for  velocity  vs  frequency - 

CALL  GRID1 V( 1 .FRQSML .FRQLG .ORDSML »ORD LG *1.0»1.0*l»l»l»l»-2»-2 ) 

C - TO  FIND  LARGEST  1X10  TO  THE  PTH  POWER  LINE  FOR  X - 

omst'rt=?.*pi*frosml 

TEMP=AL0G10(0RDLG/0MSTRT I 
IF(TEMP)627.628.628 

627  LOP*TEMP-l. 

GO  TO  629 

628  LOP«TEMP 
679  01=10. **LOP 

V1=0MSTRT*D1 

WIG*10.*ORDSML 

XMARGN=NXV ( FRQLG  )-N XV ( FRQSML ) 

YMARGN=NYV(0RDLG)-NYV(WIG) 

SLOPF=YMARGN/XMARGN 

IXC0R*4+NXV(FRQSML) 

CALL  PRINTVI-1 » 1HD » IXCOR  »N YVC  VI ) ) 

CALL  PR  INTV ( -2  »2HD=  *900*1 7 ) 

CALL  LABLV ( D1 *924 .17.-2.1.3 ) 

C - DRAW  LINES  UP  FROM  Vl - 

VIS*0. 

DO  632  1*1 .9 
VIS«VIS+V1 

I F(ORDLG-VIS >647*630.630 
630  XTCH«NYV(ORDLG)-NYVIVTS) 

I TCH *X TCH /SLOPE 
I TCH*NXV ( FRQSML ) ♦ I TCH 
XITCH*NYVI0RDL6) 

CALL  LINEV  <NXV<  FRQSML ) »NYV < Vi S ) » I TCH.K I TCH ) 

632  CONTINUE 

C - DRAW  LINES  FROM  VI  DOWN - 

647  DO  638  J«l*20 
LAL-LOP-J 
lan-lal+i 
dlan«io.**lan 

DLAL*10,*»LAL 
DO  638  1*1*9 
FFT1*I 

VIS«OMSTRT*DLAN-OMSTRT*DLAL*EFTl 
IFIVIS-ORDSML >639.633.633 
633 ,KITCH*NYV(0RDLG1-NYV(VIS> 

IYMARG*YMARGN 

I F I  I YMARGHCITCH)635» 634.634 
634  Z I TCH=K I TCH 
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ITCH*Z ITCH/SLOPE 
GO  TO  637 

635  ITCH*XMARGN 
KITCH*YMARGN 

637  ITCH* I TCH+NXV ( FRQSML ) 

K I TCH*K I TCH+NYV (VIS) 

CALL  LINEV  ( NXV! FRQSML ) »NYV ( V I  5 ) » I TCH »K I TCH ) 

638  CONTINUE 

C - DRAW  FROM  THE  RIGHT  ORDINATE  THE  REMAINING  X-LINES - 

639  OMLAST=2.*PI*FRQLG 
LEFT  2* J 
LFFT1*FFT1 

DO  641  J*LEFT2.30 

lal*lop-j 

lan*lal+i 

dlan*io.#*lan 

dlal*io.**lal 

DO  640  I*LEFT1 .9 
FFT1*I 

VIS*0MLAST*DLAN-0MLAST»DLAL»EFT1 
I F ( V I S-ORDSML 1642*636*636 

636  KITCH=NYV!VlS)-NYV<ORDSML) 

CAPT  *K I TCH 

ITCH*CAPT /SLOPE 
ITCH*NXV(FR0LG)-ITCH 
K I TCH=NYV ( ORDSML ) 

CALL  LINEV  <NXV< FROLG » .NYV! VIS > ♦  ITCH*  KITCH) 

640  CONTINUE 
I.PFTI*! 

641  CONTINUE 

642  CONTINUE 

C - .—PLOT  POINTS - 

CALL  APLOTVINOFREO.OSCIL.VEL.1.1  .1.44.LIT) 

IF!LIT)643.644.643 

643  WRITE!6.922)  L  I T  .ALPHA ( I  X } 

644  GO  TO  1648.646) » I ALPHA 

648  CALL.  APLOTVINOFREQ.OSCIL.RESIO* 1.1. 1.38. KIT) 

CALL  PRINTVI-28.28HOOO  *RE5 IDUAL  SHOCK  SPECTRUM.NXVI FROSML ) . 1006 ) 
IF(KIT) 64 5. 646.645 

645  WRITEI6.923)  K I T  .  ALPHA < I X > 

646  CONTINUF 

CALL  PRINTVI-15.15HFRFQUENCY  ICPS).452.6) 

CALL  APRNTVI0.-14.-8.8HVEL0CITY .9.576 ) 

CALL  PRINTVI-8.8HALPHA*  .20.17) 

CALL  LABLVI ALPHA! IX ) . 72. 17 .6 . 1 . 1 ) 

CALL  PRINTV(-19»19H*»*  ■ SHOCK  SPECTRUM. NXV! FRQSML ). 1015 ) 

- »  GENERATION  OF  ACCELERATION  GRID - 

-  TO  FIND  A 1 - 

TEMP«ALOG10IORDSML*OMSTRT) 

IF! TEMP  1850.851.851 

850  LOP*TEMP 
GO  TO  853 

851  LOP-TEMP+1, 

853  A1»10.**LOP 

Vl *A1 /OMSTRT 

CALL  PRINTV!-1»1HA. IXCOR.NYV!Vl ) ) 

CALL  PRINTVl-2 .2HA* .900.4 ) 

CALL  LABLV!A1 .924.4.-2.1  .3) 

C - DRAW  LINES  DOWN  FROM  A1 - 

K I TCH*NYV ! ORDSML ) 
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DO  856  f=l«9 
FFT1*I-1 

VlS»Vl*(t.-.1*EFTl  I 
IF  (ORDSML-Vl  5)855  *855  ♦857. 

855  XTCH*NYV(VIS)-NYV<ORDSML) 

ITCH«XTCH/SLOPE 

I TCH* I TCH+NXV ( FROSML ) 

CALL  LINEV  ( NXV< FROSML ) tNYV ( V I S )♦ 1 TCH ♦  KITCH) 

856  CONTINUE 

C - DRAW  LINFS  FROM  A1  UP - 

857  DO  870  J=1#20 
lan«lop+j-i 

ALAN*10.*#LAN 
DO  870  I=?tlO 
FFT] *1 

VI5«ALAN/0MSTRT*EFT1 
I F( V I S-ORDLG )858t858t871 

858  KITCH*NYV(VIS)-NYV(ORDSML) 

TYMARG=YMARGN 

I F ( I YMARG-K ITCH)862i061t861 

861  Z I TCH=K I TCH 

I TCH*Z I TCH /SLOPE 
KITCH=NYV(0RDSML> 

GO  TO  868 

862  ITCH«XMARGN 

K I TCH*NYV (VIS) -I YMARG 
868  I TCH* I TCH+NXV I FROSML) 

CALL  LlNFV  < NXVI FRQSVL I tNYV I VI S ) ♦ I TCHf  KITCH) 

870  CONTINUF 

C - -DRAW  FROM  THE  RIGHT  ORDINATE  ^REMAINING  A-LINES 

871  LFFT?=J 
LFFT1=EFT1 
K!TCH=NYV(OROLG) 

DO  875  J*LFFT?t80 
LAN*L0P+J-1 

alan*io.»*lan 

DO  874  I*LFFT1.10 
EFT1-I 

vis«alan/omlast*efti 

I F ( V  I S-ORDLG )872t872#876 
87?  CAPT*NYV(0RDLG)-NYVIVIS) 

I TCH *C APT /SLOPE 
ITCH*NXV(FROLG)-ITCH 

CALL  LINEV  ( NXVI FRQLG ) tNYVl VIS  > ♦  iTCHtKITCH) 

874  CONTINUE 
LFFT1*2 

875  CONTINUF 

876  CONTINUF 
GO  TO  214 

- # —  SHORT  FORM - - — 

650  CONTINUF 

I F  (  ABS ( FREO— FRE01 ) -• 1 F-04 )651f651t65? 

651  WRITE(8t91*)  ALPHA(IX) 

65?  GO  TO(658t655) tlALPHA 
655  X0MAX*AMAX1 (ABS(XMAX)t ABS(XMIN) ) 

GO  TO ( 668 ♦ 669 ) ♦ I FREO 
658  RMAX*SQRT(XOMEGA(KH)*»2+XDOT(KH)»«2) 

X0MAX*AMAX1 ( ABSI XMAX ) ♦ ABS ( XM IN ) +ABS I RMAX ) ) 

GO  T0(681t678)#IFRE0 
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668  WR I TFI 8  *906 )  FRFQ.XOMAX 
go  TO  215 

669  WRITF( 8.905 1  FREO.XOMAX 
GO  TO  215 

678  WRITE! 8*907 )  FRFQ. XOMAX »RMAX 
GO  TO  215 

681  WRITE<8.908)  FREO.XOMAX. RMAX 
GO  TO  215 

- long  FORM - 

700  CONTINUF  * 

IF( ARS! FREQ-FREQ1 1 -.IE-04) 701 *701*702 

701  WR ITF! 6 .91 5  1  ALPHA!  IX ) 

70?  GO  TO! 704.705 1 .1 ALPHA 

704  RmAX*SQRT!X0MFGA!KH1#*2+XD0T!XH)*»2) 

GO  TO! 745.742 1 .IFRFO 

705  GO  TO(7?6»727)»IFREQ 

7?6  WRITE! 6*916 1  FREQ.XMAX.XMIN 
GO  TO  785 

727  WRITE16.910)  FREQ.XMAX.XMIN 
GO  TO  785 

742  WRITE!6«910)  FREO»XMAX  »XM IN  *RMAX 
GO  TO  785 

745  WRITE!6,911>  FREQ «XM AX .XMIN • RMAX 
GO  TO  785 

- PRELIMINARY  CALCULATIONS  FOR  MIN.  AND  MAX.  RESPONSE 

756  GO  TO  1758.7571. I ALPHA 

757  NX«N 

GO  TO  759 

758  NX*KH 

759  XMAX-XOMFGA11  ) 

<MX«1 

XmIN*X0MEGA!1  1 
<MN»1 

DO  760  I >2 .NX 

IF1XMAX-XOMEGA! 1 1  )761 .762.762 

761  XMAX«XOMEGA! I  1 
<MX«I 

762  1F1XMIN-XOMEGA! 1 1  1760.765.765 
765  XMlN«XOMEGA! I ) 

XMN»I 

760  CONTINUE 

I F ( KMX-1 1771*770.771 

770  KMX*KMX+1 

771  IF1KMX-NXI773. 772.771 

772  <MX«<MX-i 

775  1F0CMN-11775.774.775 

774  ICMN«XMN+1 

775  1FIXMN-NX1777.776.777 

776  <MN«XMN-1 

777  CONTINUF 

!F(XDOT!XMX-l) *XDOT ! XMX 1 1778.778.779 

778  XOOX«-l 
GO  TO  781 

779  XOOX-O 
781  CONTINUF 

I F ( XDOT 1 KMN-1 1 *XDOT l KMN) 1783.783.784 
786  XICX»-1 
GO  TO  950 
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784  <  FfK=0 

950  WTNs? 

GO  TO  707 

- EITHER  RESTART  WITH  NEW  FREQ.  OR  RESTART  WITH  NEW  ALPHA  OR  STOP 

785  CONTINUE 

I F ( ABS ( FRE0-FRE02 )  -  ,  1 F-03 )800»755*755 

755  FREQ=FREQ+OFLTAF 
GO  TO  ?07 

800  CONTINUE 

GO  TO  190  ^ 

801  WRITE! 6*909 ) 

CALL  framfv 

STOP 

900  FORMAT 1 1  HI »49X  .22HSHOCK  SPECTRUM  PR0GRAM//50X  *2 1HAPPL I  ED  MATH  LAB 
10TMR,///) 

901  FORMAT ( 40X » 38HFOUNDAT ION  ACCELERATION  VS  TIME! INPUT ) // 1 8X  » 23HFOUND 
1ATI0N  ACCELERATION, 36X.15HTIME  IN  SECONDS) 

90?  FORMAT !  42X  »35HFOUNOAT ION  VELOCITY  VS  TIME  ! I NPUT  ) //20X » 19HF0UNDAT I 
ION  VELOCITY.35X.15HTIME  IN  SECONDS) 

903  FORMAT! IX. F35. 9.  44X.F20.9) 

904  FORMAT ( 47H  PLOTTING  ERROR  NUMBER  OF  POINTS  OUT  OF  RANGE*. 12//) 

905  FORMAT !15X.F1C.5»29X.F12.7) 

906  FORMAT (10X.11HFIRST  FREQ* . F  10.5 . 7X . 13HMAX  RESPONSE* »F 12. 7/ ) 

907  FORMAT! 15X.F10.5.2! 29X.F12.7) ) 

908  FORMAT!  10X.11HFIRST  FREO* » F10 .5 »7X » 13HMAX  RFSPONSE=  .F12. 7 . 12X . 
117HMAX  RES  RESPONSE*. F12. 7) 

909  FORMAT ( // /38X . 25HPR0GRAM  RAN  TO  COMPLETION) 

910  FORMAT !10X.F10.5»3!18X.F12*7) ) 

911  FORMAT I5X.11HFIRST  FREO* .F10.5.9X.9HMAX  RESP* .F 12. 7 » 10X .9HMIN  RESP 
1*.F12.7,7X. 15HMAX  RESID  RESP*.F12.7) 

912  FORMAT! 1H1.55X.10HSHORT  FORM//45X.20HDAMPING  COEFFICIENT*. F10. 8/// 
1 /9X.22HFREQUENCY ( CYCLES/ SEC .)» 19X .20HMAX  RESPONSE  X  OMEGA, 18X.26HM 
2 AX  RESID  RESPONSE  X  OMEGA) 

916  FORMAT ! 59H  PLOTTING  ERROR  FOR  X  VS  T.  NUMBER  OF  POINTS  OUT  OF  RAN 
1GE* . 12 »6H  FREO*. F10. 7//) 

914  F0RMATdX.I2.43H  POINTS  OUT  OF  RANGE  FOR  XDOT  VS  T IME  .FREO* .F10. 7 
1//) 

915  FORMAT! 1H1.55X.9HL0NG  F0RM//45X.20HDAMPING  COEFFICIENT*. F10. 8/7// 
17X. 16HFREQ I  CYCLES /SEC ) ♦ 12X .20HMAX  RESPONSE  X  OMEGA. 10X.20HMIN  RESP 
20NSE  X  OMEGA. 8X.22HMAX  RESID  RESP  X  OMEGA) 

916  FORMAT (4X.11HFIRST  FREO«»F10.5»10X.9HMAX  RESP*.F12.7.10X,9HMIN  RES 
1P*.F12,7I 

917  FORMAT I1H1.18X.1HI.38X.4HS!  I ) »34X .7HS2ND! I )// 
1(17X.!4,34X»F12.7,28X.F12.7) ) 

918  FORMAT! 1H1 .50X ♦ 19H INTERMED I  ATE  VALUES/ 19X, 1H I .35X.9HXOMEGA! I ) , 
131X.7HXD0T! I )// 

2(17X,I3,34X.F1?.7.28X»F1?.7) ) 

919  FORMAT (1X.I2.43H  POINTS  OUT  OF  RANGE  FOR  Y  VS  TIME»FRE0*»F10.7 
1//) 

920  FORMAT ! 3F?0. 8 ) 

921  FORMAT (6F1 5 .6 ) 

92?  FORMAT! IX, 12. 49H  POINTS  OUT  OF  RANGE  FOR  VEL  VS.  FREO  WITH  ALPHA*, 
1F10.8) 

926  FORMAT! IX, 12. 50H  POINTS  OUT  OF  RANGE  FOR  RESIDUAL  PLOT  WITH  ALPHA* 
1.  F10.8) 

END 
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Listing  oi  me  Alternate  rrognm  \nu  o  -v<  w«u  miuucu; 


SEXFCUTF  IRJOR 

SIRJOR  MAP  #F IOCS  MELODIA 

SIRFTC  SNFRD  LIST  .NODFCK 

C  PROGRAMMER - ANTHONY  MELODIA - - — 

C -  APPLIED  MATHEMATICS  LABORATORY - ' - 

C - DAVID  TAYLOR  MODEL  BASIN*  WASHINGTON  D.C. - 

C - TELEPHONE  NUMBER - AREA  CODE  301*  995-1514 - 

C 

C 


C 

C 


DIMENSION  Z (2500 ) • T C 2500 » .5(2500) .S2NDI 2500) *XOMEGA(2500) ♦ 

ixdot (2500 ) .alpha ( 25 ) 

NAMELIST/ZAZA1/Z/ZAZA2/T.TZER0.  H.  XO.N. 

1XDOTO. ISZ/Z AZA3/ IP4. IP5/ZAZA4 /ALPHA. N ALPHA 
WRITE! 6.900 ) 

PI=3. 1415927 


100 

104 

105 

106 


RFAD!5.ZAZA1  )f- 
READ(5.ZAZA? ) 
READI5.7AZA3 ) 
READ  ( 5‘.  Z  AZ  A4  ) 


delete  if  2's  are  to  be  generated  by  the  program 


■ 

)  PLO 


-  WRITE  AND 

DO  100  1*2. N 
T(I)«T(I-n+H 
IF(ISZ)104.105.104 
WRITE! 6.901 ) 

GO  TO  106 
WRITE! 6.902  I 
CONTINUF 

WR  I TE( 6 .903 )  ( Z (  I  )  *T (  I  )  .  I  =  1 .N ) 


#- 


PLOT  INPUT  DATA - 

insert  loop  here  if  Z's  are  to  be  calibrated  or 
generated  by  the  program,  see  listings  of  sample 
problems  1  &  2 


C - COMPUTE  SIN)  AND  S2ND<N)  FOR  EACH  T I  ME  »  T - 

SI1)«Z!2)-Z!1) 

S2NDU )=ZI3)-2.*Z!2)+Z<1) 

K=N-1 

DO  200  1*2. K 
5ll|iZII+l)-Z(I| 

200  S2ND !  I  )=Z! 1+1 )-2.*Z< I)+Z! 1-1) 

WRITE!  6. 91 7)  I  I. SI  I  )  .S2NDIH  .  1*1  .N) 

C 

190  READI5.920)  FREOl .FRE02 .DELTAF 

I F ( ARS ( FREQ1-FRE02 ) -• 1 E-20 )  801.801.191 

191  CONTINUE 

C - DAMPED  OR  UNDAMPED. IF  UNDAMPED  FIND  NUMBER  OF  RESIDUAL  POINTS 

DO  800  I X* 1 .NALPHA 
N  IN*  1 

FREQ*FREQ1 

I F ( ABS ! ALPHA ( IX) )-. IE-09  )201 .201 .202 

201  I  ALPHA* 1 
GO  TO  222 

202  IALPHA*? 

222  GO  TO  !  203. 207). I ALPHA 

203  TT*TZERO-T! 1 ) 

IF(AMOD!TT»H)-H/2. ) 205.204. 204 

204  KH-TT/H+1. 000001 
GO  TO  207 

205  KH*TT/H+. 000001 
207  CONTINUE 

C 

C - COMPUTE  CONSTANTS - 

OMEGA*FREQ»2.*PI 
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RADCAL=*ORT ( 1 .-ALPHA( IX  I **2 ) 

OH=0MFf,A*H 

p=omfga#radcal 

FX=FXP(-ALPHA( IX)«0H) 

COSPH*COS(P*H ) 

S I NPH=  SIN ( P*H ) 

A20H«2.#ALPHA{ IX) /OH 
RFCOH=l./OH 
ALPHA2*ALPHA( IX)**2 
X12*1.-2.*ALPHA2 
FACT1=FX»SINPH/RADCAL 
FACT?=FX»fOSPH+FACTl*ALPHA( IX) 

!F(APMFRFQ)-.lF-2n)??Ct  220*221 

220  f FRFQ= 1 
GO  TO  22? 

221  IFRFO=? 

221  IF!  ISZHOO. 200*100 

208  GO  TO  ( 400 *209 ) * I FRFQ 

209  CONTINIJF 

- COMPUTEXOMEGA(N) *XD0T  FOR  VELOCITY  INPUT - 

GO  TO! 224*225 ) .NlN 

224  XOMEGA ( 1 ) *OHEGA*X0 
FACT1«REC0H»! 1.-FACT2) 

FACT4*RFC0H*( .5-2.*ALPHA( l X )*REC0H+FXM  (  .5+2**ALPHA| IXJ»RECOH) 
1*C0SPH-!X1  2*RFCOH-ALPHAI  IX ) /2.  )*SlNPH/RADCAL ) ) 
FACT5*EX»COSPH-ALPHA( IX)*FACT1 
FACT6*RFC0H«FACT1 

FACT7»!RECOH*«2)*I 1 .-EX* ( COSPH+f ALPHA!  IXI+OH/2* ) *S INPH/RADCAL ) ) 
XOOTI 1 )»XDOTO 
DO  210  I *2  *N 

XOMEGA ( I )*XOMfc'GA( I -1) *FACT2*XDOT  < I  —  1 » *FACT1“S(  I “1 ) *FACT3-S2ND I  1-1) 
1*FACT4 

210  XDOT 1 1 ) «-XOMEGA ( 1-1 I *FACT 1+XDOT !  1-1 )*FACT5-S!  1-1 )*FACT6-S2ND! I-D* 
1FACT7 

GO  TO  -756 

22 5  CONTINUE 
KRO«KMX+XOOK 
KRA«KMN+XICK 
DO  215  I«2*5 
VA* I -1 
TA»VA*.2*H 
COSPT»COS(P*TA) 

SINPT*SIN! P*T  A ) 

EXT*EXPI-ALPHA( IX)*0MEGA#TA) 

FACTR1*EXT* ( CO SPT+ ALPHA! IX ) *S I NPT/RADCAL ) 

FACTR2*FXT*SI NPT/R ADCAL 

FACTR1* ( 1 .-FACTR 1 ) »RFC0H 

FACTR4=(TA/H-.5-A20H+( ALPHA! I X  ) / 2*-Xl?#RECOH > *F ACTR2+ < • 5+A20H ) 
l*EXT»COSPT)*RFCOH 

XMAXO  «XOMEGA(<RO)*FACTRl+XDOT!KRO)»FACTR2-S(KRO)»FACTR3 
1-S2NDIKR0)*FACTR4 
IFIXMAX-XMAXO  >229*210*230 

229  XMAX«XMAXO 

230  XMINO  «X0MEGA( KRA )*FACTR1+XD0T ( KRA) »FACTR2-S( KRA )*FACTR3 

1-S2ND! KRA )*FACTR4 

I F ( XMl N-XM I  NO  >235*235*232 

232  XMIN«XMIN0 
235  CONTINUE 
NIN»1 
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214  IFIIP4)650.215.650 

215  IF! IP5J700* 785*700 
C 

C - COMPUTE  XOMEGAIN)*  XDOTIN).  FOR  ACCELRAT ION  INPUT - 

300  CONTINUE 

GO  T0<350*30«>)  .IFREQ 
305  CONTINUE 

GO  T0I308*325) »NIN 
308  XOMEGA ( 1 ) =0MEGA*X0 
E COS  Is 1 »-EX*COSPH 
DFACT3® ( 1.-FACT2 ) /OMFGA 

0FACT4=(1.-A20H*EC0S1-X12*F ACT  1*REC0H> /OMEGA 

DFACT5=<-4.«ALPHA( IX)-(2.*< l.-4.*ALPHA2)/OH-2.*ALPHA( IX) )»EC0S1+ 
l’l  X12+?»#ALPHA  (  IX)*(3.-4.*ALPHA2) /0H)*FACT1 )*RECOH/(2.*OMECA) 
0FACT6=(1.-FC0S1 ) -FACT1 *ALPHA I  IX) 

DFACT7=FACT1 /OMFGA 
DFACT8= ( 1.-FACT2 ) *RFCOH/OMEGA 

DFACT9=(2.-(1.+4.*ALPHA(  I X ) +RECOH ) +EC0S1- (  2 .*X 1 2*REC0H-ALPHA ( I X ) ) 
1*FACT1 )*RECOH/(2.*OMEGA) 

XOOT ( 1 ) =XDOTO 
DO  310  I =2  »N 

XOMEGA ( I )*XOMEGA< I -1 ) +FACT2+XD0T ( I- 1 ) *FACT 1-Z (I -1 ) +DFACT3-SI 1-1 ) 
1*DFACT4-S2ND(  1-1 )*DFACT5 

310  XDOT ( I ) =-XOMEGA( I -1 ) «FACT1 +XDOT I  1-1 )*DFACT6-Z( 1-1 >*DFACT7 
1  -SI  1-1 )*DFACT8-S2ND» !-l >»0FACT9 

GO  TO  756 
325  CONTINUF 

KRO«KMX+KOOK 
KRA*KMN+K I CK 
H2«H**2 
DO  335  1«2#5 
VA* I -1 
TA«VA#.2»H 
COSPT=COS(P*TA) 

SINPT*S!N(P»TA) 

EXT»EXP( -ALPHA! IX) *OMEGA*TA  ) 

FACTR1=FXT»IC0SPT+ALPHA( IX1+SINPT/RADCAL) 

FACTR2=EXT*SINPT/RADCAL 
FACTR3M  1.-FACTR1 ) /OMEGA 

FACTR4=(TA/H-A20H*(1.-EXT*C0SPT)-X12*FACTR2*REC0H)/0MEGA 
FACTR5=(TA**2/H2-TA/H-(2.*I l,-4.*ALPHA2)*REC0H**2-A20H)# 
l ( l.-EXT*COSPT  »  +  ( X12*REC0H+2.#ALPHA( IX )»| 3.-4.+ALPHA2 )*REC0H**2 ) 
2*FACTR2)/(2.»OMEGA) 

XMAXO«XOMEGA I KRO) *FACTR 1+XDOT C  KRO) +FACTR2-Z ( KRO )*FACTR3 
l-SI KRO )*FACTR4-52ND I KRO ) *FACTR  5 
I F ( XMAX-XMAXO  )329. 330*330 

329  XMAX«XMAXO 

330  XMINO-XOMEGA <KRA)*FACTR 1+XDOT <KRA)*FACTR2-Z(KRA)*FACTR3 
1— S ( KRA )*FACTR4-S2ND  <XRA)*FACTR5 

IF (XMIN-XMINO  >335*335.332 
332  XM!N*XMINO 
335  CONTINUE 
N  IN*  1 
GO  TO  214 

- FREQl*ZERO  FOR  FOUNDATION  ACCELERATION - 

350  GO  TOI 352*365 ) »NIN 
352  XOMEGA (1)«X0 
XDOT 1 1 )*XD0T0 
H2«H**2 
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iju  sou  i-£»n 

XOMEGA(I)»XOMEGA( 1-1 )+H*XDOT( 1-1 l-ZII-1 )*H2/2.-S( I-l)*H2/6.+ 
lS2ND(I-l)«H2/24. 

360  XDOTtI )«-Z(I-l)*H-S(  1-1 >*H/2.+S2ND(  I-1)*H/12. 

GO  TO  756 
365  CONTINUE 
NIN-1 
GO  TO  214 

- COMPUTE  X»  XDOT  FOR  OMEGA-ZERO - VELOCITY  INPUT 

400  GO  T0(402»425 ) »NIN 
402  XOMEGAt 1 )-X0 
XDOT ( 1 l-XDOTO 
DO  410  I-2.N 

XOMEGAt I )«X0MEGA(  1-1 )+XDOT( 1-1 )*H-S( 1-1 )»H#.5+S2ND( 1-1 )*H/12. 
410  XDOT ( I )*XD0T (I-l)-S(I-l) 

GO  TO  756 
425  CONTINUE 
NIN-1 
GO  TO  214 

-  SHORT  FORM - 

650  CONTINUE 

IF ( ABSt  FREQ-FREQl )-.lE-04)651 .651 *652 

651  WRITE(8.912)  ALPHA(IX) 

652  GO  TO(658#655 ) » I  ALPHA 

655  XOMAX-AMAX1  t ABSt  XMAX ) »ABS( XMIN) ) 

GO  TO (668 *669 ) t IFREQ 
658  RMAX-SORT ( XOMEGAt  KH ) **2+XDOT ( KH ) **2 ) 

XOMAX-AMAX1 ( ABSt  XMAX ) .ABSt  XMIN ) .ABSt  RMAX ) ) 

GO  T0(681.678) • IFREQ 

668  WRITE(8*906)  FREQ.XOMAX 
GO  TO  215 

669  WRITE (8 *905)  FREQ.XOMAX 
GO  TO  '215 

678  WRITE(8.907)  FREQ.XOMAX. RMAX 
GO  TO  215 

681  WRITE(8.908)  FREQ.XOMAX. RMAX 
GO  TO  215 

- -LONG  FORM - — 

700  CONTINUE 

IF(ABS(FREQ-FREQl)-.lE-04) 701 .701.702 

701  WRITE(6.915»  ALPHA(IX) 

702  GO  TO (704. 7051. I ALPHA 

704  RMAX-SQRT t X OMEGA ( KH ) »*2+XDOT ( KH ) **2 ) 

GO  70(743.742). IFREQ 

705  GO  TO(726. 727). IFREQ 

726  WR!TE(6.916)  FREQ. XMAX. XMIN 
GO  TO  785 

727  WRITE(6.910>  FREQ. XMAX. XMIN 
GO  TO  785 

742  WRITE(6.910)  FREQ. XMAX. XMIN. RMAX 
GO  TO  785 

743  WRITE(6«911)  FREQ. XMAX. XMIN. RMAX 
GO  TO  785 

- PRELIMINARY  CALCULATIONS  FOR  MIN.  AND  MAX.  RESPONSE - 

756  GO  TO  (758.757) .IALPHA 

757  NX-N 
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1.0  IU 

758  NX=KH 

759  XMAX*X0MEGA(1) 

KMX*  1 

XMIN*X0MEGA(1 ) 

KMN*1 

00  760  1*2. NX 

I F { XMAX— XOMEGA ( 1 1)761.762.762 

761  XMAX*XOMEGA( I ) 

KMX*  I 

762  I F ( XMI N-XOMEGA( I )) 760.765.765 
765  XMlN*XOMEGA ( I ) 

KMN*  I 

760  CONTINUE 

I F ( KMX-1 ) 771 .770.771 

770  KMX*KMX+1 

771  IF(KMX-NX)773. 772.773 

772  KMX«KMX-1 

773  IF (KMN-1 >775.774.775 

774  KMN*KMN+1 

775  I F ( KMN-NX J777. 776.777 

776  KMN*KMN-1 

777  CONTINUE 

IF(XDOT( KMX-1 )«XDOT (KMX) >778.778.779 

778  KOOK*-l 
GO  TO  781 

779  KOOK*0 
781  CONTINUE 

I F ( XOOT ( KMN-1 ) *XDOT ( KMN ) >783.783.784 

783  KICK*-1 
GO  TO  950 

784  KICK*0 
950  NIN*2 

GO  TO  207 

- EITHER  RESTART  WITH  NEW  FREQ.  OR  RESTART  WITH  NEW  ALPHA  OR  STOP 

785  CONTINUE 

IF(ABS(FREQ-FREQ2)-.lE-03) 800.755.755 
755  FREQ*FREQ+OELTAF 
GO  TO  207 

800  CONTINUE 
GO  TO  190 

801  .WRITE!  6.909) 

STOP 

900  FORMAT! 1H1.49X.22HSH0CK  SPECTRUM  PROGRAM//50X.21HAPPLIED  MATH  LAB 
1DTMB.///) 

901  FORMAT (40X.38HF0UN0AT ION  ACCELERATION  VS  TIME! INPUT ) //18X .23HF0UND 
1ATI0N  ACCELERATION.36X. 15HTIME  IN  SECONDS) 

902  FORMAT ( 42X.35HF0UNDATI0N  VELOCITY  VS  TIME  ( INPUT) //20X.19HF0UNDATI 
ION  VELOCITY. 35X.15HTIME  IN  SECONDS) 

903  FORMAT! 1X.F35. 9.  44X.F20.9) 

905  FORMAT! 15X.F10.5.29X.F12.7 ) 

906  FORMAT! 10X.11HFIRST  FREQ-.F10.5.7X.13HMAX  RESPONSE* »F12. 7/ ) 

907  FORMAT (15X.F10.5»2(29X»F12.7) ) 

908  FORMAT! 10X.11HFIRST  FREQ«»F10.5»7X.13HMAX  RESPONSE* »F12. 7. 12X. 
117HMAX  RES  RESPONSE*. F12. 7 ) 

909  FORMAT! ///38X.25HPROGRAM  RAN  TO  COMPLETION) 

910  FORMAT! 10X»F10.5»3(18X»F12.7) ) 

911  FORMAT! 5X.11HFIRST  FREQ«.F10.5»9X»9HMAX  RESP«.F12.7.10X»9HMIN  RESP 
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y  X  £  runn*  I  \  1  ni9?7A»  iunsnvn  i  i  wnnr  r  »c  w  .  —  - 

1 /9X# 22HFREQUENCY ( CYCLES/SEC • ) ♦  19X  *20HMAX  RESPONSE  X  OMEGA* 18X.26HM 
2AX  RESIO  RESPONSE  X  OMEGA) 

915  FORMAT ( 1H1 *55X*9HL0NG  FORM//45X*20HDAMPING  COEFFI C I ENT= . F 10.8//// 
17X • 16HFREQ) CYCLES/SEC ) » 12X *20HMAX  RESPONSE  X  OMEGA* 10X .20HMIN  RESP 
20NSE  X  OMEGA *8X.22HMAX  RESID  RESP  X  OMEGA) 

916  FORMAT (4X.11HFIRST  FREQ= *F 10. 5 ♦ 10X . 9HMAX  RESP= » F12. 7* 10X *9HMI N  RES 
1P=*F12.7) 

917  FORMAT <1H1*18X»1HI»38X*4HS( I ) »34X*7HS2ND( I)// 
1(17X*I4*34X*F12»7*28X»F12.7) ) 

918  FORMAT ( 1H1 *50X»19HINTERMEDI ATE  VALUES/19X.1HI *35X*9HXOMEGA( I ) * 
131X*7HXOOT( I )// 

2(17X*I3,34X*F12.7,28X*F12.7) ) 

920  FORMAT ( 3F20.8 ) 

921  FORMAT ( 6F15.6 ) 

END 
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